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ABSTRACT 

The method developed by Osterberg, Sarkar and Kruck 
fOr Obtaining the concentrations of free ions in solutions of 
complex equilibria from oH titration data has been studied, 
to establish the conditions under which the method gives 
accurate values for free ligand and free metal concentrations 
in Systems containing a variety of complexes, including 
protonated, hydroxy, mixed ligand and polynuclear complexes. 
Simulated titration data has been used so that the true values 
of free ligand and free metal concentrations would be known. 
Calculation procedures are described for each step in the data 
€valuation, including procedures for extracting information 
about the stoichiometry of the complexes from the unique 
information provided by this method, The effect of various 
systematic and random errors is also considered. 

An automated gravimetric titration system was developed 
for the purpose of collecting highly accurate potentiometric 
data suitable for use with the Sarkar-Kruck method. In 
order that none of the theoretical requirements of the method 
be violated, the system was also designed to maintain total 
ligand and metal concentrations constant employing computer 
controlled addition of concentrated reagents from an auxiliary 
DULeL., 

The Sarkar-Kruck method and the gravimetric titration 
apparatus were tested in a study of the complexation of zinc 
by aspartic acid. A system having simple complexation 


chemistry was chosen in order to avoid problems in the 
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assignment of a model to the metal ligand reactions. Many 

of the calculation procedures developed on the basis of 
simulated titration data were tested and found to be valid 
With experimentals data. Values for jthe acid’ dissociation 
constants of aspartic acid and the formation constants of 
complexes between zinc and aspartic acid agree well with those 
reported in the literature. 

The complexation of zinc, cadmium and’ lead by) glutathione 
was studied using the method and apparatus described above. 
Pontraryato: somesceports Gin tthe; duterature. itjwas*ifound chat 
at low pH the metals bind to glutathione only through the 
sulfhydryl group with minor complexation occurring through 
the glycine carboxyl group. The glutamyl amino group remains 
maotonated ..cinvaddition, there wassevidencesfor poly-= 
merization reactions occurring through the glycine carboxyl 
at low pH, and through the glutamyl terminal groups at 
somewhat higher pH. These reactions have been suggested by 
spectroscopic data but other workers using potentiometry have 


avoided consideration of these complicating processes. 
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CHAPTER 1 


INTRODUCTION 


The chemistry of complexation reactions has been the 
subject of considerable research in the past because of their 
importance in analytical methods (1), and more recently 
because of interest in the role which metal ions play in 
biochemical systems. The methods which have been used to 
determine formation constants include a variety of spectro- 
scopic methods, calorimetry, reaction kinetics, liquid- 
liguid partition, solubility measurements, On exchange, 
potentiometry, polarography, conductivity, light scattering, 
osmosis, electrophoresis and many more (2,3,4). Of these, 
potentiometric titrations with a pH electrode have probably 
been the most frequently used, and are generally considered 
to give good results. The reasons for this popularity 
include: (1) the applicability of the pH titration method 
to a wide variety of complexes, (2) the relative simplicity 
of the experiment and the reliability of the data, and (3) 
the straightforward relationship between the experimentally 
measured quantities and the chemistry in the solution. A 
carefully designed and executed experiment will provide a 
large amount of data, collected under controlled conditions: 
the concentrations of metals and ligands in solution and 
the quantity of titrant are accurately known, and the ionic 


strength and temperature are held constant. The major 
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limitation of the method is that it does not provide any 
information about the complexes at the molecular level. 

This means that a model of complexation must be derived 
intuitively, or by trial and error, or must rely on spectro- 
scopic or other data which gives Ee eens ape information 
at the molecular level. 

The focus of this thesis research has been to determine 
ene formation, constants for the zinc, cadmium and lead 
complexes of the tripeptide, glutathione. Glutathione is 
an important component in the red blood cell and the 
formation constants of its complexes are of interest for 
the Ange esbeicd ars ofthe: distribution and transport of 
toxic metals in the body. These complexes have been 
studied by several investigators, and formation constants 
reported (5-8). However, there has been no agreement about 
the model of complexation. 

A pH titration method recently proposed by Sarkar and 
Kruck (9) has the potential for providing information not 
directly available by any other pH titration method, 
namely, it derives the concentrations of the free ligand and 
the free metal from the experimental data without any pre- 
conceived model for complexation. The work reported in 
the wtirst partverethis thesis deals with ay study and 
further development of this method, particularly with 
respect to the design of an experiment which will give the 
most reliable data. The results of this study indicate 


that’ large volumes of highlyaccurate*data aremneeded’ for 
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the best application of the method. For this reason an 
automatic titration system based on a Sartorius electronic 
balance’ anda PDP-11/10 minicomputer was designed and built. 
With this apparatus and the data handling routines developed 
for the Sarkar-Kruck method, the complexation chemistry of 


glutathione was studied. 


A. The Evaluation of Complex Formation Constants from 


Potentiometric Titration Data. 


The pH titration is a convenient technigue with which 
to collect data reflecting the competition between protons 
and metal ions for complexation sites on the ligand 
throughout a wide range of conditions (2,10). Titrations 
are performed with sample solutions containing several 
combinations of total metal and total ligand concentrations, 
and through a wide range of pH. The experimental data is 
related to the formation constants of the complexes by the 


mass balance expressions 


a = TH] - [087] +) Boar P1200 (1) 
amie +) gQeqe Pra F007 (2) 
CenamiMl e+ » Boge Pra Fn? (3) 


where Ci is the total concentration of titratable proton and 
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is calculated from the amounts of acids and bases added to 
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the -solution,, and C C,, are the total ligand and metal 


Le 


concentrations. The formation constant for the complex 


hes 1s Ble. where 
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ism = De em (4) 
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L! Cu and the hydrogen ion concentration 


(from the pH) are known at each data point while the 


The values of Cue C 


concentrations of ligand and metal, [L] and [M], and the 
formation constants are unknown. The evaluation of the 
constants from a set of experimental data may be very simple 
if only one or two complexes form, but at the other extreme 
may be very complex if protonated and polynuclear species 
exist'in solutron. 

A variety of approaches for the handling of the simpler 
Situations have been reported (11,12). These generally 
are only applicable to very select situations, and very 
often require additional measurements, for example with an 
ion selective electrode. The most widely known and 
applied is an approach described by Bjerrum (13). The n 
method is’ based upon the calculation of the average number 
of ligands complexed per metal directly from the experimental 
data. The method fails, however, if protonated, hydroxy 
Or polynuclear complexes form. 

Methods which handle more complex situations fall into 
two categories, those based on linear statistics and those 


based on non-linear statistics (2,14). If the data can be 
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arranged into a set of equations such as 
N 
y =)) x, (5) 


where x are measured quantities, and es are a set of N 
unknowns, the constants may be found by a linear least 


squares calculation to minimize the residual U. 
xs 2 
Liss hile Daye (6) 


This method is not generally applicable, and usually 
requires several extra measurements, aside from pH, to 
provide suitable values of xX. 

Non-linear methods are very complex, require consider- 
able computer time, and in addition, are reportedly not 
very reliable (14). The technique, in general, involves 
systematic adjustment of the formation constants until 
they provide the best fit to the data. A sequence of 
calculations will involve: 

(a) estimation of the formation constants 

(b) with these constants, solve equations (2) and 

(3)for thesunknowns;,. [hij -andeiMic «(it must be 
noted at this point that ‘the, valuesvort? fu), and 

[M] are dependent on the model for the complexes 
and the estimated constants for those species, and 


may be quite different than the [L] and [M] 
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(c) © calculatefithe | expected’ value ‘oft an experimentally 
measured quantity, for example Cin GAUSS) (15/16), 

(d) adjust the values of the constants in the 
direction, and with the proper magnitude to reduce 
EhetreswdualsUse (boqnces/) ah Returntois tepe (cb) 


U = 5 (coXP i ccale, 2 (7) 
He H 


antidmohe residual. Winisemnimaimiged: 


Steps (b) and (d) may be executed in several ways, and 
these differences have led to the major data processing 
computer programs. For example, two approaches have been 
used to adjust the constants, step (d). The first is 


based on the first derivative 


Ove 


dQ. = residual (8) 
O2n 


OY LS ae Dveuaets 


and is used in GAUSS (15)., SCOGS “CI7) > “and: (MINTOUAD 18)" 
The second is the so-called pit-mapping approach of programs 
LETAGROP (19)" and LETAGROP VRID "(20)" and 1s based"on a 
Telacion sami lal tO BON, oO, DOUt i). Cemis -O.1, Ele 
fecond derivative. in addition, programs have been 
written which use both approaches, LEAST (21), and others 
using a completely different method, STEW (22). 

Recently a pH titration method was reported by Sarkar 


and Kruck (9) which is capable of providing the 
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concentrations of free ligand and free metal throughout the 
pH range, using only the experimental data from an 
appropriately designed set of titrations. This is in clear 
contrast to the methods described above, where [L] and [M] 
must be calculated using a hypothetical model, and guesses 
Porsche :Otmatloneyconstants.; With. the values. of. [Db] “and 
[M] derived experimentally the non-linear problem of 
solving for formation constants reduces to a linear one, 
even for systems which include hydroxy and polynuclear 
complexes. The formation constants may be extracted from 
Equations 1, 2 or 3 using a linear least squares calculation. 
Application of this method to a number of metal-amino acid 
systems has been reported (23-26). Nonetheless it appears 
that the method has not been proven applicable to poly- 
nuclear complexes, and that the exact consequences of error 
originating in the experimental data, and in the data 
handling steps have not been described. In addition, the 
reports published by Sarkar and Kruck did not describe the 
experimental method in detail. From these reports it 
appears that no attempt was made to determine the conditions 
which would give the best experimental results, and as a 
result it seems that the experimental methods used in these 
reports may have violated the theoretical requirements of 
the method. 

The main thrust of the work reported here was to 
simplify the calculations as muchas possible, to prove 


that the method is applicable to the more exotic complexes, 
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and finally to establish the experimental conditions which 
will give the most reliable data possible. The results of 


this study are described in Chapter II. 


B. Automatic Titration Apparatus. 


Of the many automatic solution handling methods 
available, the most popular are those based on a stepping 
motor driven piston buret. Several units are available 
commercially, including the Mettler DV1l and the Radiometer 
ABU13, and many workers have devised their own (27-34), for 
example using a commercially available micrometer driven 
buret such as the Gilmont ultra precision buret, powered 
by a stepping motor. These types of systems are very easy 
to automate; a delivery of a known volume only requires a 
series of electrical pulses. They are however, high in 
BEILCG, and usually low; in precisyvon. (iis) lack of saccuracy 
usually stems from imperfection in the mechanical drive 
from the stepping motor to the piston, and the usual 
problems with volumetric methods, including temperature 
erlects,, Gatiol calibration dlassware, leaks, and so on. 

Only a few alternatives to the piston buret have been 
reported. Coulometric generation of titrant is an old, and 
potentially very precise method(35). Its general application 
to high precision acid-base titrations has been limited by 
two factors; first, the titrant generation reaction must be 


very efficient and allow high current densities and, second, 


Ef th 
‘By ® wii eon: “he at 
ri 
I ;' 
; } t 
be “ 
4 
Whe 
a 6 
‘ “ he. 
" Pr 
F dee , 
6h r * i th ib 
ie Fs 
uf 
ein harem tie) ar ee Te eit OY Ge 
j 
ea ae eee Ps } ci, te H 
LAVAL OO) it es Sy nk = ae PS i 
rn 9, Il i i 
y bal i ’ 
nA ? Ad . u ate yh as. Ba ia ba et it Ce q 
RSH: FRE OER alka. Wh 
: na ; e . 
« - ” | 
ee hy Are fi POMEL, os 
Pi ' 
y * , 
‘ aaa La 55 its it - i 4 i 
; re 
i) | Sad ; v ‘ ‘ i 
Vo paeu ytcmyd tore ties 
¥ ph w ied , i 
Heo Vee ee: 4 a ks ae ‘ 
[ ¥ 
' t ‘ 
ee na i if da Hi eit Lene 3S 
Q a x wuge 1 
ian Pe ” : 
i. fi ire Lal. see pone £0 bh 
ay Sian r i "edi penne, : A C 
ai : a “ Oo” wl ne | @] ; A a 
baat o ? & ne Sak a ey SUR ey eR ae 
' Lan) yy y i ay : . 
; f De ce wii u -, : ter a ¥ 
“awd fas" eat flit: FE MIT RT 
Rie OBee Ai ‘ uae ¥ 


Paw? ry wa rs aA a 7 any 


eT RL AE TO y ea ya 
i bm \ @ ¥ fs f Rit id 4 Whee, 


or ; ‘ 
ey an | ie fan) On 
V 


aw, i 0 abies Sans 5 nt 


Pha’ f Bs arr A, 
8 kes Bs eis ine at Wet! i = ; 
un re ao baer: date : 


(Los, med Aaron ee: ua 


ni oe 4 a ri : 2 

TRUM Hcp ee ae ae | I 
i *) b - | Pat “Ang: 

ri y ‘ i: phe oe fa i j 


mie 4 ae 1h y ‘ tie, . Lb oy Ay 


er RL Sa pee 0 > Ba SRE, Maree a ta Pier ig 3 fh 


= 


See: 


if 
ie 


oh 
abi 


f a hy J uy Fac) 


aaa oh), ted cnet Te aa riiye Te ne a ea Fae: 


Srey ae iia ihe ’ 7? 
asi: oer eae fond ei ag bios ,. 


re 
i yi ¥ or, ' 
cl 


= iv 


ee i ice n ut M6 


jae? had 
i, 


AN i i yea i es ; A : 
ey “¢ Ly 
ee stp aa Fe. 


‘dian Bt te fab a6, 68 ee at?) Mote 5 heer PaaS a aaa 


» sae 
ayn aoa oh none oF . Per Higa me tow a asanioo 


PF Saray | erhdtgote ie tae NE 


Neeson 5 ay 5%, 


ee de yf i” A ou 
a 


antl ee sete rb ak ccna a ona er te ibe 


iv . on 
f oan 
i a ny: 


the working electrode must be isolated from the sample 
solution if the sample contains electroactive species (36). 
Coulometric methods can be easily automated’ (37). 

A unique form of real time titrator has been reported 
(38), in which the sample and an added colour indicator is 
passed through a narrow dialysis tube within a bath of 
titrant. The titrant diffuses into the tube and the point 
where the sample is neutralized is marked by the colour 
change of the indicator. 

A titrator based on the delivery of a series of sub- 
mvcrolitre droplets (39,40) was. reported by Hieftje (41,42), In 
this system, titrant was forced through a capillary which 
was vibrated with a piezoelectric crystal. This dispersed 
the liquid into individual droplets, which the authors 
ingeniously counted, and directed into a sample solution 
using electrostatic forces. Another variation was reported 
(43) in which the time during which the droplets were 
allowed into the sample was measured instead of the number. 
Both of these systems suffered from very poor precision. 

The handling of solutions on the basis of weight has 
long been regarded as superior to any type of volumetric 
technique. However, due to the slowness of the weighing 
Operation, and the virtual impossibility of automation, 
gravimetric procedures have not gained much popularity. 

The advantages of gravimetric titrations over volumetric 


methods include: 


(a) elimination of drainage, reading, and temperature 
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change error inherent to volumetric equipment, 


10 


(b) elimination of time consuming gravimetric calibration 


of volumetric equipment, 
(c) ability to measure and dispense accurately known 
weights of volatile or highly viscous liquids, 
(d) ability to measure with high precision small 
quantities of solution for microchemical work, or 


with costly or scarce materials. 


Manual methods have been described for special situations 


(44,45), including work with non-aqueous solvents, and 
where particularly good precision was desired (46). 
Malmstadt (47) has reported an automatic system in which 
the weight of reagent delivered is determined by weighing 
the container into which the material is dispensed. This 
method is not suitable for a titration in the usual manner 
because no provision is made, or seems possible, for 
monitoring a sample property, for example the pH. 

In this thesis a completely automated gravimetric 
titration system is described, based on a Sartorius 
electronic balance and a PDP-11/10 mini computer. The 
delivery of reagent from a reservoir located on the 
balance pan was initiated by an optical device. After 
delivery, the weight of titrant delivered was known by 
the decrease in weight of the unit on the balance pan. 
This system combined the accuracy and precision of a 


technique based on a weighing procedure, with the 


speed and reliability of automation. Gravimetric titrations 
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were run with high precision, and in a fraction of the 


time needed for equivalent manual Sperations, This system 


will be described.in considerable detail in Chapter IV. 


C. The Complexation Chemistry of Glutathione. 


Glutathione is a tripeptide composed of the amino 


acids glutamic acid, cysteine, and glycine. 


S 2 ° " 
HO—C—CH—CH —CH—C-NH-CH—8—NH-CH,—C—On 
Bs 1 
NH : CH, 
Ysa 
glu cys gly 


The binding to metal ions may occur through the carboxyl 
and amino groups of the glutamyl terminal, through the 
Suifhydry!l group, and: through ithe glycine carboxy! group. 
In addition, it has been suggested that some binding may 
occur through protonated and deprotonated peptide linkages 
(5,6). With this number of binding sites, and the number 
of groups which are protonated at various pH ranges, the 
complexation chemistry of glutathione is not simple. 

The earliest attempts to unravel the complexation 
of glutathione (8,48-50) relied upon relatively simple pH 
ered experiments, and on polarography, but the data 
handling methods were not sufficiently advanced to handle 


protonated and polynuclear complexes. The development of 
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Spectroscopic techniques such as NMR allowed the study of 
metal-ligand interaction at the molecular level, for 
example indicating which groups were bound to the metal. 
The sulfhydryl group is known to be the most active 
in the complexation of glutathione to a number of metals, 
binding very strongly to "soft" metals such as mercury, lead 
and cadmium and less strongly to zinc, and other transition 
elements: (7,50,51,52)... The interaction: of mercury with 
Bie sulthydryl rs particularly strong, and the metal ion 
can only be displaced by protons at very low pH. Since 
the carboxyl and amino groups do not bind as strongly, 
several protonated complexes exist at low pH (51). Other 
mercury complexes, containing as many as three metal ions 
and two glutathione molecules have been shown to exist (50, 
Sls). 

Studies of the binding of glutathione to other metals, 
Miecluding cu,;.Ni, berand, Pd.(>55-55) “have undicated thatthe 
sulfur is again the primary binding site. Complexation 
does occur through the glutamyl and glycine residues as 
well, but it has not been possible to determine if those 


groups are bound to the same metal ion as the sulfhydryl, 


or to another. It has been pointed out that the complexation 


of the glutamyl and sulfhydryl groups to the same metal 
would form an unstable ten membered ring, indicating that 
at least in some cases, a polymerization process is likely 


Eo voccur. (7). 


The development of non-linear methods of treating pH 
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titration data made possible a more detailed analysis of 
the composition of a solution containing a mixture of 

metal ions and glutathione (5,6). The derivation of a model 
of complexation has, however, been largely a matter of 
intuition and guesswork, and the results of studies 
involving zinc, cadmium and lead do not agree with each 
other nor do they agree with the results of NMR experiments 
(7). Since pH titration methods do not, in general, 

provide any information about the types of complexes 
existing in solution, the choice of a model must be guided 
by some other evidence. The Sarkar-Kruck method provides, 
unlike other methods, experimental values for [L] and [M] 
which, as is discussed in Chapter II, allows some rational 
deduction of a model of complexation. In addition, the 
number of models that will fit the data is considerably 
reduced. For these reasons it was hoped that the Sarkar- 
Kruck method would produce more meaningful results about 
the complexation of glutathione than have been thus far 
reported. The results of the metal-glutathione experiments 
and the discussion of models for complexation appear in 


Chapter VI. 
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CHAPTER II 


A STUDY OF THE OSTERBERG-SARKAR-KRUCK METHOD OF EVALUATING 
FREE METAL AND FREE LIGAND CONCENTRATIONS IN SOLUTIONS OF 
COMPLEX EQUILIBRIA 


A. Introduction. 


Sarkar and Kruck recently described a pH titration 
method in which the free metal ion and the free ligand can 
be calculated simultaneously throughout the central titra- 
tion of an appropriately designed set of titrations. This 
method is appealing in that the concentrations are obtained 
directly from the experimental] data, without assuming any 
hypothesized collection of species for that system. The 
calculation of stability constants then reduces to a linear 
least squares problem using any of the mass balance 
equations. If the stoichiometries of the complexes are not 
known, it still is necessary to postulate a model. It is 
to be expected, however, that fewer models will fit the 
data because the free ligand and the free metal concentra- 
tions are known and are not derived from the model. 

In this chapter a critical examination of this 
approach is described. The emphasis of this examination 
has been to establish the conditions under which the 
method gives accurate values for the free metal and the 


free ligand concentrations in systems containing a variety 
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of complexes, including protonated, hydroxy, polynuclear and 
mixed ligand complexes. The approach has been to analyze 
Simulated titration data, so that the: true values of free 
metal and free ligand concentration would be known. In 
addition, since the method provides information not 
available in the study of complex systems by other data 
evaluation procedures, ways in which this information can 
be used to provide guidance in the selection of a model are 
considered. 

This chapter refers exclusively to simulated titration 
data based on several hypothetical complexation systems. 
The complexation systems and the parameters used in the data 
Simulation are summarized in Tables 1-5. Series A and B 
(Table i) are.-for: a poltyprotic acid, Series 'C (Table, 2) is 


£Or a triprotic Jigandyand the complexes HML, HML. and) ML, 


2 
and Series D (Table 3) is for a diprotic ligand and the 


Conplexes Ml, ML, Mijgand HM. . Serves Em (Table 4) 


2% 5 1 

includes polynuclear species and Series F (Table 5) includes 
mixed ligand complexes. In the following discussions, 
examples of calculation methods and the effects of error 
will be based on these series of titrations. The computer 


program used to simulate the titrations of Series A through 


F is shown in Appendix A. 


Overview of the Method 


This method provides the free concentrations of ligand 
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Table 1. Simulation Parameters: Series A and B 


Complexes Formation Constants 8B 
HL ita) Steals che 
H3L OS rhe 
HL 10 eam Oue 
HL Toten Oe 
Jvesreauly Intepal pH Titrant 
Ligand pH Interval Concentration, 
Eitration Cy, M 
A, Bee ie 2 ee s00 0.05 1.00 
A, DeNhe TO". 3.00 0.05 1.00 
A, ue Or 3.00 0.05 1.00 
A, fae Toe | HOO 0.05 1.00 
As WoIs dps) We S00 0.05 1.00 
B, aeosdipe Sea 0.075 0.1800 
By SOM Om fees 206 0.066 O. 2013 
B, Olea ore) | 3.399 0.071 0.220 
B, iedeaiGe 4 tacnat 0.056 olay 
B hase Ghar = Ayes 2a 0.068 0.2500 


Table 2. Simulation Parameters: 


Series?®c 


Complexes Formation Constants 
H3L 10 x) 10" 
H5L USO re 3 iif 
HL ibaa Saye 
HML 0a) 10m 
HML. lity oe 
ML» 1.0 x 10°8 
DH Ranger 3.0) — L020 
PHeiintenval ss 10.1) (7 iudaran points) 
initial Volume: 300.m1 
Titrant Concentration: ‘no titrant 
Middle Titration: Cc, = SOEs 10a. Cr = 
Titration Cy, 
oA 2.99 x 10°? 1.00 
e. Se 0paxe 108 1.00 
e. 35 Doe ean 1.00 
Si, Said os eae 9.9 
C. 3.00 x 10> 1 Ou 


Cc' Series: With Dilution 


As above except 


Titrant Concentration: 0.4 M NaOH 


x 


x 


x 


Ly, 


| 
| “giod 
A 
Cede bod akiet, ea eae Re is ere 
tae Bhi ¢ a nih i. Lite 
iain fe betel WA tra ane He. ie if, tonsa) $i 
i ih 8 Wee 3 ot al oh: VE ROO. a Pe oe : ive eri £ iy « 
y } ; A pe 
; vn Ae) fe nomi ras Aumann y é ae o 


Tats +00 . i) Pe Ty Abe a ie mi 


Le net Ai a th ee | 
: Le Hee (9! . bi acl | 


hihi: a Pane a ACF Oka pte 
Te oh) e 


Table ws 


Complexes 


HL 


HL 


ML 


pH Range: 


pH Interval: 


Initial Volume: 


SUE —peLU sO 


5.00 sm 


Titrant Concentration: 


Middle Titration: 


Pitration 


Cy ea 


no titrant 
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Simulation Parameters: 


x 


x 


Series D 


0. (71 data wpoants) 


Formation Constants 


10 


TO! oe el 


90) X16 


9.5 x*104 
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Table 4 . Simulation Parameters: Series E 


Complexes Formation Constants 
Hal  0ae 
HL 1.0 x 1026 
10 
HL Le Ser) 
HML 1 ax 102” 
HM, | 1.0 x 1028 
Ld 
ML, is On se LO 


pe eh aera hove (oly Wt OE iE ehesa 6 
DHe intervals 0.1 (71 data moints) 
Initial Volume: 300 ml 


Titrane Concentratrone no, CLerant 


Middle Titration: C, = 2.5 x 1077; ¢, = 1.0 x 10° 

Titration om ce 
E, S.0 x 107" 1.0 x 107" 
E. x 10 1.0 x 107° 
E Aux eee 1 OSX noe 
E, 5 2110 1.0 x 10 
E. 2.6 x 107 1.0 x he 
E¢ 8 x 107 1.0 x 1077 
E. .0 x 10 1.0 x 10 
FE ho T= S.0 Ronde 
E. .5 xX 10" 9.0 x 107" 
Es 2.5 x 10" 9.5 x 10 if 
E so. B0 tes 0 Saescae 0) 

11 3 ye 
E12 Pangea ae 4 a oboe hee 
E 2.5 x 10 1.20 x 10 
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Table 5). 


Complexes 


HL 
HoL 


HL 


ML ' 


HMLL' 
MLL' 


pH Range: 


Du cincervial.: 


Initial Volume: 
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Simulation Parameters: Series F 


Formation. (Constants 


Pitrants Concentration: no Eltvantr 


Middle Titration: 


Preratcion 
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iro Soe? 
10% APO” 
0. 10? 
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i 10° 
BeOS OL 
108 oe 
ie 1014 
3. 0— 10.0 
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sopemt 
7 = 4 = 
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and metal in solutions of complex equilibria, and for 
Simplicity in this discussion will be called ETCGS;* £Orer ree 
fon Concentrations in Solution. The basic relationships 


in this method are 


pH 
pL] = pIL], - {ac dpi (9) 
dc 
PHI, L Cyr PH 
PH) ac 
p[M] = PIM]. - H. dpH (10) 
pu! Cy C, PH 


where Cy is the concentration of titratable proton, C. and Cy 


L 1 


are the total analytical concentrations of ligand and metal 
respectively, and [L] and [M] are their free concentrations. 
Alternate derivations of Equations 9 and 10 may be found in 
references 9 and 56. An abbreviated example may be eave 


in Appendix B. 
For systems containing more components, equations 


analogousmto Sand <L0Mcan) be written’ for each component 
in the solution. For simplicity let us only consider the 
component L, remembering that the discussion is applicable 
to every component in the system. 

The quantity PIL], is the negative logarithm of the 
concentration of free L at some initial pH,. A pH, is 


usually chosen at which no complex is formed, and 
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28 
consequently [L] is equal to (or may be calculated from) 
C 


The integral of (dc. /dc,) is evaluated from PH, 


1 ae mr PH 
BO; Che pH iats~which pil] is. to’ be evaluated. 
One experimental meth pe ini 

p ethod for obtaining Sa Minae Kober: 
as a function of pH is with a series of titrations where the 
EOLals Concentration. of Lis’ varied,Awhile all athe, other 
components are at constant concentration. At a given pH, 
Cu is}extracted. from.ecach,of,the, titrations: « The. derivative 


SS he is the change in Cy with Che measured at C 


,PH L 


of the middle titration. In exactly this way the derivative 
at each pH throughout the titration is calculated. The 


values of (dc, /dc,) o are integrated from PH, to a given 


m’PH 
PH, ‘and.p[L| <atithetipH is) then, calculated from) Equations9. 
The experiment is deSigned in such a way that, given a 
system of composition Ch and Cur referred to as the 
composition of the middle titration, each component is in 
turn selectively varied. The FICS method will then provide 
fuleand (M].dsherfunctionvot spl, sboreawsolution,ofMehat 
Overall composition. 

The.stabi laity «constants, of,the complexess are scalculated 
bVeawdinear) Least squanestit of.the-JL)eands[Mi,deataginto 
the mass balance equations. Any one mass balance is com- 
pletely sufficient unless a given complex does not contain 
any of that particular component. 


In the following section, the calculation procedures 


at each step in the method are considered. 
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B. Calculation Procedures in the FICS Method. 


Cy Calculations 


The extraction of useful information from the pH 
titration requires accurate knowledge of several parameters. 
The FICS method requires the total ligand concentration, 


Cc the total metal concentration, C the total concen- 


re, 


tration of mineral acid added, C 


M’ 
at the pH and the concen- 


tration of added base, C All of these quantities change 


B° 
during a conventional titration and careful work is needed 
to provide accurate knowledge of their values. 


The first quantity extracted from the raw data is the 


total concentration of titratable proton, Cue 


x Cy) He ect (ids) 


where Nop is the number of disposable protons on the ligand. 


In order that the FICS method be used, Ch must be known at 


closely spaced pH values throughout the titrations, and 
moreover, it must be known at the same pH's in each of a 
series of titrations. Since this is difficult experimental- 
ly, some manipulation of the original data is necessary to 
find these C.. values. Several possibilities exist, 


H 


including the fitting of the pH-C,, data to a mathematical 


function, perhaps a type of polynomial, and in subsequent 


steps extracting Ch at any needed pH. This is the 
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approach taken by Sarkar and Kruck (9). However, in our 
hands a program of this type obtained from other workers 

did not function very well. It was found that linear 
interpolation between existing data points was sufficient 

to provide Cu aCispecific pH" si throughout the itittrakion. 
Efficient use of experimental data means fixing a pH interval 
approximately equal to that of the available data. 

The linear interpolation is limited in two ways; first 
it does not use the entire body of data as a unit, and 
therefore does not have the smoothing potential which other 
procedures may have, and second it will bias the data in 
some cases. A bias to large C,, will result if the Ch vs 


H 


pH curve is concave upward, and to low C's if the curve 
is concave downward. In most cases this will give very 
small errors) at, pH intervals of about 0.1 units. 

The effect of linear interpolation on the final results 
Of FICS is illustrated in*Tabbles6. Series A: uses C's 
directly from the simulation program, and Series B uses 
C,'s calculated by interpolation on data generated at a 
varying set of pH's. Possible error was encouraged in the 
second set by increasing the pH intervals from 0.05 to 0.1 


and by increasing the dilution by titrant (lowered NaOH 


concentration). The simulation parameters for Series A 


and B appear in Table 1. 
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Table 6. 


simulated Titration’ Data, Series A and B 
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The Derivative Mast rules: 


The second calculation in the FICS procedure is the 


evaluation of the derivative (de/dce) The quantity 


71 PH” 


Cy has been discussed. Cy will be used to’ denote the total 


analytical concentration which is varied. This may be Che 


the ligand concentration, or C etc. The term Cy 


mM! 
represents all the other total concentrations which are 
fixed. 

The experimental determination of this derivative is 
the most challenging step of the FICS procedure. The 
straightforward experiment is to perform several titrations 


varying one C perhaps Che while maintaining Cu and any 


baud 
other concentrations constant. The simulation parameters 
for Series D, shown in Table 3 illustrates this more clearly. 
At some pH, for example 4.0, Cy is calculated from each 


er the titrativons. This qives,the variation 41n Cy with 
CL She geal coal OPen obs onto cas 4.0° 
The accuracy of this quantity depends on several 


factors: 


(1) the accuracy of Cy and Che 


(De theucurvature OL) the C, vs pH function, 

(3) the number of titrations, and therefore the 
number of C's available, 

(4) the degree to which Cy (the other components in 


solution) remains constant. 
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The function Cy Vor and the effects of problems (1) to 


es 
(4) must be clearly understood in order to obtain the best 


results from the FICS method. 


The variation of Ch with Ch for Series D, which 


consists of a diprotic ligand and the complexes ML, ML., 
ML. and H_|M, Lstshown vin Figure], = The: total imetal 
concentration was-Lixed at 1.0 x Tue y and curves drawn at 
several pH values. The simulation parameters for Series D 
May be found in Table 3. 

Several characteristics of the Cr vs Cr functioniican 


be seen in Figure 1, including severe curvature in some 


regions, and reasonable linearity at excesses of C 
3 


L (Cy was 


i O-x 10 2. antadawtivton it. shouldsatso bel neted that.c 


H 
may vary Over a wide range of concentrations, and tends to 
be small at low ligand concentrations. Several practical 
considerations can be based on these observations. Titra- 
tion conditions must be chosen to avoid areas of curvature. 


Each point on the C,, vs Ch plot represents a data point 


H 
from a separate titration, and the number of titrations may 
be limited. A reasonable estimate of AC, /4C, requires 
therefore that Cy vary linearly with Che and: >, un®aaddition, 
there is some advantage to having ac, of a reasonable 
Magnitude. If the change in Cu with Ch iss linear, only 

two titrations are really necessary, and more will improve 
the situation. Reasonable random error is acceptable since 


allinear least squares fit will tend: to give better 


derivatives than the limit of accuracy of Cu might 
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Figure l. Cy vs C 


2.0 3.0 
C, x10 (M) 


L for a system containing a diprotic 
ligand and complexes ML, ML,, ML, and H_|ML,. The 
Overall ligand protonation constants are 1.0 x iat 


ancien) On x 107°. The formation, cons tantsior’ the 


complexes are 1.0 x Oa Pe Oe Hite AREY Ue ‘os and 


ee Ox ron respectively. Cu was. constant) at 1. 0ux 410” 


M. (See Table 3.) 
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indicate possible. 


A least squares polynomial fit of Ch vs Ch data will 


not improve the situation. First, if the function of Cy 

vs Cr, is linear and has some random error, the polynomial 
will try to fit even those errors, and second, if the 
function of Cu vs Ch is nonlinear it is not likely to be as 
simple as a polynomial of low degree and a systematic 
distortion of the actual relation will occur. 

Consideration must also be given to problems which may 
arise in titrations where the Cy may become small. This 
will happen if the ligand concentration is low, or if the 
metal concentration is high. Since Cu is calculated from 
the difference in two large numbers, the total acid added 
to the solution and the total base added, it may lose 
Significance in some parts of the titration curves. No 
problem is encountered if only one out of several Ch 
values at different c,'s is small, which is a common 
situation, but difficulty is encountered if the C's 
corresponding to several Ch values are small simultaneously. 
Each Cy has limited significance, and therefore the aC, in 
the (ac, /ac, ) quantity has low precision. In general, small 
values of (dc,,/dc,) and (dc,,/dc,)) will always have low 
precision, but since they also contribute less information 
about complexation in solution they should be avoided. 

In summary, with data of limited accuracy and in 


limited numbers, the best results are obtained if the C's 


are of reasonable magnitude, if (dc,,/dc,) is of reasonable 
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Magnitude and if the change in C.. with Cy is linear. 


H 


Figure 1, mentioned above, is the variation in Ch with 


2! ML. and 


H_jML..- At pH 4.0 very little complexation is Occurring 


Ch for Series D; /ancluding the complexes ML, ML 


and Cu increases continuously with increasing ligand 


eoncentration. Each increase of CL adds protonated ligand 


Species to the solution, and therefore the concentration of 


titratable protons,C,,,increases. Above pH 5.0 complexation 


becomes appreciable. At Ch concentrations below 2.0 x One 


M, changes in the ligand concentration have very little 


effect on Cue This is because the ligand is bound as ML 


and ML. complexes, which have no titratable protons. Near 


pH 6.0 this situation is severe, and the problems discussed 


above may become significant; Cu is nearly zero, and 


changes only slightly with C This problem is avoided if 


L° 
the titrations are performed at an excess of ligand, in 


this case about 4.0 x 107° M. 


The curves above pH 7.0 have two regions. Below a Ch 


rer e210 3 fon M, the changes in Cr with Ch are governed by 
the complex H M15: Increasing Ch frome. 0n x ror 4 to 
2.0 x 107° will increase the concentration of H_,ML,, and 


since this is a hydroxy species, it effectively contributes 
a negative number of titratable protons, and Cy becomes more 


negative as the concentration of H_|ML, inereases Barbeyond 


a Cr renigy WALN OMD< 107° 


metal is completely coordinated. Addition ofiiagand 


no further H_,ML, may form since the 


contributes protonated HL species to the solution, and Cy 
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increases. At pH 10.0 the ligand is completely deprotonated 


and Ch no longer increases with CL beyond 2.0 x ike M. 


Figure 2 shows the C_ vs C,, plots for a system similar 


H 


to’ Series C, composed of HL, H.L, HL,HML, HML., and ML... 


3 2 2 Z 


The details of the simulation parameters for these titra- 


Erons are Shown in Table 7. At pH 4.0 only slight 


complexation is qccurring ; . increasing Ch adds protonated 
ligand to the solution, and Cy rises continously. At: pH 7.0 
however, complexation is occurring and the change in C, with 


Ch is more complicated. At Cc, below 1.0 x 10> M, Cu rises 


with increasing Cy because each addition of CL increases 


the concentration of the complex HML. No other complex 


may form because there is a large excess of metal in those 


regions (C,, = 1.0 x 10°} When C, becomes larger than 


x. moe? each addition of ligand converts HML species to 


HML., and ML, - The only reason HML existed was because 


there was insufficient ligand to form HML. and ML. - The 


has no protons and therefore Cu falls with 
3 to 2.0 x 107° M of 


species ML. 


increasing CL in the range 1.0 x 10. 


-3 
ligand concentration. Beyond a Ch Oise .Ous ee Mesene 


2 


addition of ligand.in excess of that. needed to coordinate: the 


metal adds protonated ligand species to the solution, and 


Cy increases continuously. At pH 9.5 no HML species will 
-3 
form, and the maximum in CH ata Ch Ota Orc ee has 


disappeared. In this case all of the added ligand is 


immediately converted to ML. in spite loft the vexcess meral- 


Beyond the complete coordination of the metal, the C,, rises 
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Table 7. Simulation Parameters: Variation of Cry with C 
and Cy, for Figures 2 and 3 
Complexes Formation Constants 
HL 1.0°% lor, 
H5L tia Oiaeyl Oe 
HL fh oe ee 
HML LO Ux 107° 
HML., im 3 ne 
ML. nels a uate s 
pH Range: 3.0 = 10.0 
pa eintervalie JO. Ci) data points) 
Initial Volume: 300 ml 
Titrant Concentration: no titrant 
Component Figure 2 Figure 3 
Ge varied Asis Om 
Cc a ase 107° M varied 
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Figure 2. Cy vs Cc. for a system containing a triprotic 
ligand and complexes HML, HML. , and ML... The overall 
LiGanda protonation) Constants are 1.0 x TOs s SE § Mio 101® 


and sO to." The formation constants of the complexes 
are 1.0 x Tyee. des 0a 3% Lone anded Ox noe respectively. 
Cu was constant at 1:0 x 10° °M. 
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due to the addition of protonated ligand to the solution. 

The Cy vs Cu plots for the same system as shown in 
Figure 2 are shown in Figure 3. The total ligand concentra- 
fron is fixed at, 4.0 4x 107°, Sor thar the farerigncunand 
SroewOL the figure is ai:1 Ligand to metal ratio. “Atle 
4.0 the values of Cu decrease continuously as Cur is 
increased. This occurs because each addition of metal con- 
verts some of the ligand, perhaps present as HL GAG HoL, to 
the complex HML. The average number of protons available 
to each ligand therefore decreases and Cis decreases as 
Weide “At Yoo” 7.0 the: curve: hasitwo segions. Below a metal 
ecoucentrabtion: of \2.,0)\x 10°> there is a large excess of 
ligand (CL walsiv4s. 0 ox nO ak) and each addition of metal 
favours the formation of ML,. Ligands which were in the HoL 
and HL form are converted to ML., and again the average 
number of protons available to each ligand decreases. Be- 
yond a metal concentration of 2.0 x Nays each addition of 
metal favours the production of HML, and liberates ligands 


from the ML, complexes, and Cy beginssbo .increase with 


2 
addition of metal... .At pH’ 9.S,. asi meéntionedwabove,- the only 
complex which is possible is ML... At a low metal concentra- 
tion each addition of the metal converts free protonated 
ligand species to ML, and C,, falls. Beyond a C,, or 

3. 


20 x 10) Miall of tehe ligand is “held” as ML. and each 


addition of metal only adds free uncomplexed metal to the 


solution, which has no effect. on Cue 


These plots indicate that there are regions where the 
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the legend of Figure 2 . 


2.0 330 4.0 
CX 10° (M) 


C, was constant at 4.0 x LO 
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Cy vs Cu for the system of complexes defined in 
M. 
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Ch vs Ch or Ch vs Cu functions may have severe curvature, 
and though such areas are of interest to investigate the 
complexes forming in solution, and potentially have a large 
amount of information in them, difficulty is encountered 

in experiments designed to give (aC,,/ AC.) and (dc,,/dc,,) - 


In these cases linear regions are preferred. Cry vs Cc, and 


Ch vs C,, will be linear if the ligand concentration is 


M 


sufficiently high that the metal coordination is completely 


Satisfied at all pH's. In other words, the Cc. to Cur ratio 


should be larger than the L to M ratio within the complexes 
predominating in any pH range. 


In summary, (dc,,/ac may be difficult to evaluate 


x) CDH 
with reliability throughout wide pH ranges. The most 


severe problems are avoided if: 


(1) C.,, is maintained at a reasonable absolute magnitude 


H 
(2) (dc,,/dc,) are of a reasonable absolute magnitude 


(3) vegaons of cumvatine an C, vs Cy, are avoided. 


Bane -integrat 


The values of (dc,,/dc,) usually fall within +n 


C, PH 


where n is the number of protons on the fully protonated 
ligand. There are no regions of extreme curvature as seen 


on: the Cy vs Cy plots. Mathematical treatment of these 


curves is therefore straightforward. Figure 4 shows the 


values of (dc,,/dc, ) and (dc,,/dc,,) for Series D (Table: 3) 
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(dC, /dC,,).(dC,, /AC, } 


40 6.0 8.0 10.0 
PH 
Pugure. 4 .  (dG../dc_) > ee conch (ae ZdGe) 1, Be asia 
H L Cur pH H M Cl. PH 
function of pH for the system of complexes defined in the 
legend of Figure l. Cu = 130 x 107 °M, Cy = 4.0 x 107 °M. 


(See“Tabler 34.) 
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is (dc, /dc,) and (dc,/dc,,) for Series E (Table 4) where 


3 3 


Sor 100 ehdy Cc. ee 1 One, 


L 
The curves shown in Figures 4 and 5 represent the 
derivatives of the p[L] and p[M] functions with the pH 


(eon, 12). 


dp[L] 
dpH 


dC, 


= (a2) 
dc, 


Cyr PH 


A simple series manipulation of this equation leads to a 
relation which allows the calculation of [L] and [M] from 
the experimental values of (dc,,/dc, ) and (dc,,/dc,,) - 


Integrate both sides of Equation 12 from PH, ee pela 


pH pH 
adp[L]\ dpH dc,, dpH (13) 
dpH i dc 

pH) pH LARC pe 

Simplify the left hand quantity. 

PH pH | 

EhS Ne ea ar dpH (14) 
dc 
PH, H L! Cy, PH 


Integrate dp[L] from PH, to pHs to get PIL], at PH, and p[L] 


ai pHi. 
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dpH (15) 
Cyr PH 


p[L] = Se MON 


PH, 


The integral may be evaluated in one of two ways. 
First , as discussed. above, the entire set of (dc, /dc,)'s 
may be fit to a polynomial, and the integral evaluated by 
simple calculus. Second, it may be evaluated numerically 
using Simpson's parabolic rule, or the trapezoid rule. 

In this work the trapezoid rule was used to evaluate the 
integral directly from the tabulated Gatos kucue DH 
function. It was found that this procedure, used at 0.1 


pH increments, gave. as) good accuracy as is in most cases 


justified by the accuracy of the experimental data. 
Formation Constants 


The FICS method provides free ligand and free metal 
concentrations as a function of pH. Formation constants 
can be extracted from a linear least squares fitto the 
mass balance of each component in solution. In sucha 
procedure, the sum of the square of the residuals, U, is 


minimized. 


+ yy (a) °(119[mj]7) se ae 
c= c. - ( cee (16) 
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In all the examples presented in this thesis a Gaussian 
elimination (57) followed by matrix inversion has been 
used to solve the simultaneous equations. The mass balances 
for the ligand and metal were treated individually, or in 
some cases together, while the mass balance for protons 
Was NOt Used to calculate formation constants!" "In" this 
respect the FICS method has considerable advantage over 
Other pH titration methods. Several problems are inherent 
to the use of Cu in the calculation of formation constants, 
and these are avoided with the FICS method. The free 
hydrogen ion concentration is generally difficult to deter- 
mine experimentally, involving either an estimation of the 
hydrogen ion activity coefficient, },+ , or an electrode 
calibration in terms of (Hal. At low and high pH's [H’] 
and [OH ] constitute a large portion /of Cyr and if there is 
uncertainty in their values this will produce a large 
uncertainty in the concentration of bound protons. In 
comparison however, the FICS method does not require [H™] 
at any stage. The electrodes may be standardized with 
NBS reference solutions, and an estimate of ess is only 
required if the mixed activity-concentration constants are 
converted to concentration constants. 

For simplicity in this work, the notation pe will be 
used if the constants are evaluated from the ligand mass 
balance, > M if they are evaluated from the metal mass 


balance, and >> ML if both mass balances are used simul- 


taneously. 
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Some estimate of the accuracy of the formation constants 
obtained by the FICS method is necessary. The simplest, 
though not the most reliable quantity, is the goodness of 
the fit of the model to the known total analytical 
concentrations of ligands and metals. The experimental 
values of [L] and [M] and the formation constants are 
inserted into Equation 16 and the standard deviation of the 
tates tyrone ee where n is the number of experimental 
points. The relative standard deviation of the fit, in 
percent, has been shown for several examples in this work. 

The standard deviation calculation is limited in that 
it does not give any indication of the accuracy of the [L] 
and [M] data, rather it only estimates the fit of the 


hypothesized model to C_ and Cy: A better measure of the 


L 
reliability of the model and the formation constants can 
be had by simulating titration curves and by comparing 
them to the original data. This was done as a final check 
with the metal-glutathione complexation study reported in 
Chapter VI. However, it is impractical to use routinely 
while attempting to find a model for complexation because 
the simulation procedures are slow and, in addition, may 
run into difficulties specific to each model proposed to 
describe the complexation. 

The suitability of a model describing the composition 
of complex species existing in solution may be checked in 


several ways. It is expected that a good hypothesis will 


give a low standard deviation of the fit and that the 
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constants derived from one mass balance will agree with 
those derived from another. The constants which are known, 
for example ligand protonation constants, must be derived 
correctly from the data. The existence of some species 
may be questioned if small or negative formation constants 
are produced, or if the agreement of values based on 
different mass balances is poor. 

The next section will deal with the effects of various 
types of errors, and the problem of assigning a model to 


the system. 


Se. The Effect -oL Error on the FICS Method, 


This discussion will be limited to error introduced 
through the experimental procedure. The systems of 
simulated titrations discussed to this point have been 
designed to represent feasible experiments, but without the 
error associated with the actual experiment. In this 
section on errors, several systematic and random errors 
inherent to a real experiment will be considered. The 
error in [L] and [M] will then include uncertainty arising 
from both the FICS data handling and the imposed experi- 


mental error. 
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Badiaition Error 


Three separate types of error are introduced by 
dilution. First, during the series of titrations necessary 
al 
to evaluate Som onl AORyDH the values of Cy are not constant 
Betevary according tothe quantity of titrant added to, that 


PH. This is easily compensated for by using C. exactly as 


L 
it is found at that pH, based upon volume of titrant added 
eee that point... Each, pH.will have a unique: set of c.'s 
through which the derivative is evaluated. Series A and B 
in Table 1 were simulated including significant dilution, 
and the values for [L] shown in Table 6 agree very well 
with the expected values. 

The second and third dilution problems are related to 
dilution of components other than Che and cannot be 
corrected in the data handling process. During a set of 
titrations designed to give the derivative ake, Sora ce pH’ 
the Cur Will vary! sthrough, dilution. Assuming that, all the 
titrations, have equal dilukion to a given pH, the derivative 
will be evaluated ata Cu different from Cu at PH: 

The. third typecot ernortis) implied above... All ithe 
titrations do not have equal dilution to a given pH, and 


consequently, at each pH, Cur will be different for the 


ait perent titrations... [he wterm (dC,,/ AC, explicitly 
Cy Oe 
M 


requires that the change in Cu with CL be measured at 


constant pH and Cur 


These problems contribute significant error to the 
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free ion concentrations. Table g shows [L] and [M] for 
Series C titrations with and without dilution. With a total 
difution of 7.201 ml titrant into, 300 ml, the free ligand 
concentration at pH 10.0 had about 10% error and the free 
metal concentration about 25%. Without dilution the error 
tee tt] and) [M] ‘was only 0:5%.5 Table) 9 summarizes 
formation constants derived from these two sets of free ion 
concentration data. 

The second and third types of dilution error, discussed 
above, cannot be corrected once the data has been taken. 
On the other hand, it is possible to experimentally ensure 
that Cu remains constant throughout’ a set of “Citrations, 
Cither by eliminating dilution by use of coulometric 
titration, or by simultaneous addition of titrant and a 
concentrated metal cota om Computer automation of the 
latter is a practical alternative, and was used for the 
studies of complexation reported in this thesis. Details 
of the manner in which this was carried out experimentally 


are reported in Chapter IV. 


Pecors in Cy, and Cu 

The experimental values of Che Cur Cy and pH and so on 
may have systematic and random errors associated with them. 
Studies with simulated titration data show that systematic 
errors give rise to the most serious problems with the 


FICS method. The nature of the calculation allows small 
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Table 9. A Summary of the Formation Constants 


Calculated by FICS from Series C and Cc! 


Formation Constants 


Complexes Known No scbateries abt 

HL 1.00 t? . PeOgOs sR Lote an Vio eae nae 
HL Pete. 1.00620! 10” 1,073%x 107° 
HL 1.0 x 1029 9.9996 x 10? 1.215 x 1019 
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errors to propagate and give rise to substantial error in 
the final [L] and [M]. The error increases, in most cases, 
from PH, to the final pH at a rate dependent on the size of 
the error on a point to point basis. The magnitude of the 
error seldom increases, then decreases at some further DH. 
The FICS calculations appear to have a tendency to 
accumulate error even though the error appears to be random 
in nature. 

Systematic error in quantities such as Ch and Cu may 
arise in two ways. First, there is an error associated 
with the concentration of the stock solution from which the 
samples are prepared and, second, an error associated with 
the Se Oe of the samples. The latter is usually 
considered to be a random error, but in the FICS method it 
will have the same effect as a systematic error. This is 
because the derivatives (dc,,/dc, ) are repeatedly calculated 
through the pH range using for example five titrations of 


varying C If some error in (dc,,/dc, ) arises because of 


L° 
an error in the sample preparation, that same error will 
appear in (dc,,/dc, ) at every pH where the calculations are 
G@arrica out. TO minimize this source of exrxrOn, as, many. 
titrations as possible must be performed, the samples must 
be prepared with utmost care, and measures must be taken to 
ensure as much randomization of errors in the sample 
preparation as possible. For all the experimental parts of 


this thesis, the samples were prepared in the OLIGSE S712, 


4, 1 and 5, where 1 would be the lowest concentration of Cr 
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and 5 the highest. In this way, the titrations performed 
first were in the middle of the Ch range and errors in 


sample preparation had somewhat less effect on (dc,,/dc,). 
Error in pH Measurement 


Unlike other pH titration methods, the FICS procedure 
does not require that the hydrogen ion concentration or the 
hydrogen ion activity be known to get accurate free ion 
concentrations for ligands and metals. The pH is merely a 
reference scale on which to first evaluate ton eoSn ce, BH 
and then to integrate this quantity from PH, to the 
various pH values. Its numerical value is not required in 
any stage of the Seni ion of [L] and [M]. However, the 
PH does enter into the calculation of formation constants. 

If the electrodes are standardized with NBS activity 
references, then the formation constants derived are mixed 
activity-concentration constants, and if the electrodes are 
calibrated for concentration, the constants are concentration 
constants. 

The FICS method requires only that the pH electrode 
standardization be internally consistent, and serious 
systematic error is avoided. Random errors in individual 
measurements will exist however. These will give rise to 
exactly the same type of error as random errors in Cyr and 


will therefore be discussed in the following section. 
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Errors in Cy 


The effects of errors in Cy are quite unlike those of 


exrors “in Ch or Cure 


data point, being the combination of an error in the 


Eryvors "in Cy Vary Lrom data porne, co 


quantity of titrant added to that point, and an error in 

the pH measurement. Unlike the results of an error in Che 
where every pH is affected more or less equally, the errors 
in Cy will be different at every pH value in every titration. 


Therefore, if errors in C, give rise to a (dc,,/dac, ) too 


L 
Harge by 0.05 at pH 5.0, the same error will appear at pH 
G20, and so on. The. errors in (dc,/dc, ) caused by 
Gpcertainties an Cu will be different at every pH. 

In order to evaluate the effects of error on the FICS 
procedure, random error was imposed on simulated titration 
data to give some qualitative idea of the order of 
magnitude of errors the method may endure. However, since 
real error may not behave exactly as assumed here, 
conclusions must be drawn with some care. 

Two different approaches are considered. The first is 
random error based on the magnitude of Cue Chateis,* the 
error is larger at larger C,.- THIS tS SilliadaefeO,d Constant 
sized error in pH. The second approach is a random error 
of constant magnitude. This sets the effective sensitivity 
of the Cy eee rena As Cy becomes small it loses 


significance, being the difference between two large 


numbers, total acid and total base concentrations. 
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Table 10 summarizes results of the FICS method on Series 
D data, including the free ligand and free metal concentra- 
tions, the constants calculated using the free ligand and 
free metal concentrations, and the standard deviation of 
Poe -Lit,-invpercent. Random errors of 5%, 2%/ 12% and 0252 
Of Ch have been imposed upon the Cu data which was used to 
give the results in Table 10, and the effect of these errors 
is shown in Table 11. The free ligand and free metal 
concentrations are accurate only if the error was less than 
0.5%. The constants calculated using the [L] and [M] values 
from Table 11 are summarized in Table 12. It appears that 
errors, OL’ 0:53" in Cy will give very unsatisfactory results. 
These are random errors, however, and a systematic error of 
Dy. 5% in Cry will give much worse results. 


Table 13 presents results of imposition of absolute 


errors of various magnitudes on C,. data from Series D. The 


H 
SBrror is meported in parts per thousand of Cy where Cr 
fe) O 
: -2 
is the value of Cy at PH. - ae Cry 2 Sa erecr LO M, then 
fe) 


the random error is 1.0 x apiee or less at the 1 ppt level. 


This means that a Cy Vaineroxr )1..0lix Toe. has only one 
significant figure of accuracy, and: aL. Og ‘Mem has virtually 
none. The results indicate that even up to 4 ppt error, 
significant problems do not arise in free ligand and free 
metal concentrations. Table 14 indicates that the errors 

in formation constants are moderate at the 4 ppt level, and 


acceptable at 2 ppt. 


In summary then, the two types of error imposed on Cy 
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Table 10. The Results of the FICS Calculations on the 


Simulated Titration Data, Series D 


Free. Ligand 


pH Known 


a0 = «1.22045 
Bee L250) 5x 
0s 5.7920 
15 SEAN EES Opa Ni 
Beor «1.422% 
Be. , 1. S44ix 
oe 1. 980"'x 


Complexes 
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HL 
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2, 


Ve 20 2e.x 
tet 46 Xx 
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~ mw mM MK 


Known 


120 
10 
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x 102° 


sc. 0 


x mw KM MK 
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Calculated 


Or 
AOD 
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3.414 x 10 
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ALSi 4, x2 108 
Tadd: (x-TOe 
9.959 x 10 
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Free Metal 


Formation Constants 


Sum over C 


9.983 &x 
i Toate bes ay 
yep oak 
SPARES) 
ge MBS 
1.004 


x Ke Mm 


L 


10° 


Calculated 

poms 419 sone 
20746 menton 

Paha 389 ee Ton 
10 7.262 x 10719 
rae b Sane bis bee x de Sle 
T1469 .474 x07) 
12 9.510 x 10722 


Sum) over sC 


1.004 
9.948 
1 OS 
OZ 


M 


Sy) 


si ee lh Si 


p ad ee 
* ie y : yl “ee 
" ryere)* 
‘ i ‘ , ala oat eat ins 
Ley) co teneNe “’ ; a, i : 
pes! ey ftir a wame 
ea reaveeare weed ae " heb 
mi te 
a : bo | 
Un m me ae (1 ‘ i ah 
wt a i # uM g | a a . 
Bede ula ry eee 
Da fe ae Pa we li, 
| if tH hiths a, 
weg wi | 
v.16 i, 


i 

iyi 

: : 

7 | U 
eb} i 


X ne ' i } : de’ 
¢ 4 jaf j ‘ 7) : ee ae i 
» ‘ oer : f pareeariges Wee nara 
r . } . 
ee 
Pi Le 
oe : 
“ie ie " i 


ft Lae at | 


. rh ie ci 


ie aed 


* 7; 
ie i) if 
ss } oe Rey 
mm j 
Frei 
y re ae | 
% 


. oa 
¢ é Ptr " 


if mi secs aot aera 


we) i te 


lla 


pa win ia iN 


& igh mi 


er 


Wig ha. i” ie : 


a 
ak 


mn 
AN 


Pik 


Pry yee mo) 
jib ee a a, 
arate | Mistted 
her hae 


ray 


pas Py 


“ure 


ne ie ie 


dhs tigi we A} 


my) 


ayy 


ok" 
(he 


Pat 
"et, 


- e % i . 
aos ee 
ie Re ae Poe 
mbit lt Lay n ; 
=9@ 7 7 a F ig : 
a A ss 4 x ey Pe ix) 1 


Table ll. The Effects of Random Errors in C. on the 
Free Ligand and Free Metal Concentrations 


Calculated by FICS® 


Free Ligand 


pH Known +5b +22 +13 +0.53 
We) 1204 )% 107° Meas ea tig 2) eed 
pmo) 18250) x 101” 1 G2 1.38 ge 7 
pmo) 58792) 107" 6.80 6.4 6.00 BEG 
mmoe 1 f055rx 107 1.24 peg 1.09 1207 
emoe 12422 ¥ 107° 1.68 el 40 1.45 
omon eda x doo sae ec 1.92 1.88 
homo) 1.980% 100° baer po 2.06 2208 
Free Metal 
4.0 weetate 10)” 8.08 4.82 3.40 3.23 
emu 22722 x 1077 6.80 B85 94 252 
pmo @ 40314 x 10.- 8.34 BuGT 4.99 4.17 
pote Test 10729" 12.20 9.50 eee 6.95 
pmoee 21959 x 10 7° 3,82 Py IS pA 1.88 
Pee h4c3 x 10s- id, 91 3225 AEG rome, 
10.0 DRAGS ee OR A 5.40 Bi 2.70 2.47 


“The random error in Ce. is proportional to the magnitude of Cy. 


Dohe size of the random error, in percent of Cue 


“same order of Magnitude as known ligand concentration. 
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Table 13. The Effect of Random Error in Ch on the 


Free Ligand and Free Metal Concentrations 


Calculated by FICS 
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indicate that the FICS method requires accuracy of Ch of the 
order of 2 ppt or less for most reliable results. It must 
also be kept in mind that the ligand to metal ratio, and 

so On, used in the Series D data simulation were chosen 
with all those factors in Part B of this chapter given 
careful consideration. These conclusions cannot be applied 
to situations where C,, is unnecessarily small or where 


H 


there is severe curvature in the function of Cy Vs Ch 


The Effects of Errors in the Initial Free Ligand and Free 


Metal Concentrations 


Since the FICS procedure requires free ligand and free 
metal at some starting point, this imposes restrictions 
which contradict many of the experimental requirements 
discussed above. Free ligand can be calculated from the 
total ligand in the presence of metal only at low pH, where 
protons compete successfully for ligand. However, a large 
excess of protons means that Cry will be large, and Cy will 
therefore lose significance at 5) PSOE ree oe higher values, 
bringing about the problems discussed in an earlier section. 

Table 15 illustrates the effects of error in PIL], and 
PIM]. on the free ligand and free metal values derived from 
Series D data. The exact [L] and [M] may be found in 
Table 10. 


The far left column in Table 15 is the set of [L] and 


[M] calculated using PIL], of 6.439, PIM], of 3.000, based 
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Table 15. The Effect of Systematic Error in pX, on 
Free Ligand and Free Metal Concentrations 


Calculated by FICS 
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63 
on the assumption that no complexation of metal has 
occurred, while in fact about 25% of the metal ahd 
coordinated to ligands. There is consequently Significant 
error in’ [lL] and [M]. This problem may be*corrected* with 
the following algorithm: 

(1) assume no complexation, calculate PIL], and PIM] o, 

(2) calculate [L] and [M] using PIL], and PIM] | from 
C44) (for example Column)1)- Table 15)? 

(3) calculate formation constants based upon these [L] 
and [M] values (Column 1, Table 16); the entire 
set of [L] and [M] may not necessarily be used, 

(4) calculate new values for PIL], and PIM], based on 
the estimates of constants derived: best estimates 
for protonation constants or any known constants 
should be used, 

(5) determine a new set of [L] and [M] based upon the 
new PIL]. and PIM]. (Column 2, Table 16), and 

(6) determine a new set of constants, and a new 
estimate of PILI, and PIM], (Column 2, Table 16). 

As the results in Tables 15 .and 16 ‘show, the [i] and 

[Minitaprove; until they anemone ox lessvexacts As» the 
PIL] , and PIM] | are improved, the constants are more 
reliable, and change less with each iteration. 

This method will fail if too much metal or ligand 


are complexed at PH. : If the types of complexes at low 


pH are not known, the procedure Wall LA eSOUwaad . 


Values of PIL], and PIM], may also be found by a trial 
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64 
and error approach in which the standard deviation of the 
fit is minimized. Each p(L], and PIM], is used to get a set 
of formation constants and these constants, with the [L] 
and [M] data are fit to the total ligand and metal mass 
balances. The residual defined by Equation 16 will reach 
a minimum when PIL], and PIM] , provide a set of [L] and [M] 
most closely matching that of the model of complexation. 
This procedure will not be meaningful if the model for 
complexation is incorrect. This method has the advantage 
that it choses the best PIL], and PIM] | in terms of the 
entire range of pH, but the disadvantage that a model of 
complexation over a wide range of pH must be known. Both 
of these approaches have been used with experimental data 


and give good results. 


Titration methods do not in general provide information 
from which the model, or the types of complexes in a system, 
can be derived. In this respect, FICS has an advantage over 
other methods, in that the FICS results can be used to 
provide some guidance in the selection of a model. Using 
(hl, (Ul -andeim; thetconcentrations of proton, , ligand and 
metal whose chemical situation has not been determined 
(AC, Ac, and ACcy) can be calculated. 

For example, given the total ligand concentration Che 
and the free ligand concentration [L], AC, equals (CL = ELS 


If the pH and the protonation constants are known, 
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65 
Ac, = ¢, - (inl + Aytwlin) + 8,187 (1) (17) 
L L us 2 


If the formation constant for HML is known: 


ac, = ¢, - (inl + Aymin + Boon7tn + Boom iim) 
(18) 


Figures 6 and 7 are combined plots of ACs AC, and 
ACy for the Series E data (see Table 4). These quantities 


were calculated from Cur Cc and the pH, [L] and [M] 


i Sy Cy 
from Table 17. Three distinct areas can be seen: 

(1) below pH 3.5 there is a region where AC,, is about 
twice as large as AC, or ACy- The plot indicates 
that this complex reaches a concentration of about 
Ue Ogos 10-7 M, and then is replaced by another 
appearing in larger quantities. 

(2) the predominant complexes near pH 6.0 are composed 
of more metal than protons or ligands. The ratio 
Of Metad..to protons or ligands is about two’ to one. 

(3) above pH 6.0 the complexes which predominate 
contain only metal and ligands, in equal proportions, 
perhaps ML or MoL.. 

Figure 7 is an expanded part of Figure 6 and shows that 

the first complex to form at low pH must be H,ML. Using 
the, (i) and (M]" data. from pHs fo05 the formation constant 
18 


for this complex was found to be 4.0 x 10. 


The ACur AC, and ACy values were calculated using this 
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Figure 6. ACuy oO; AC, 4; and AC 

complexes defined in Figure 5. 
He cr 
and [M] values in Table l7. 


Calculated directly from C 


for the system of 


Meo: 
These quantities were 


and Cur and the: pH pa3| Gi] 
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69 
estimated constant. The AC, plots are shown in Figure 8, 


for several guesses of that constant. The first guess, 

PaO x MO. was too low, since AC, remained large. A guess 
oe. OX Tere proved tobe” too~large; giving negative AC, 's. 
ane guess of 1.5 x sre was also too large. The value 

mee xX ton was good, giving a AC, of nearly zero from pH 
Boo, to 3.5: 

Figure 9 is a combined plot of AC. AC, and AC, using 
pene initial guess of 1.2 x FON" for HoML. The complex near 
PH 5.0 is almost exclusively of a HM,L type stoichiometry. 
The free ligand and free metal data up to pH 6.0 was used to 
get estimates of formation constants for the two complexes, 
H.ML and HMoL. These were found to be about 1.02 x io}? 
and 9.98 x rhage respectively. 

Once again a combined plot of AC. AC. and AC, was 
prepared using the estimated constants for H,ML and HM5L. 
The curves are shown in Figure 10. The estimates were 
very good, giving small ACh AC, and ACy Values. up to pH 6.0. 
The final complexes clearly contained no protons, and were 
composed of equal amounts of metal and ligand. The free 


ligand and free metal data for the entire pH range was 


used to get formation constants based upon the model H ML, 


eye. 
HML and MoL, (see Table ) 


Refinement of the Model 


Given that some model of the system is proposed, the 
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final procedure must determine the minor constituents in 
solution and distinguish among species having the same 
stoichiometry. , Series E, discussed above, will be used as 
an example. 

Table 18 summarizes formation constants derived from 
the same [L] and [M] set, based upon the series Eb.) Column 
I is a model based on HoML, HM.L and MoL., and Columns II 
through VI are similar except one or more additional species 
have been included in attempts to determine which, if any, 
minor constituents exist in solution. 

Column II includes ML, ae the standard deviation of 
the fit somewhat improved. However, the formation constant 
for ML is small, with the polanueiesy MoL, present in large 
excess. In addition, the formation constants for ML based 
on the mass balance over CL () L) and over Cu () M) are not 
in as good agreement as for the other species. These 
considerations suggest that if ML exists in solution, it is 
in small quantities and the constant has low reliability. 

Column III is for a model which includes homologs of 
HoML, such as HML and H2ML. The constant for HML is 
negative in the ™ evaluation. The constants for HML are 
not very consistent, based upon the)L and >M. in, addition; 
standard deviations are not improved. In Column IV the 
Species H,ML was rejected, and HML remains. The constants 
disagree, and no improvement in the standard deviation has 
resulted. 


Column V includes the complex MoL, and Column VI 


re 


age es ae ee a 


ona ‘Bea esi Lis Be oe 


Do 

a ‘ eet _ 

rick r r 
Se Ten 


peer Bh aw 


7 


ek Sec a 
ee nd 


ve ; “" vis niet \ at ty 6a areca: 4 as a mite ‘ Jeon 
Lip gah dua ‘eae rel ‘heue esa 
eh aramid Iles) 's aka me ane yt ie 


‘ 2a. 4 Legs oe bate era a a) NST rn “eal ud ie i) “ 


ey eS, areas PER Is | fide PAA aes Be eit! mers pee inh 
’ pon cn 


Choo) oe ae if a, Sa hua ne 


er ue oY aed Rhee dale | lea a ra tng A | wilt <8 ce a i fh hots 
Mey eis) Qa tat Sa ae 1 <0) Wieeren hy “a »4 her ae " 


= Ne i) "by; 4 ‘ a 59g I ry, 7 & 3 Lat va re ‘ Hi) “ea tantly tin 


me 


baleiest tod? ain ae Rae Hee sai by Wit | ane ei ic 
te 2 ie Me M 4} pe? J "oe SAA ty i ae (r) Oe hics ite ha ty bissahia 


1? Vaprap ine nels : ni sah bana » ong fuss 


\ 


it Deeb aes kts mina snag wg aa sag 
: bie: “hei bd aOR an pedi f Ceylevoe dh jp at ve ote 1 
‘e i, igh 10% | “yous ar: i; xg ay ry pe ‘atta, +. 
salt me af Sault ibs ume — ; ¢ GM tuk hen a gota) hi : 
pnace ness ol Ro Bu Ge tga ‘tapes enantio 
eh We oped toa, oe. igawaga Son: oe ecirar 
a aryrnbe At fee bee, peat bir panels any 


sjpie! wee ao i patna wie ‘ad i, snstheteieont cian 


Sam ene as oe et a 
: ae l i Mi “in vs cf f 


74 


S0€0°0 
gy Ol xX 8S8°6 
gol x 000°T 
Byte O00} 
x e 
oA SLE 
Wd 


s9T0°O 


gp Ot Xx GL6°6 


Se @ 
g Ol OcO'T 


erut * 196°6 


Il 


$€60°0 $810°0 Wa Jo ‘a's 
Cr Cig 
OT X SOO'T OI X ZT0°T Coy 
LT LT 
: Z 
xX e 
g OT ¥ ¥00°T ot * TE0 4 7° WH 
oH 
my 
: : Z 
gt0T X £86°6 g TOT X 826°6 IW H 
WH 
TW 
501 ¥ 7L6°6 TH 
e ‘A 
rid: Saneore oH 
rz0T ¥ 00°T qe 
WX rT SaxeTduioa 


STSpoW ZUSASTITGd AOF ejzeq 


TeIOW SS7T pue puebT] eerq worZ sque jasuos uoTzJeWIOT “gt eTqey 


ys. 


eat 


as 


6T 


tL 


0S0°0 


Jt S00 7 


OT X 7S6°6 


GE = LLO-T 


OT X P7S*6 


WX 


sLT0°O 


iyOl x tto tf 


xX 6 
gt Ot GcO Tt 


gi Ol x GL6°6 


OT X P69°E 


OT X 796°6 
GOT * ¢LG 6 


O12 4 00U.T 


of 


xX ® 
izot SO UmeL 


izot xX 776°6 


0c 


xX e 
60! OcO'T 


xX @ 
zzot S054 


W 


OT *.00 Tt. o= 


SCeOnO 


OT ASE UO a: 


HE x=6605-. 


OT X LLE*T 
OT X S80°T 


UL A=OTE 6 


"IWH 


SexeTduoa 


(penuTzUOD) BT eTqeL 


76 


$Z700°0 
xX e 
tol SO0O°T 
pole cO0-T 
xX @ 
6 Ot SO0O°T 
Wx 


a ES 


$LT0°0 
Tz0T X Ovs*T 
x< e 
poe xeeLOt 
gy7Ol x 00" T 
x c 
TOL XLT Bt 
g0T X L96°6 
ctOT ¥ 7L6°6 
xX e 
tz0T * 00°T 
ae 


%S70°0 
pot ¥ S00°T 
gt OT X 200°T 

gOT X Z9T'6 
xX @ 
peor * 900° 
We 


$8T0°O 


URE Ses eye] BO ey § 


gr OE * O£0°T 


OT * L6S°F 


gol x 61T6°6 


tH 


SaexetTduos 


(penut{UuOd) g7T eTWQeL 


V7 
includes HM.L.- Both of these complexes are rejected. The 


best model in terms of standard deviation is that of H.ML, 


ac 


HM.L, ML and M.L As discussed above, the existence of ML 


2 Dees as 
cannot be conclusively proven. 

The original Series E data did not contain the species 
ML. This illustrates the possible consequences of assuming 
that no error exists in the data. The simulated titrations 
did not themselves contain errors, however, slight curvature 
mir the function. of Ch Vs Ch and so on have given rise to a 
small, but, finite error in the [L] and: [M]. values. The 
only defence against this problem is good understanding of 


the FICS method and of the pitfalls inherent to its 


application to real data. 


E. Summary. 


The FICS method has several important advantages over 
most pH titration methods for the study of complex equili- 
bria. The method provides free ion concentrations directly 
from the experiment. Moreover, they are derived for each 
component individually, that is, independent of the other 
components,and independent of the numerical value of the 
hydrogen ion concentration. In contrast to most other 
methods, FICS derives free ligand and free metal concen- 
trations completely independent of any preconceived model 
of interactions which may be occurring in solution. 


Estimates of [L] and [M] have been determined 
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even if a model for the system is never found. 

The FICS procedure was tested extensively on simulated 
titration data, and found to be applicable to systems 
containing a variety of species. Examples. include protonated 
complexes (Table 9 ), hydroxy complexes (Table 10), poly- 
nuclear complexes (Table 17) and mixed ligand complexes 
(Table 19). At this time no systems have been found for 
which the FICS method cannot be used in exactly the manner 
Gescribed in this work: No modifications ofsany kind, 
based on types or number of complexes or magnitudes of 
formation constants, were needed. 

No attempt has been made to show that FICS will give 
completely exact [L] and [M]; rather, attention has been 
focused on each step of the procedure to find the most 
practical and reliable experimental method to use with FICS. 
In this respect it was found that situations where Cy is 
small and where the curvature of C,, vs C, is: lange must be 
avoided. Systems containing an excess of ligand appear to 
be most favourable. Under these conditions a series of 
relatively simple manipulations of the experimental data 
will give accurate free ion concentrations. 

The availability of free ligand and free metal 
concentrations presents excellent possibilities for the 
logical deduction of the complexes existing in solution. 
For this purpose, the quantities AC» AC, and AC, were 
introduced and applied to the deduction of a model from 


simulated titration data. Considerable information 
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these simulations are shown in Table 5. 
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concerning the model is also potentially available from 
plots such as shown in Figures 1-5. 

The effect of systematic and random errors of various 
kinds has been considered. A study of the effect of 
dilution revealed that significant error can result due to 
the systematic violation of theoretical requirements of the 
FICS method. The effect of error on PIL], and PIM], was 
considered, and a method suggested which will make 
calculation of their values possible even when a small and 
unknown amount of complexation has occurred. Random errors 
of various magnitudes were applied to titration data, and 
FICS found to be reliable with up to about 2 ppt random 
error in Cie 

The major advantage of the FICS method over other pH 
titration methods is that it provides an experimental value 
for {%), and. [M]\«which, <in,.turn,,aids in the ,selection ‘of 
a model and considerably reduces the possibility of having 
to choose among equivalent models. The major disadvantages 
Of the method are the stringent requirements placed on the 
data aquisition, including its accuracy and problems 
involving dilution. The FICS method may not provide the 
best possible values for formation constants, but with 
good data the model of complexation can be deduced with 
some confidence. 

Chapters III and IV will describe the methods and 
apparatus designed specifically for the collection of high 


Precision,data. suitable for.use with the FICS data handling 


methods. 
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CHAPTER, Iii 


EXPERIMENTAL 


A. Materials. 


All solutions were prepared with oxygen and carbon 
dioxide free distilled water. Removal of the gases was 
accomplished either by vigorous boiling, or, if the solution 
was acidic, by bubbling with CO. and 0. Pree vnitrogen’. 
Boiled solutions were cooled in the CO. and 0. free atomos- 
phere of a glove box. The ionic strength of the solutions 
was adjusted to 0.3 moles/kg of solution with sodium 
perchlorate (G. Fredrick Smith Chemical Company, Columbus, 
Ohio). 

Carbonate-free sodium hydroxide was prepared by 
dilution of a saturated NaOH solution (AnalaR analytical 
reagent NaOH, BDH Chemicals, Toronto, Ont.). The saturated 
NaOH was stored a minimum of one month to allow the fine 
particles of Na.CO, to settles; All alkaline solutions 
were stored in, and dispensed from, polyethylene containers 
ina co. and 05 free atmosphere of a glove box. 

Potassium hydrogen phthalate (Primary Standard, J.T. 
Baker Chemical Co.) was dried one hour at 120°C and cooled 
before sample preparation. Individual samples for titration 


were weighed out by difference on an analytical balance 


(Gram-atic, Mettler Co., Zurich, Switzerland). Sufficient 
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quantities were used to ensure a precision of 2 parts per 
ten thousand. 

Nitric acid solutions were prepared by dilution of 
the concentrated acid (J.T. Baker Chemical Co.). Indivi- 
dual samples for titration were weighed out by difference 
using a syringe and an analytical balance. 

Aspartic acid was recrystallized from cold ethanolic 
solution, washed thoroughly with ethanol and dried in air. 
A stock solution was prepared by dissolving the crystals 
and was standardized by titration with NaOH. 

Glutathione was recrystallized from a cold 50% ethanol 
solution, washed thoroughly with absolute ethanol and dried 
in air. The crystals were refrigerated until use. A 
short time before the experiments, a stock solution was 
prepared by dissolution of the crystals. This stock 


solution was stored under argon and in ice at all times. 


The recrystallization and these other precautions were taken 
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in order to minimize the concentration of oxidized glutathione 


in the sample solutions. The glutathione concentration was 
determined by titration with NaOH, Tables 20, 2) sandal? 
Stock solutions of zinc, cadmium, and lead were 
prepared by dissolution of the metals in HNO , with 
subsequent neutralization of the excess acid by carbonate 
free NaOH. These steps were carried out with small volumes 
in order that the final contribution of nitrate and excess 


acid to the sample solutions was negligible (approximately 


6 


0.005 moles/kg of nitrate and 1.0x10 ~ moles/kg of excess 
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Table 20. Standardization of Sodium Hydroxide used in all 


the Metal-Glutathione Complexation Experiments. 


Consumption Calculated 
titration Wawerght KHE, of NaOH, NaOH Concentration 
gm gm moles/kg 
1 On3679 Oot, O53 .L49 
2 023806 SSG) Ooi g 
3 0.4377 Samy ea § OV3S175 
4 0.3184 4.911 One Aa 
5 0.2742 4.230 0.3174 
Average 0.3176 


Relative Standard Deviation 


0.08% 
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Table 21. Standardization of Glutathione for the Zinc- 


Glutathione Experiments. 


Consumed 
Weight of 4 Weight gf Weight of Concentration 
Sample Glutathione, Solvent, WaOH, of Glutathione 
gm gm gm moles/kg 
1 202002 190.002 Eyes WN) 0.09640 
Z 8.000 192.003 ZEO0 0.09647 
S 4.000 195.998 1.284 0.09623 
4 15) 00). 189.002 3.205 0.09631 
5) 6.003 194.004 Tao o2 0.09640 


Average 0.09636 


Relative Standard Deviation 0.103 


Astock solution 


Dcolvent contained 1.159 x ioe moles/kg of HNO 3. 
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maple... otandard: zation of Glutathione forthe Cadmiun— 
Glutathione and the Lead-Glutathione Experiments. 
Consumed 
Weight of Weight of Weight of Concentration 
Sample Glutathione, Solvent, NaOH, of Glutathione 
Leoes gm gm gm moles/kg 
cL 8.440 Lo 7 62 2.426 0.08678 
2 AO. 183.485 2ELAG On08i67 1 
3 7.646 LTT 2956 2.204 0.08681 
4 8.053 199 153 2320 0.08676 
5 8.345 180.094 2.4393 0.08670 
Average 0.08675 


Relative Standard Deviation 0.05% 


“me Solvent contained 1.159 x vor: moles/kg of HNO 
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acid ). Sodium perchlorate was added to bring the ionic 
strength to 0.3 moles/kg. A solution density of 1021.8 
gm/litre (58,59) was used to convert metal concentrations 
from moles/litre to moles/kg of solution. 

A solution which served as a dae or, solvent...1n 
the preparation of samples was composed of 0.3 moles/kg 


NaClO, and a small excess of nitric acid. This excess was 


4 
determined by titration with NaOH. The solvent contained a 
sight. quantity of nitric acid in order that it’ might serve 
as a pH standard at the beginning of each titration. This 
will be considered in more detail in the section on the 
experimental method. 

The solutions used to compensate for dilution by 
titrant during the course of each titration were aliquots 
of the stock solutions used to prepare the samples. The 
samples, in general, contained about 10 gm of the stock 
solution of «each component in a total of 200 gm. The 
quantity needed to compensate for each addition of titrant 
was therefore about 1/20 of the weight of added titrant. 

Ati tration consuming 5 -om of titrant therefore required van 
addition of 0.25 gm of the stock solution of that component 
to maintain ite-at’ constant concentration: 

All of the solutions, with the exception of the Zinc, 
cadmium and lead, were prepared and standardized on the 
basis of weight. The concentrations Chr Cu [Lbs 3 Le) 


and so on, discussed in later sections, are reported in 


terms of moles/kg of solution. 
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Boot rtLation Apparatus. 


The FICS method requires a considerable quancityveot 
haghniy reliable data, In addition, to satisfy the 


theoretical requirements of the method, C,, must be constant 


M 
while evaluating SS SHgvan ene isis and Cc, must be constant 
while calculating et OOM! Ce epEr This has been done by 


adding appropriate quantities of concentrated components 
M and L to compensate for dilution caused by the titrant. 
The system designed and constructed for this thesis 
work is shown schematically in Figure 11. It included a 
PDP-11/10 computer, a Fisher 520 digital pH meter, a 
Sartorius 3015 electronic balance and its related titrant 
delivery unit, and an auxiliary delivery unit composed of a 
Mettler DVli,.digital, piston drive: andi ml\buret fox 
adding concentrated metal or ligand and a multiplexer. The 
PDP-11/10 computer had cassette tape and floppy disc options 
for storage of data and a Laboratory Peripheral System 
(LPS) option to facilitate input/output (1/0)operations. The 
Parts of the LPS used in this work were the LPSDR option 
consisting of a sixteen bit buffered digital I/O, and the 
LPSKW option including a real time clock and Schmitt 
Ehlggers. 
Titrations were performed under computer control, and 
included the following operations: delivery of titrant from 
a gravimetric buret located on the pan of the SAT COLLIS 930.5 


balance, delivery of concentrated reagents from the 
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METTLER DV11 auxiliary buret, sampling the digital outputs 
of the pH meter and electronic balance through the multi- 
plexer, and manipulation and storage of the data for 
further use with the FICS method. A titration suitable for 
the FICS method would include a sequence of operations 
beginning with weighing of the unit on the balance, delivery 
of a portion of titrant, and reweighing. The weight of 
titrant delivered was used to calculate the amount of 
reagent to be added from the auxiliary buret, and the 
delivery was then executed. The electrode potential was 
measured and tested for equilibrium. When the electrodes 
had reached equilibrium, the potential was recorded and 
converted to pH on the basis of a calibration with NBS 
standards. This completed the operations required to 
collect a Single data point, and the process was repeated 
every 0.1 pH unit to the end of the titration. Each data 
point On cern of the following: informations »tiearora 
weight of titrant delivered to that point, the total weight 
of reagents added through the auxiliary delivery unite, sche 
equilibrium electrode potential after those additions, the 
pH calculated from the electrode potential and the total 
length of time since the titration began. This information 
was printed out for permanent record and stored on floppy 
disc for future manipulation by the FICS method. 

Complete descriptions of the gravimetric Citracaon 
system and details of the software control of each of the 


operations described above appears in. Chapter IV. (In the 
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next section the method used to perform the titrations 


designed for the FICS calculations is described. 


C. Method. 


Ghe Tierations were pertormed at. 25.0 +70. 1) °Cyunder 
an inert atmosphere of CO, and 0. DEecwnrerogenmorsangon:. 
The room temperature was maintained at 23 + 1 °C and the 
titrations were performed in a room having constant 
artificial illumination. (The titration procedure consisted 
of delivering 0.3 moles/kg NaOH titrant to 200 gm of sample 
solution with a gravimetric buret. The pH of the solution 
was measured with a Corning glass electrode and a Ag/AgCl, 
Na,SO, (saturated) reference electrode. Addition of 
concentrated reagents, as required by the FICS method, was 
Carried outiwith wa,digital piston bunet. 


The glass electrode was standardized in terms of 


hydrogen ion activity with phthalate (pH 4.008), phosphate 


(pH 7.413) @md carbonate (pH 10.012) NBS reference’ standards. 


This consisted of measuring the electrode potential in each 
solution and calculating a least squares linear regression 
between pH and electrode potential (mV). (See program 
BUFFER, Appendix C.) The electrode potential at each 


experimental point was converted to pH with the relation 


pH = slope(mv) + intercept CL9) 
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where the Slope and intercept See eae by the linear 
least squares fit to the NBS standards. 

The solvent containing a slight excess of nitric acid 
Served as a secondary pH standard. The hydrogen ion 
activity of this solution was measured in terms of the NBS 
Standards and was used as a reference starting point for 
all of the metal-ligand titrations. A weighed portion of 
the solvent was added to a titration cell and brought to 
25.0 °C, and the electrode potential was measured. The 
other components were then added by weight to make up the 
total weight of 200 gm. This eliminated electrode 
handling between the standardization and the beginning of 
each titration, an especially important consideration when 
sang the FiCS method. “As descrabed in, Chapter, 11, the pH 
serves as a reference scale for the comparison of Ci from 
eitration. ter,titratdon.. Lb.as essentralethar thas) scale 
be ene tee especially at low pH. The solvent,containing 


about 1072 moles/kg of HNO made this possible. The pH 


3" 
of this solvent was assumed to be constant and provided a 
fixed initial point from which to start each titration, 
without any electrode handling between standardization with 
the solvent and addition of the components of the sample. 
The solvent was, therefore, in effect simultaneously a 
single point pH standard and a part of the sample solution. 
As discussed in part C, Chapter II, as much care as 


possible must be taken to avoid systematic error in the 


preparation of sample solutions. Since it seemed likely 
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that errors would be largest in the first solution prepared, 
the order of the titrations was chosen so that these errors 
would have a minimum effect on (dc,,/dc, ) and (dc, /dc,,) - 

The middle titrations were performed first, and those at the 
extremes of Ch and Cu last. The usual order was C.{3), 

C, (2), C, (4), C,(1) and Cc, (5) where Ca Die G2 Ca) 


CL (4)< cc, (5). 


De. Activity Coefficients, 


The hydrogen ion activity served as a reference scale 
for use with the FICS method in the work described in this 
thesis. As described in the previous section, the pH of 
data points was calculated on the basis of the activity 
values of NBS standard solutions. In this way, all 
experimental values of [L] and [M] were measured at pH 
values. The stablity constants reported in this work are 
therefore mixed concentration-activity constants. Because 
the mass balance for Ch was not used in the calculation of 
formation constants, the hydrogen ion concentration was not 
required at any stage of the FICS calculations. However, 
the use of the sample diluent, or solvent, described in the 
method section, provides a convenient estimate of the 
hydrogen ion activity coefficient, (ha 

Theractaivaty coeriicient Vopaea a solution of 
0.3 moles/kg NaClO, was estimated using the measured a,,+ 


4 
based on NBS reference standards and the known concentration 
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of nitric acid, determined by Ceeracion with NaOH. A 
Solvent with an acid concentration of 1.16 x TOs moles/kg 
had a measured pH OF 4.006, cglving a ies OLm Ooo. FAs 
repeat of this procedure with a second solvent yielded an 
estimated Vat Or 0.91) Bothyorethese values are wider 
than an expected {+ of 0.7 in (0.3: mohes/ kg Naclo, (58). 
The formation constants reported in this work may be 
converted to concentration constants by converting the 
hydrogen activity terms to hydrogen ion concentrations in 
the species that include protons. This step requires no 


modification in any parts of the FICS calculations. 
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CHAPTER IV 


THE GRAVIMETRIC TITRATION SYSTEM 


Ree ntroauction . 


A titration system was designed for the specific 
purpose of collecting data suitable for the FICS evaluation 
Procedure. Based on the discussions in Chapter II, such a 
system must be accurate and must be able to maintain one of 
the components in the sample at constant concentration. 

From a practical point of view, the system must also be fast 
and reliable. These requirements have been met by a 
titration system shown schematically in Figure 1l. It 
consists of a digital pH meter, an electronic balance and 
the associated titrant delivery system, an auxiliary metal 
or ligand delivery system, a multiplexer, and a PDP-11/10 
computer. The computer controlled the delivery of reagents 
from each of the titration units and, through the multi- 
plexer, collected electrode potential and weight information 
from the pH meter and the electronic balance. This data 
was stored for future use with the FICS data handling 
method. The auxiliary delivery system, under computer 
instructions, maintained the metal or ligand concentration 
constant throughout the titration by compensating for 


Ride On sy, ie cant. 


Two designs of gravimetric delivery units were built 
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and tested. The first was based on a simple mechanical 
couple which opened and closed a pinch valve that controlled 
the titrant delivery, and the second employed an optical 
couple in which a beam of light triggered the opening and 
flossing of theiwvalve!that.controlled the!’ £low of titrant. 

In both cases, the system was designed so that the weight 
of the delivered titrant was known. 

Each part of the automatic gravimetric titration 
system will be considered in detail and the results of its 
application to simple acid-base titrations discussed in 
terms of the accuracy, precision, speed and reliability 


of the system. 


B. Instrumentation. 
The pH Measurement 


The pH of the sample solution was measured with a 
Corning glass electrode and a Ag/AgCl reference electrode. 
The potential between these electrodes was measured by a 
Fisher 520 digital pH meter and the value of that potential, 
invmiliivolts (mV), was fed as a ‘parallel binary coded 
decimal (BCD) signal to the multiplexer. This value was 
sampled by the PDP-11/10 under program control. 

After each delivery of titrant a series of electrode 
potential measurements were taken and the average, the 


standard deviation and the linear regression of the points 
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calculated. If the standard deviation and the drift in 
electrode potential were within the limits defined as 
equilibrium, the average of the set of measurements was 
taken as the electrode potential for that data point. This 
potential was converted to pH on the basis of electrode 


standardization with NBS buffers. 
A Titrant Delivery System Based on a Mechanical Couple 


A schematic diagram of the mechanical couple delivery 
system is shown in Figure 12. The titrant was delivered 
by gravity from a reservoir located on the balance pan of 
a Sartorius 3015 electronic balance, through:a normally 
closed pinch valve, and to the sample cell. 

The system consisted of four main components. The 
first was a polyethylene reservoir having an opening in the 
bottom oe delivery of titrant and an opening at the top 
leading to a CO. trap. The titrant was delivered to the 
sample cell through Teflon tubing. 

The second component of the system was a small, 
normally closed pinch valve. The mechanical energy necessary 
to open it was supplied by a motor driven device not 
actually sitting on the balance pan. With this design, a 
mechanical couple existed to the balance pan during the 
process of titrant delivery but not during the weighing. 
The mechanical couple provided the means for simultaneously 


transferring the information to start and stop the titrant 
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flow, and for transferring the mechanical energy necessary 
to open the valve. In the next section a device based on an 
optical couple will be discussed in which the information, 
but not the energy necessary to operate the valve was sent 
to the unit on the balance pan. 

The third major component of the mechanical couple 
System was the motorized drive. This consisted of a small 
1.5 V motor which drove a gear mechanism. The device was 
designed so that an arm pulled down to initiate titrant 
delivery and returned up when the motor was reversed. In 
the up position no contact was made with any part actually 
on the balance pan. When the motor pulled the arm down 
contact was made with the pinch valve, opening it and 
allowing the titrant to flow. 

The power to the motor was supplied under program 
Sewtrols..Ausimple Circuit Involving, two, relays, anda 
battery allowed two bits of digital information to drive 
Hae motor .in eather, direction .or to stop it... The binary 
code 0,0 and 1,1 supplied no voltage to the motor, while 
Orleranathewmotor sin Ones di uectionpandy,.0uran Ue einetie 
other. Delivery of titrant occurred through the following 
Series of events : (a) 0,1 began the opening of the valve, 
(b) once the valve was open a microswitch signalled that the 
motor should be stopped, and the bits set to 0,0, (c) there 
was a delay while the titrant was delivered, (d) 1,0 fora 
period of time allowed the valve to close, and (e) 0,0 


stopped the motor. At this point the routines designed to 
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sample electrode potentials and balance reading were set 
into operation. 

The time necessary for a delivery of titrant included 
delays at steps (a) and (d) during which the motor opened 
and closed the valve and a delay in step (c) while the 
titrant flowed. Steps (a) and (d) required 4.5 seconds 
each, and the delay-in step (c) varied from several milli- 
seconds to 20 seconds. The overall delivery time was a 
minimum of 9 seconds. The total time for each delivery in 
eieactual titration also included titrant flow time, about 
5 seconds for the balance pan to come to rest and a time 
for the electrodes to reach equilibrium. The total time 
between titrant additions then came to about 30 seconds. 

The final major component of the system was a housing 
Pitt around the titrator so that air currents would not 
affect the balance reading. Two housings were 
tested, The first covered all of the parts on the balance 
pan and, in addition, had a tube leading to the sample cell 
and coupling with it. This housing made the atmosphere 
over the sample cell continuous with that over the balance 
pan and nitrogen could be pumped over both simultaneously. 
The second housing, and the preferred one, also covered 
the balance pan but did not lead to the sample cell. A 
hole was provided for the titrant delivery tube, but the 
atmosphere over the balance was completely independent 
from that over the sample cell. Several problems were 


encountered with the first design. The balance was very 
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sensitive to changes in humidity and temperature caused by 
the flowing nitrogen and, in addition, some buoyancy effects 
caused by the change in atmospheric components were observed. 

The procedure used for the collection of data from the 
gravimetric delivery system was very similar to that for the 
Piemeler. The BCD value of the weight .of the unit on) the 
balance pan was fed to the multiplexer. This value was 
sampled under computer control and the average of the series 
of weights taken as the best measure of the balance reading. 
Unlike the treatment of values from the pH meter, a standard 
deviation and slope were not taken. It was therefore 
necessary to delay the sampling a sufficient period of time 


to allow the balance pan to come to rest. 


A Titrant Delivery System Based on an Optical Couple 


The optical couple gravimetric delivery system is 
shown schematically in Figure 13. The titrant was delivered 
from a spring loaded syringe, through a solenoid valve, and 
EOrthe titration -cell. 

The optical couple system consisted of four major 
components. The first was the titrant reservoir and 
delivery system. This consisted of a 10 ml Hamilton Gas- 
Tight syringe with a plunger modified to make it lighter,and 
Teflon tubing connected to the syringe by a tight Luer- 
type fitting and to the solenoid valve by regular high 


pressure liquid chromatography forcings. The: tubing 
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leading into the titration cell was drawn to a fine jes go) 
so that titrant was delivered as a tiny jet of Liquid. 
This eliminated the problem of evaporation from a hanging 
drop, which may be serious with volatile solvents. 

The second major component of the Be was the valve 
and its power supply. A General ValveCorporation 2 way, 
12 V DC, Iso-latch miniature solenoid valve was used to 
control the solution flow. The valve was opened by a single 
12 V pulse and closed by another 12 V pulse. Current drain 
during each was about 700 ma, lasting approximately 150 msec. 
The power supply was a series of ten Eveready B225T, 1.25 V 
NiCd batteries. These batteries formed most of the weight 
Or chesdevicetsitting, onthe balance: pans 

The third major component of the system was the 
ancormation transfer unit, in this case an optical couple. 
A Single bit of digital information was used to control a 
meant. © Thuasear0vsignalled the laghtrtoscurnyonsandgayd 
Horned. it off. A. phototransistorvilocatedson they balance 
acted as a switch and controlled power to the solenoid 
waive. Therefore, “the anstructions toxopen and close the 
valve were supplied by the optical couple and the mechanical 
energy was supplied from the battery. This system had two 
major advantages; first, there was no mechanical interaction 
with the balance pan, making possible the simultaneous 
delivery of titrant and weight measurement and, second, the 
system was fast. The mechanical couple system required 4.5 


seconds for each of the valve opening and valve closing 
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steps, while each pulse to the optical couple required 150 
msec. In addition, the mechanical couple system required 
approximately 5 seconds for the balance pan to come to rest, 
and less than one second for the optical system. Each 
titrant addition,therefore in theory, required 13 seconds 
less for the optical than the mechanical couple system 
fapact from actual titrant flow time, etc.) 

The final component of the optical coupled system was 
the housing over the balance. This minimized drift due to 
err Currents, and provided a convenient location to hold the 
light source as close as possible to the phototransistor. 

As discussed earlier, the preferred housing did not couple 
Marectly to.the titration cell. A streamuof N, could be 
blown over the sample solution, but not directed into the 
balance. 

The collection of data from the optically coupled 
titrator was identical to that described for the mechanical 
coupled titrator. The only difference) was “that a shorter 
delay time was needed between the end of delivery and the 
sampling of the electronic balance output. 

Fach of the titrator designs had certain advantages 
and disadvantages. The optical system was capable of 
delivering titrant within short periods of time, saving 
about 13 seconds per delivery compared to the mechanical 
system. In fact however, not all of this time saving can be 


realized during a titration. The time of mixing of the 


sample solution and the rate of equilibration of the 
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electrodes must also be taken into consideration. The rate 
at which the two systems are capable of Operating is compared 
in more detail in the section dealing with the application 

of these units to simple acid-base titrations. 

Another important difference between the titrators was 
the quantity of reagent held by each. The solenoid valve 
and the battery of the optical couple system were heavy, and 
aoma result only 10 om of titrant could be held in the 
Syringe, suitable for one, or at the most two titrations 
(5-10 gm each). The mechanical system was simple and light, 
and held over 100 gm of titrant, enough to complete over ten 
titrations without refilling the reservoir. 

The quantity of titrant delivered by both of these 
systems was controlled by the length of time the valve was 
held open, measured in 0.01 second intervals. The number 
Ofesuch intervals, or the time which the’ valve was held 
open, was manipulated from data point to data point according 
roethe effect of the =previovws addition of titrant on fhe 
electrode potential. If the potential changed by less than 
a"minimum"value, then the next delivery would be twice as 
long. If the change was greater than a"maximum"limit, the 
delivery was halved, and if the change was greater than a 
"dangerous"level,then the delivery time was cut by five. 
With appropriately chosen levels, the systems could reliably 
handle even the abrupt end point of a titration of 4 moles/ 
kg HNO, with 3 moles/kg NaOH. The systems did not once 


3 
fail to get sufficient data to graphically determine the 
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equivalence point. 


The Multiplexer 


A computer-experiment interface is designed to 
coordinate the transfer of digital information from one 
device to another. For example, a short dialogue might 
occur between the computer and the device supplying 
experimental data; first, the computer requests the data; 
second, the device prepares the information, perhaps doing 
an A/D conversion; third, the device sends a data ready 
Signal; fourth, the computer reads in the data; and finally, 
the computer sends out a data accept signal. The role which 
the interface may play in this dialogue is the translation 
of the signals from one device into signals compatible with 
the next. This often means delaying of pulses, and changing 
of voltage levels or polarities. 

The multiplexer used in the titration system served as 
the interface between the digital outputs of the balance 
and’ pH meter“and the input’ of the LPS’ unit of the PDP=1%: 
The circuitry of the multiplexer will be considered first 
and then the manner in which the sampling is done under 
program control. 

The multiplexer made possible the transfer of data 
from each of several devices into a single input of the 
computer. The device employed for the titration system used 


16 parallel 4 to 1 multiplexers, one for each bit of a 4 
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digit BCD number. Bit 00 from each of 4 different devices 
was fed to the first multiplexer, bit 01 from each into the 
next and so on. Two bits of information from the digital 
output of the LPS was used to select the device whose output 
was tO be-sampled. For example, a.0,0 signalled bit 00 from 
device 1 to appear on the output of the first /multiplexer, 
bit 01 from device 1 to appear on the next output and so on. 
The outputs of the 16 channels were equivalent to the output 
Of device 1. If the computer signalled 0,1, then the outputs 
of the 16 channels equalled the output of the second device, 
and so on. This system allowed sequential sampling of four 
different sets of BCD data. In this way it was possible to 
sample the output of the pH meter and the electronic balance. 
Due to the number of digits produced by the balance, its 
output was treated as the equivalent to two devices. The 
final channel was filled with zeros for debugging purposes. 
The digital input of the LPS unit was designed for two 
Sieferent modes of input. The first was’ the bit stimulus 
mode, in which the only bits recognized by the computer were 
those eee risa the 0 to 1 transition during the sampled 
period. This could be accomplished with the multiplexer 
by feeding the channel with all zeros to the LPS, then 
switching to the channels with the pH meter or balance. 
However, this turned out to be unreliable. To avoid this 
problem, the word stimulus mode was used, allowing the entire 
BCD code to be read directly. Unfortunately, the hardware 


of the LPS was not compatible with the aquisition routines 
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of the language BASIC. For this reason the program running 
the titration, called TTR1 (see Appendix D), was written in 
FORTRAN. Because of its great flexibility, especially in 
the debugging and editing of programs, all other programs 
were written in BASIC. The data collected in FORTRAN was 
converted to a form compatible with BASIC with a short BASIC 


program called FCOPY (see Appendix E). 


C. Software Control of the Automatic Weight Titrator System. 


The titrator system operated entirely under the control 
of a single FORTRAN program called TTR1. While this program 
was being executed, the following operations would occur: 

(a) sampling of electrode potential and electronic balance 
Output through the multiplexer, (b) anticipation of quantity 
Steritrant anarther delivery, procedure, ({c). calculation’ o£ 
quantities of solution to be added through auxiliary delivery 
systems (Mettler DV11) and subsequent delivery and, (d) the 
algebraic manipulation of the data (e.g. conversion of mV 

to pH) and the storage of that data on magnetic tape or 
floppy disc. Before considering the most important 
Operations separately, an overview of the sequencing and 
timing of the entire program is essential. This is best 
described in flow chart form, Figure 14. 

The initialization of variable parameters allowed the 
Operator to fix the parameters at values suitable for the 


particular titration. These parameters included; (a) the 
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Figure 14. FLOW CHART OF PROGRAM TTRI 
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Slope and intercept determined by the calibration of the pH 
electrode (see program BUFFER, Appendix C) used to convert 
electrode potentials to pH, (b) the duration of the first 
delivery of titrant, Ce) the"minimum; "maximum"and"dangerous" 
changes in the electrode potentials used to adjust the 
quantity of titrant delivered from point to point during 

the titration, (d) the number of readings of electrode 
potential to be taken in each sampling sequence, (e) the 
maximum allowed standard deviation and drift in electrode 
potential of the series sampled in (d), (f) the electrode 
potential which signalled the end of the titration, (g) the 
number of auxiliary burets to be operated (zero, one or two) 
and, (h) the concentrations of reagents in the auxiliary 
burets and the concentrations of the corresponding components 
in the sample cell, quantities which are needed in order to 
compensate for dilution. 

With the values of these parameters fixed, TTR1 sampled 
the outputs of the pH meter and electronic balance and added 
the first titrant. The sampling was repeated, and the weight 
of delivered titrant calculated. If it was required by the 
PICS method, the quantity of reagent needed to compensate 
for dilution caused by this delivery of titrant was 
calculated ee the necessary reagent added through the 
auxiliary delivery system. After this delivery, the sampling 
procedure was again repeated to determine the electrode 
potential. The standard deviation and drift in the electrode 


potential was calculated and, 4fo the. Soluriom was. note 
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equilibrium as defined by the operator, the sampling process 
was repeated. In this way the system collected equilibrium 
data after known additions of titrant and other reagents. 
Moreover, the conditions defined as equilibrium were 
identical from point to point, something not possible if the 
titration was run manually. 

If the electrodes had reached equilibrium, TTR1 
calculated the amount of titrant to be added in the next 
delivery, based on the effects of the previous addition. 

If the electrode potential changed by less than a"minimum" 
value, usually 4 or 5 mV, then the quantity of titrant was 
doubled; if the electrode potential had changed more than a 
"maximum"value,8 or 9 mV, then the quantity of titrant was 
cut in half; and if the change was over a"dangerous"limit, 
Bout t2 my, che, quantity Of titrant was cut by «five. Jeihe 
delivery was then made by sending the appropriate signals 
to the gravimetric delivery system. After the delivery the 
weight and electrode potential measurements were carried 
Out once again. This entire process continued until the 
electrode potential passed ‘thes limit defined to bem newend 
Sreathe ‘titration. 

The remainder of this discussion will deal in detail 
only with those operations unique to the present system. 
That part of the code which is familiar to, or understood 
by the FORTRAN programmer will be omitted. 

The program TTR1 communicated with external devices 


through the LPSDR option (sixteen bit buffered digital I/O) 
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of the Laboratory Peripheral System. All of the program 
timing required the programmable real time clock paricor 
the LPSKW option. the versionso£ TTR which ran the 
mechanical couple titrator also required a Schmitt trigger, 


a part of the LPSKW option. 


Subroutine DELAY (TIME) 


Those parts of TTR1 requiring delays called on sub- 
routine DELAY(TIME), where TIME was an integer number of 
0.01 second intervals. A delivery of 10 seconds then 
required a call of DELAY(1000). The FORTRAN statements 
were: 
DO. 10. 1=1, TIME 
ICMF = 0 
CALL SETR (570,1.;,ICMF) 
CALL LWAIT (ICMF,0) 
CALL SETR(-1,,,) 

10 CONTINUE 
RETURN 


The first SETR statement set the programmable real time 


clock running at 100 Hz. The clock overflowed after 0.01 
seconds and set ICMF equal to 1. The LWAIT statement made 
the computer wait until the value of ICMF no longer equalled 


0. The second SETR statement stopped the clock. The DO loop 


then repeated this process for TIME internvalseor-Os01 


seconds. 
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Subroutine COLLECT (A,B) 


Subroutine COLLECT(A,B) sampled the digital lLoutputs) of 
the pH meter and the electronic balance. It returned the 
electrode potential in A and weight reading in’B. Since the 
output of the electronic balance was six BCD digits, two 
separate samplings were necessary, with algebraic combination 
Of those to give a single value for B. 

CALE DORI (070, 5 3510) 
CALL DELAY (1) 

Ciel DIR (0, 07=170) 
CAL aeEDORS (0+ 07 83i."1) 
CALL DELAY (1) 


B IDTR (0,0 ,-1,0) 


il 


B (10000*B)+C 
CALL apo (0:0, 3.,.°2) 
CALL DELAY (1) 
Aba eae. 0. 09—1 0) 
RETURN 
The first IDOR statement set the final 2 bits of the digital 


Subputeregister to, 0,0.. These bits were the channel*®control 


of the multiplexer. The channel signalled by 0,0 was the 
last four digits of the electronic balance, which were read 
in by the IDIR statement and stored as the variable C. This 
process was repeated for the first digits from the electronic 
balance, stored as variable B, and for the output of the pH 


meter, stored as A. Values of B and C were combined to a 
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Single value in B, including all the six digits from the 
balance. The statements placing the decimal points in the 
correct place,and so on, have been omitted for simplicity. 
The DELAY statements allowed the multiplexer sufficient 


time to change channels. 


Subroutine DELIVER(T) 


The delivery of titrant from the weight titrator was 
controlled by DELIVER(I). The duration of the delivery, in 
0.01 second units, was specified by the variable I. Delivery 
using the mechanical couple or the optical couple titrators 
were almost identical and, for simplicity, the optical system 
will be considered first. 

CALL stor. (0,0,- 4,4) 
CALL DELAY (15) 

CALL TDOR (0,0, "45 "0) 
CALL DELAY (I) 

CALE IDOR) (0,0; 4"4) 
CALL DELAY (15) 
GAREeTDOR? (0,0, "4, "0) 


RETURN 
The IDOR statements controlled a single bit of the digital 


output register. The first IDOR set the bit which turned on 


the light to open the valve. The light remained on for 150 


msec and was turned off. A delay of I times 0.01 seconds 


occurred while the titrant flowed. A second flash of light 
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of 150 msec duration closed the valve. 

The mechanical couple system required that the IDOR 
statements control two bits. The first IDOR set the motor 
moving down, the second stopped it, the third set the motor 
am reverse and the. final stopped it. The first DELAY was 
designed not to wait for the clock overflow, but for a 
Signal from the Schmitt trigger. The second DELAY was 
exactly as above and the third was fixed at 450 rather than 


15 as above. 
Subroutine ADELIVER(N) and BDELIVER(N) 


These routines were designed to deliver a series of 
pulses to the Mettler DV11 automatic piston buret. The 
FORTRAN code was: 

DO. VOT y=. 1 ZN 

CALEUTDOR: (0,0,720,) 20) 

GALL. DEGAY, (5) 

CALL LOR. 0,0, 2 2040 109 

CALL DELAY (5) 

POaecONTINUE 

RETURN 
This DO loop produced a square wave with a period of 100 
msec. The bit was repeatedly set and removed by the IDOR 
statements, with 50 msec delays between each. The sub- 
routines differed only in the bit which they controlled. In 


this way, different volumes could be added from separate 
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burets, as might be required for certain types of titration 


data. 


D. Application of the Automatic Weight Titration Apparatus 


to Simple Acid-Base Determinations. 


The reliability and the precision of the mechanical and 
the optical couple systems were evaluated using simple acid- 
base determinations. This work was designed to find those 
conditions, both from a chemical and from a software point 
of view, which would give the most valuable data. Since 
the immediate application of the system was the collection 
of data suitable for the FICS method, only discrete, 
equilibrium data point collection was studied. The apparatus 
has great potential for use in a variety of modes, including 
a mixture of continuous and discrete data point collection, 
PH stat type applications, fixed end point titrations and 
sO on; however, these have not been considered. The 
apparatus is also readily adaptable to other methods of 
sample monitoring, for example spectroscopic, where digital 
information is accessible. 

Those parameters which are of most importance in 
evaluating an automatic titration system include the rate 
OE titration, the reliability of the entire system and, the 
accuracy and precision of the resulting data. “Titratryvonseor 
potassium hydrogen phthalate with sodium hydroxide and 


titrations of nitric acid solutions with sodium hydroxide 
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were carried out to evaluate the performance of both designs 


of automatic weight titrators. 
Results and Discussion 


Consideration will first be given to the precision 
attained by the two basic designs of titrators, followed by 
a comparison of the different systems to illustrate their 
respective advantages and disadvantages. 

Tables 23 and 24 are summaries of titrations of KHP and 
HNO , with NaOH using the mechanical couple system. The 
relative standard deviation of the KHP series was 0.06%, or 
mbout: 4 mg of titrant (based on’ 7 gm total. delivery). The 
relative standard deviation decreased to.0.04%, or about 1 
mg for the strong acid-strong base titrations summarized in 
maple 24. Tables 25 and 26 are similar sets ofsdata 
Pollected using the optical couple titrator~ The standard 
deviations are essentially unchanged, about 4 mg of titrant 
fecetha KHP and’ about 2 mg of titrant for the HNO. titrations. 
Since the precision of the balance is 1 mg, this may be 
expected to be the limit of precision imposed on any set of 
titrations and this appears to be the case for strong acid- 
strong base titrations. The somewhat lower precision of the 
KHP titrations is due to a combination of factors including 
the lack of a sharp endpoint in the titration of KHP and 
the susceptibility of these titrations to errors introduced 


by impurities or co. in the reagents since the equivalence 
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Table 23. 


Sample 


NO —&F WN HF 


ey, 


Precision of the Standardization of Sodium 
Hydroxide with Potassium Hydrogen Phthalate 


Using the Weight Titrator Based on a Mechanical 


Couple. 
Consumption Calculated 
Weight KHP, of NaOH, NaOH Concentration 

oe gm moles/kg 
0.4994 65657 On3673 
0.4996 6.656 043675 
0.4856 6.464 0.3678 
025056 6.733 On 3677 
Oooo 7.830 0.3674 
OO125 oe Sy 02.4608 
Average 0.3676 


Relative Standard Deviation 


0.06% 
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Table 24. Precision of the Strong Acid-Strong Base 
Titration Using the Mechanical Couple Weight 


Tactratcor ! 
Consumption Calculated 
Sample Weight HNO,, of NaOH? HNO, Concentration 
gm gm moles/kg 
1 7.9476 35500 2.3242 
2 W6ST2, | cara 1 93247 
3 5152433 2.314 1.3240 
4 5.4243 C302 3229 
5 TAOS 5 PS he, 3256 
6 ye a AS Bea Ga kee Be, Laos 0 
if 6.6570 293) M5230 
8 Tage p 2 ay PCPs | 


Average 1.3236 


Relative Standard Deviation 0.04% 


“concentration of NaOH was 3.0 moles/kg. 
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Mable 25, 


Sample 


Oo ON HD HO F&F W ND FEF 


Precision of the Standardization of Sodium 
Hydroxide with Potassium Hydrogen Phthalate 
Using the Weight Titrator Based on an Optical 
Couple. 
Consumption Calculated 
Weight KHP, of NaOH, NaOH Concentration 
gm gm moles/kg 
epee i babe, 6.427 Oi59 20 
O75 257 PPRAE STS) 053-929 
0.6801 8.482 023926 
0.6739 8.405 0.3926 
Oita el Co) 6.452 0,..39:30 
0.6720 9.2348 Ose 
0.4927 6.140 O13929 
0.4990 6.214 03932 
0.5074 Swe fie 053962 
Average 0.3929 


Relative Standard Deviation 


0.06% 
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Table 26. Precision of the Strong Acid-Strong Base 
Titration Using the Optical Couple Weight 
TLitratort 


" Consumptign Calculated 
Sample Weight HN 37 of NaOH, HNO. Concentration 
ER eee ens iy eo RUCTI Oras. ee 5) kate moles/kg 
if 62.7709 9.962 4.4139 
Z 6.6627 9.806 4.4153 
3 72938 LOnS89 4.4146 
4 Gal 30 92938 4.4149 
5 7e0ds6'6 10.328 4.4158 
6 6.3303 10.064 Ae o: 
Zz 6.8340 EO S07 4.4148 
8 ING es OM 10... 330 4.4153 
9 Las Z5 TORS Gu 4.4168 


Average 4.4152 


Relative Standard Deviation 0.018% 


@concentration of NaOH was 3.0 moles/kg. 
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point is above pH 7.0. Error may also have been introduced 
in the sample preparation stage requiring the weighing of 
small quantities of solid KHP. 

The titrators were equivalent from the point of view of 
accuracy and precision, both devices being limited by the 
same factors, namely either the chemistry of the system, or 
by the sensitivity of the balance. The mechanical couple 
system required a minimum of 9 seconds to deliver titrant, 
plus 5 seconds for the balance pan to come to rest and 6 
seconds to sample the electrode potential. This resulted 
in a basic delivery time of 20 seconds. The optical couple 
required only 7 seconds for the same operations, saving 13 
seconds. These were not the only procedures which required 
considerable time however. The titrant delivery consumed a 
large amount of time. In order that a minimum delivery of 
5 mg could be attained, the same flow rate delivered 0.5 gm 
in 30 seconds. Both the optical and mechanical couple 
devices were limited to delivery at that rate. A,vsignificant 
amount of time was also required for the electrodes to reach 
equilibrium near the equivalence point. This was in part 
due to the rate of mixing of reagents and could be improved 
with better cell design. 

The delay during which the electrodes came to 


equilibrium removed much of the speed advantage of the 


optical couple titrator. After the addition#ot tituant both 


systems required several seconds of delay before sampling 


the electrode potential; however, because of its slowness 
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the mechanical couple titrator needed a slightly shorter 

delay. If the delay time used with the optical system was 
shortened, the electrodes would not be at equilibrium and 

the sampling step would have to be repeated, resulting ina 
longer effective delivery time than that of the mechanical 
couple. Since about 60 data points were collected for the 
FICS method, a few seconds per delivery could considerably 


change the total time required for the titration. 


jhe total time for’ titration using the optical ‘couple 
ranged from 10 minutes, collecting 30 data points suitable 
for equivalence point measurement, to about 40 minutes to 
collect 50 data points suitable for use with the FICS method 
which also included delivery from one extra buret. The 
mechanical couple system required from 20 minutes to 60 
minutes for the same experiments. The optical couple system 
required about 20 minutes to run each of the titrations shown 
in Tables 25 and 26 (35-40 data points). 

In addition to the actual running time of the titrator, 
time was consumed in manual operations. This included 
removing the titrated sample, rinsing the electrodes and 
reassembling the apparatus with a new sample. The optical 
couple titrator also required that the syringe be refilled 
between experiments. Futher refinements of this system 
should allow, as did the mechanical couple titrator, multiple 
titrations with a single reservoir of reagent. This is 


especially true if the apparatus is to be applied to routine 


gravimetric reagent handling. The total time necessary to 
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perform the nine titrations in Table 26 was 3.5 hours. This 
is slow compared to manual, colourimetric end point 
titrations with a buret, about equal in time for the same 
titrations using an automatic piston buret, but much more 
rapid than the equivalent manual gravimetric titrations. 
This increase in rate over manual gravimetric titrations 

was attained with no sacrifice in precision. 

Both titrators were found to be very reliable, both 
from the point of view of mechanical and electrical operation 
and also in the ability of the entire unit to respond to 
very sharp equivalence points. For example in the titrations 
shown in Table 26 approximately 35 data points were collected, 
wMetel2 of those were within +0.4 gmiof the end point (total 
titrant delivery, 10 gm). This allowed graphical 
determination of the end point using a significant number of 
the total data points. The system was able to respond to 
the approach of the sharp change of potential without either 
extreme caution, that is, without giving an excessively 
large number of points, or uncontrollably large additions 
which would make graphical determination of the end point 
impossible. The system never failed to collect sufficient 
data to allow graphical end point determination. 

Several potential problems and limitations: of ihe 
system were encountered. These were involved primarily with 
the long term behavior of the electronic balance. It was 


found that a very considerable rate of drift could be induced 


by a change in humidity. An increase in humidity brought 
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about an apparent decrease in the weight reading of as much 
as 5 mg in 15 minutes, while a decrease in humidity brought 
an apparent increase in weight. The effect of temperature 
was not specifically studied but an is certainly expected 
to have similar effects. 

These problems introduced difficulty in the design of 
a controlled atmosphere titration cell. First, there must 
not be any mechanical connections between the delivery tube 
and the cell. This precludes a tight seal around the 
Opening where the titrant is delivered to the cell. The 
alternatives are therefore to either pump a sufficient 
quantity of inert gas over the cell to minimize entry of 
atmospheric components into the cell or to contain the 
titrant delivery system and the cell in a common housing. 
The former alternative was taken after work with a common 
housing revealed that very elaborate means would have to be 
taken to control the atmosphere sufficiently closely to stop 
balance drift. In certain applications, for example with 
volatile or reactive substances the latter method may be 
employed if the titration is carried out sufficiently 
rapidly that drift error is negligible for the desired 


precision of results. 


Summary 


The nature of the FICS data handling procedure has 


placed several requirements on a titration system. The 
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requirements for accuracy have been met by the design of 
two forms of automated gravimetric titrators and by the 
complete automation of data handling. The system was fast 
and reliable, characteristic of computer. controlled 
experiments, and made decisions concerning equilibrium of 
ene electrode potentials and quantities of material to be 
delivered by all the burets with a speed and reliability 
impossible manually. 

The automated gravimetric titrator systems based on 
optical and mechanical couples were tested on simple acid- 
base titrations. From these experiments optimum values of 
delay intervals, criteria for equilibria, and other important 
parameters were evaluated. Several important instrumental 
problems including the design of the housing around the 
titrator were also resolved. The system based on the 
optical couple performed better than did the mechanical 
couple, especially in terms of speed, and was used in all 
later experiments. 

In the next chapter a study of the complexation of zinc 
with aspartic acid is described. The data was collected 
using the optical couple titrator system and the treatment 


of data was carried out using the FICS method described 


in Chapter II. 
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CHAPTER V 


APPLICATION OF THE FICS METHOD TO THE DETERMINATION OF THE 
STABILITY CONSTANTS OF COMPLEXES BETWEEN ASPARTIC ACID AND 
ZINC | 


Aeeantroduction. 


Aspartic acid is a potentially tridentate ligand with 


the structure shown below. 


0 Q 
HO-C—CH ,—CH-C—OH 
nat 
3 


Binding to metal ions may occur through the 3 carboxyl 
group, and the C& carboxyl and CA amino groups. At low 
PH all of the groups are protonated, and at pH values near 
PeaAsands9, the Cc carboxyl, B carboxyl, and amino protons, 
Bespectively, are removed., To besmore precise) this acid 
dissociation scheme must be described in terms of a series 
of microscopic acid dissociations. For aspartic acid the 


first proton removed may be either the Wor iS carboxyl 


protons. The structures of these are shown below. 


9 
fe) * 1 oe. 
—~CH-C—OH HO-C-CH—CH-€-0 
| + 
NHS NH, 
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Both of these species have the stoichiometry HoL. Removal 


of a second proton will lead to a single form of HL. 


ie 0 
o—t—cH 9—CH-C~07 
| +t 


NH, 


While there are three possible forms of HL, due to the low 
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acidity of the amino group this form will by far predominate. 


In a solution of metal ions, the metal ion may dis- 
Peace Carboxyl or amino protons and bind to those groups. 
If all three protons have been displaced the molecule may 
easily bind to the metal through the three groups simul- 
taneously. The structure of aspartic acid is such that it 


very neatly surrounds the metal on three adjacent sides. 


Two ligands may bind simultaneously, completing the 
octahedral coordination sphere preferred by a number of 
metal ions. . The two strongest complexes of aspartic acid 
with most metal ions are therefore ML and ML. . 

The zinc-aspartic acid system was chosen for the 
study and application of several of the procedures developed 
as part of the FICS method. A relatively simple chemical 
system was necessary in order to avoid, as much as possible, 
difficulty in assigning a model to the behavior of the 


system. The information this study was intended to provide 
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included: (a) the accuracy of experimental results which 
could be obtained with the system described in Chapter IV, 
(b) the severity of curvature of the (dc,,/dc, ) and (dc,,/dc,,) 
functions at various ligand to metal ratios, (c) the 
applicability of the methods previously discussed with 
respect to the determination of PIL], and PIM] o, and (d) 

the reliability of the constants derived with the available 
data collection methods, and the experimental modifications 
necessary for the study of the chemistry of more complex 


systems. 


B. Determination of the Acid Dissociation Constants of 


Aspartic Acid 

A series of solutions containing from 2.643 x Ton 
to 5.908 "x 108: moles/kg of aspartic acid were titrated 
meen NaOH from pH 3.5 to Ti.0°” Ehe*rilCS method can moose 
eacily be applied to the titration of an acidic” ligand 
Since: (a) no acid needs to be added to the solution to 
Siang 1t to PHos making the solution composition data more 
accurate, (b) the concentration of the acid in each sample 
can be determined from the equivalence point following 
the consumption of the carboxyl proton, also avoiding the 
errors normally expected in the solution preparation, and 
(c) the range of C, can be unlimited since (dC,,/dC,) is 
linear. This final factor is of considerable importance 


since the precision of (dC,,/dC,) or alternatively 4C,/4C) 
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is very much dependent on the magnitudes of ac. and 4C,. 
Both must have good accuracy before (dc,,/dc,) is known 
with some reliability. 

The raw titration data was handled in exactly the 
manner described in Chapter II, and involved the calculation 


Of the total titratable proton concentration, C at each 


jab 
data point and interpolation between those points to deter- 
mine Cy at exact pH values. The raw data was collected at 


iieervals,of 0.1 to,0.2 pH units, and the Cu values 
required for evaluation of (dc,,/dc, ) were obtained by 
interpolation at exact 0.1 pH intervals. Figure 15 shows 
the variation in Cu with change in the total ligand 
concentration. As expected, the plots are linear and 
indicate that the experimental methods used can provide very 
accurate data since within the sensitivity of the graphical 
method, all of the data points fell along the lines. More 
exact mathematical treatment revealed that the fit of the 
Points to va straight line was less precise at extremes of pH. 
This is to be expected because at extremes of pH large 
quantities of titrant are required to change the pH slightly, 
Andi as al result, smallocerrorsan pHi produce increasingly 
large errors in Cir: This introduces uncertainty in the 
derivative (dc,,/dc, ) at pH values below 4 and above 10. 
Figure 16 summarizes the derivatives (dc,,/dc, ) asa 
function of pH. With those values the concentration of free 


aspartic acid, [L], was calculated throughout the pH range 


BRCheeavailabile data. The PIL], was first chosen assuming 
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mrgquyre 5"; Cy vs C, at various pH for Aspartic Acid. The 


values of the derivatives, (dc,,/dc,), are shown in 


Figure lo. 
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no protonation of the ligand at pH 11.0. The acid 
dissociation constants were calculated on the basis of that 
p{L], and from the new set of constants a new PIL], was 


Baiculated. This process required three stages, in which 


PIL] | was 3.000, 3.006 and 3.020, and the standard deviation 


132 


Se eecie fitsvaried as3%, 2.3¢%)and 1.83 srespectively. Using 


eepib) of 3.020,-a new p[L]_.of 3.020 was calculated. This 
re) fe) 


procedure worked as predicted by the studies using simulated 


Eitration curves. The .resulting protonation constants are 


enowh in Table 27. 


C. The Formation of Complexes between Zinc and Aspartic Acid. 


Several aspects of the FICS procedure discussed with 
respect to simulated titration curves were investigated 
using the zinc-aspartic acid system. It was of interest to 
determine the nature of the problems of curvature of the 
function of Cy with change in Ch for this metal-ligand 
System. A series of titrations was performed with the 
Soncentration Of zinc fixed at 9.370 x Ore moles/kg and 
the concentration of ligand varying from 1.1 x Mee to 


ie 0 aye moles/kg. The concentration of total titratable 


Proton, € was calculated; the variation of Cry at several 


H’ 
PH values is shown in Figure 17. There were no regions of 
sharp curvature, however, some curvature in the region of 


2:1 ligand to metal ratio developed in the pH range above 


7.0. This is expected because the complex ML, predominates 
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Table 27. Determination of the Protonation Constants of 


Aspartic Acid. 


Species tog ()/3)) 


pL] 5= 3.000 P{L] = 3.006 PIL] = o5020 
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in that pH range. At higher pH precipitation occurred in 
all of the titrations; the effects of this can be seen in 
Figure 17. The curve for pH 8.0 behaves as expected’ for a 
system containing the species ML and ML. (see Figure l, 
pH 5.0). However, at pH 9.0 the effect of precipitation 
at low Ch ends to turn othe Cu curve towards lower Cue 
This curvature becomes more severe as the pH is increased, 
and at pH 10.5 the curve has a pronounced downward curvature. 
These results are reasonable if the precipitation of metal 
hydroxide is considered to be equivalent to the formation 
of soluble hydroxy complexes which, as discussed in Chapter 
II, contribute a negative value to the concentration of 
micratable rotons. Ch therefore becomes more negative. In 


addition, it should be noted that the effect of precipitation 


on the Co vs C. curves begins at low Ch at pH 9° 0land 


L 


progresses towards higher’ C. as the pH is raised. (This/is 


L 
because precipitation of the metal is inhibited by the 
increasing concentration of ligand. On the basis of Figure 
17, the FICS method could most effectively be used above a 
3:1 ligand to metal ratio and up to a pH (arp Gages bette ys 

Ten titrations were performed; five at a constant Zinc 
Concentration of 9.370 x sly moles/kg and five at a constant 
aspartic acid concentration of 5.308 x 1st moles/kg. 
Figures 18 and 19 are plots of C, vs Ch and C, vs Cy 
calculated from the raw titration data. The slopes of these 
lines were calculated by a linear least squares procedure. 


Figure 20 is the results of those calculations shown as a 
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Figure 18. Cy as a function of the total concentration of 


aspartic acid (C,) for the Zinc-Aspartic Acid system. 


The total zinc concentration was 9.3/0 x Lome moles/kg. 
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Figure 19. Ch as’ a function’ of the total concentration oF 
zinc (C,,) for the Zinc-Aspartic Acid system. The total 


aspartic acid concentration was 5.308 x 10 ~ moles/kg. 
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Pacure: 120, (dc,,/dc, ) and (dc,,/dc,,) asa puaai Of phe OG 
the Zinc-Aspartic Acid system. The derivatives were 


evaluated at a 5:1 iiganda to metal ratio. 
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SuncEeion Of"pH.» The) curve: for (dc,/dc, ) is very similar to 
that of aspartic acid alone, shown in Figure 16, because 
the titrations were carried out with a large excess of 
ligand; however, the curves are not identical. The curve 
for (dc,,/dc,,) indicates that slight complexation occurs in 
the range of pH 4.0 to 5.0 since (dc,,/dc,,) ToenotmexactLy 
zero, and that the strongest complexation begins at pH 5.0 
where the derivatives become much different from zero. On 
the basis of the tabulated function (dc,,/dc, ) and (dc,/dc,,) » 
and PIh] , and PIM]. calculated assuming no complexation at 
PHos the free zinc and free aspartic acid concentrations 
were calculated using Equations 9 and 10. The [L] and [M] 
values were then used to aid in the selection of a model of 
complexation. 

Aspartic acid has three coordinating groups, two 
carboxyl groups and one amino group. Near pH 5.0 the 
uncomplexed ligand contains a single proton on the amino 
group. This proton may be displaced by a metal ion to give 
the complex ML. Two of the tridentate aspartic acid 
molecules may easily surround a metal ion, forming the 
complex ML. These are the only complexes reported for zinc 
and aspartic acid. Other possible complexes, reported with 
other metals, and with glutamic acid, a ligand very similar 
to aspartic acid, include HML where the amino group remains 
protonated and ML 4 where presumably the Jor Vocareent 
group is dislodged and the ligand is bidentate. Tiesf1 rst 


step in deciding a model for the zinc-aspartic acid system 
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was to attempt to calculate Bh ay ta for all four complexes, 
HML, ML, ML. and ML3, and to compare standard deviation of 
the fits to the total ligand and metal in solution. MThis 
process is summarized in Table 28. In the first attempt 
only the complexes ML and ML. and the protonation constants 
of the ligand were considered. The standard deviation 
using the ligand and metal mass balances from pH 4.0 to 8.0 
were 0.15% and 7.22% respectively. Next, the species HML 
was added for the same pH range; the fits were 0.103 and 
3.83% respectively, a considerable improvement. This same 
model was tried for the pH range 4.0 -9.0 but the constants 
for both HML and ML were negative. Finally, the species 
ML, was included, using the pH range from 4.0 to 9.0. and 
the standard deviation of the fits were 0.22% and 3.52%. 
This indicated that all four complexes might be present in 
solution. 

The delta quantities, ACs AC, and ACy described in 
Chapter II, were used to further study the model of 
complexation. Figure 21 shows the values of these quantities 
calculated on the basis of the [L] and [M] data discussed 
above and estimated values of the ligand protonation 
constants. Below pH 5.0 there appears to be only slight 
complexation, as predicted by the (dC,,/dC,,) plot in Figure 
20. Between pH 5.0 and 7.0 AC, and Ac, are about equal, 
indicating that the complex ML predominates in that region. 
The ratio of Ac. to AC, increases from that point on to 
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Table 28. Determination of the Formation Constants of 


Complexes of Zinc and Aspartic Acid with pH, 4.0 


Species 
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Table 28 (continued) 


“pIL] = 8-068, pIM] = 3.028 at pH)= 4.0 
Range of data: pH 4.0 —- 8.0 


bas in a. Species HML has been added to the model of 
complexation. 


“as in a and b, except that the pH range has been expanded 


to 4.0 - 9.0 and species ML. included in the model of 
complexation. The models from a and b were not suitable 


for this pH. range. 
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this pH range and a complex with a higher ligand to metal 
ratio may be forming. 

The AC, curve decreases from pH 4.0 to 5.0, remains 
constant from 5.0 to 6.0 and gradually increases beyond pH 
6.5. It was expected that AC, would be zero at all 
pH ranges above 6.0, if not negative because of the formation 
of hydroxy complexes. This illustrates one of the 
limitations of using the delta quantities with experimental 
data. AC, was calculated on the basis of the ligand 
protonation constants, total ligand in solution and the 


experimentally measured C AC,, may sehereforevsutser from 


H 
Many errors. First, the [L] and [M] data have been 
calculated with a PIL], and PIM], which may be in error, 
second, the ligand protonation constants may be in error 
and, finally, the Ch Values derived from a, titration, cuave 
May be less than ideal since in the particular titration 
chosen Ch was maintained constant while Cu was allowed to 
vary through dilution. As was discussed in Chapter II, 
during the experiment only one of the components needed to 
be held constant. Cy and Cu used toycatculatce AC, and AC, 
values were fixed at the values of the middle titration, but 
the value of Cy varies through the titration and must 
therefore be taken from actual experimental data. The Ci 


values are also limited by the need to use a hydrogen ion 


concentration and, at extremes of pH this could contribute 


a large error to the AC, values .lf.the pH is, converted 


directly to ee Bae However, the delta quantities were 
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introduced for the purpose of aiding in the selection of a 
model and approximate values are generally "sutticientercr 
this purpose. 

The rise in AC, ae therends of the pl range, showman 
Figure 21, is most likely a result of the errors just 
discussed. A rise in AC, atehacghw pl eis unl ikedivs inthe 
zinc-aspartic acid system because there are no known 
protonated complexes at high pH. If some error existed in 
jue protonation constant for HL the effect of that error 
would decrease at higher pH where HL is being titrated to L. 

Figure 22 shows the AChr AC, and ACy curves 
calculated on the basis of the [L] and [M] data discussed 
above, and estimated constants for the species HoL, HL, HML 
and ML. The values of AC. AC, and AC, have essentially 
been reduced to zero from pH 4.0 to pH 6.5. The complex 
forming above pH 6.5 is ML.- 

Figure 23 shows AC yr Ac, and AC,, calculated using 
constants for HoL, Hoe ML. - There is a Slight 
mise: in Ac, and ACy above pH 8.5, indicating a complex 
ML... As discussed above, AC, has sufficient error 
associated with it that it should not be used as the basis 
of a search for protonated complexes above pH 7.0. The 


value of AC,, can be made more reliable with better sets of 
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[L] and [M] data and better values for the ligand protonation 


constants. It should be recalled that one advantage of the 
FICS method is that it does not attempt to explicwiehy? tat 


the Cy data, but vrather fitting the ligand and metal mass 
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balances is considered more important at this Pome jinvthe 
calculations. 

The next stage in the FICS calculation procedure is the 
search for good estimates of PIL], and PIM]. and the 
calculation of formation constants. Since the major com- 
plexation begins at pH 5.0 this was chosen as PHos and the 
initial free ligand and metal at this point was found. 
Jnitially, it was assumed that no complex formed, and PIL], 
was calculated to be 6.890 and PIM] , as 3.030. With these 
values and the functions (dc,,/ac, ) and (dc, /dC,,), [Ey and 
[MJ from pH 5.0 to 8.0 were calculated. Assuming the 


existence of HML, ML and ML constants were calculated from 


2! 
the ligand mass balance,the metal mass balance and both mass 
balances together. The constants on iL were unreasonable, 
and a negative constant was found for > -M. The constants 
for both balances combined, )-ML, were acceptable (good 
eonstants tor )ligand protonation) >) fhe ‘concentrations of 
HML and ML were calculated and the initial free ligand and 
free metal concentrations were adjusted accordingly. A new 
set of [L] and [M] values were calculated, also a new set of 
formation constants. After three such iterations, the value 
for HML based on mM was no longer negative. The iteration 
was Says Wea out 6 times, after which the values of PIL], 

no longer changed and PIM], only “changed “by "07002 units: 

The constants changed only slightly with the last iteration. 
During this iteration the values of p[L], and PIM], changed 


less and less with each repeat of the procedure. efor 
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example PIM], Vat bedeaSins 1050. 134.0375 s 204 oe 32048,. 338052, 
3.055 and 3.057. This sequence indicates that PIM], was 
3.057 + 0.003 at pH 5.0. Unlike the sequence described to 
give the acid dissociation constants for aspartic acid, 
this sequence does not appear likely to produce a new 
estimate of PIL], and PIM]. exactly the same as the previous 
One. in this casé the iteration can only be carried on 
until the changes are as small as can be justified by the 
expected accuracy of the formation constants. 

Table 29 shows formation constants derived with the 
data above pH 5.0. Two alternatives are shown, the first 
in which HML is not included in the model of complexation 
and the second in which it iS considered. In the first 
case iteration stopped immediately because very little ML 
existed at pH 5.0. The second set of constants is the final 
mBesult of the series of 6 iterations described above. fiiis 
Of Interest to note that the constants zor ML.and MLo, 
calculated on the basis of these two models, are almost 
identical. There is insufficient evidence to prove that 
mMbaactually exists in Solution. “\First, the formation 
constant for HML is small and, second, its inclusion does 
not significantly change either the standard deviation of 
the fit rue values of formation constants of the other 
Species in solution. ~A species distribution for the model 
including ML and ML, is shown in Figure 24. 

The FICS method would be most convenient if it was 


possible to begin at a pH, at which no complexation had 
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Table 29. Determination of the Formation Constants of 


Complexes of Zinc and Aspartic Acid with PH 50 


Species ) L Nii >? ma 


HL S700 Shi ase) 
ML | D0 SAIS, 52166 
ML. TORS EW Us iON) 10.14 
Seo tee its 0.963% 1.76% SD,_=1.08% 


L_ c. 
SDy 4.26% 


HD he ONe, eles), 
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ML D.O2 Dish De Oc 
ML. Oe ey AS nO. 0'O LOLs 
SD POty ert: 0.96% 1.49% SD_=1.02% 


L_ 2 
SDyn4- 74% 


“p{t] (= 6.890, Di Mre 3.030, on the basis of no complexation 
at PH = 5.0. grange: OL datas pH 5.05 — 32.0. 


Spies: and PIM], calculated by the iterative procedure 
described in the text. P[L] \= 6283.1, BINS Siien! sue 
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occurred and where PIL], and PIM] | could be calculated 
easily. This cannot be done because the values of (dc,,/dc, ) 
and (dc,,/dc,,) become less reliable at extremes of the pH 
range, and these are usually the only regions where 
complexation may be absent. In addition to the problems 


caused by the uncertainty in C, at extremes of pH, a second 


H 
problem became evident during the zinc-aspartic acid 
experiments. 

At pH 4.0 very little complexation occurred. This 
could not be considered an extreme pH, yet the reliability 
of (dc,/dc,,) is guite low. This is because the titration 
curves were very Similar in that pH range and, consequently, 
the AC, term in ( 4C,/ 4C,,) was small. The difference in 
the titration curves caused by the normal errors was larger 
than the effect of significance, namely the effect of 
complexation. Integration over a set of (dc,,/ac,,) in which 
no complexation is occurring only accumulates the errors in 
those derivatives and provides almost no information about 
complexation. The most reliable data appears to be obtained 
if more than 10% of the metal (conditions of large excess 
of ligand) is complexed. 

Table 30 summarizes a number of the protonation 
constants of aspartic acid and formation constants of 
complexes between aspartic acid and zinc in the literature. 
Martell and Smith (60) consider that the best reported 
values for the formation constants of the species HL and HL 


are 13.33 and 9.63 respectively. The best values reported 
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taple 30." The Protonation Constants for Aspartic Acid, and 
the Stability Constants of Zinc-Aspartic Acid 


Complexes Reported in the Literature. 
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in this work are 13.356 and 9.678 respectively, found in 
Table 27. The ligand protonation constants determined 
Simultaneously with stability constants, shown in Tables 

28 and 29, also agree well with those in the literature. 
mae constants for ML and ML. reported by Chaberek and 
Martell in reference 61 were 5.84 and 10.15. The constants 
for those species, shown in Table 28, are in good agreement 
with these values. The constants for ML appearing in 


Table 29 are about 0-2 log units smaller than those in 


are similar to those 


Table 28, while the constants for ML, 


in Table 23. 

The stability constants for the model including ML and 
ML, shown in Table 29, were used to simulate a titration 
curve which had been obtained experimentally. The curves 
are shown in Figure 25. The agreement between these curves 
is very good, especially in the region used to calculate 


the formation constants of ML and ML.- 


D. Summary. 


Many of the calculation techniques described in Chapter 


II were applied to data collected to study the complexation 


Sreaspartic acid to zanc. It was observed that there were 


no sharp discontinuities in the curves of C, vs C,. The 


onset of precipitation may have been responsible for this. 


Receries of titrations were performed at a styl aed a Norchelel Gute: 


metal ratio, and each component was varied about 10%. It 
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Pagure.252 Experimental,o, and simulated, ay, titration 
eurnves Ot a solution containing 9.370 x Noe moles/kg of 
zinc and 5.308 x 10> moles/kg of aspartic acid. The 
titration was simulated using a model of complexation 


including the species ML and ML. (See Table 29). 
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was observed that below pH 5.0 only Slight complexation of 
aspartic acid to zinc was occurring and as a result the 
values of (dc,/dc,,) became Less reliable. “This dafficulty 
can in part be avoided by using a wider range of Ch and Cur 
especially if C. > Cur The series of titrations used to 
Study the complexation of glutathione, described in the next 
Chapter, were run at less of an excess of ligand and a wider 
variation in Ch and Cur both corrections designed to keep 


AC, Significant over as wide a range of pH as possible. 
The methods used to find the best estimate of PIL], 
and PIM], were found to work in the manner described in 
Chapter II. It was found however, that there were many 
Situations where this procedure would fail. For example, 
if too wide a pH range or a pH range containing too many 
complexes was chosen, the PIL], and PIM]. could change 
continuously beyond values that would be considered 
Peasonable. This will occur if, within each lteration, 
PIL], and PIM], change in the direction of more complexation. 
As a result the values of each set of new formation 
constants is larger, and the new PIL], and PIM], must be 
calculated assuming even more complexation. If the data is 
good and the model of complexation appropriate, this problem 
should in general be avoidable. In the discussion of the 
complexation of glutathione to zinc, in the following 
chapter, an alternative method of finding p[L], and PIM] , 
based on a trial and error search for a minimum in standard 


deviation is described. There is no reason to believe one 
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method is better than the other and, if the data is good, 
both should give reasonable, if not identical, values of 
P[L], and PIM] .- 

Finally, it was found that the protonation constants 
and complex formation constants derived with the various 
mass balances and different values of plL], and PIM], were 
rather insensitive to those changes, as may be seen in the 
constants for ML. Shown in the various parts of Tables 28 
Blicmco.) es tN1iS 2S a good indication of the’ relaability of 
that constant. The constants for species found in low 
concentrations and at the extremes of pH, and for the other 
species which coexist with those in some regions of pH, 

May be somewhat less reliable. The final test, as described 
Bopeheapter iti, is the Simulation of) a titration using, che 
constants derived from the FICS method and comparing this 
titration to actual experimental results. Good agreement 


between these is expected if the model for complexation and 


the stability constants are correct. 
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CHAPTER VI 


THE COMPLEXATION OF ZINC, CADMIUM AND LEAD BY GLUTATHIONE 


A. Introduction. 


Glutathione is a tripeptide composed of Glutamic: acid, 


cysteine and glycine residues. 


9 Q 0 Q 
HO-C-CH—CH~CH —CH—-cH-t-—10Ht—-CH 6-H 
+ 
NH, CH 
SH 


The amino and carboxyl groups of the glutamic acid residue, 
the sulfhydryl group of cysteine and the carboxyl group of 
the glycine residue have ionizable protons and are the 
groups most likely to participate in coordination to metal 
ZOnS.. At pH's above. 12.0 protons on, the peptide linkages 
may also be removed. 

The acid-base behavior of glutathione must be described 
in terms of a series of microscopic ionizations. At low 
PH all four groups are protonated and a single form of HL 
exists in solution. Near pH 2.0 a carboxyl proton from 
either the glutamyl ox glycine group may be removed, giving 
two possible coexisting species with the stoichiometry HL. 
Since the glutamyl carboxyl is slightly more acidic, the 


Species deprotonated at that site predominates; however, 
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both do exist at detectable concentrations (72). 
In the pH range from 3.0 to 4.0 the second proton is 


removed and a single species HoL exists in solution. 


Q Q Q O 
O—C—-CH—CH Fea gC -NH-CH-C—-NH-CH »—C-0 
+ 
NH, CH, 
SH 


Above pH 8.0 the process repeats itself, the amino and 
scufaeeeee groups losing protons simultaneously. The 
sulfhydryl is slightly more acidic and as a result the 
species having that group deprotonated is the predominant 
form of HL. The final proton is removed above pH 9.5. 

Glutathione may bind to metal ions through the four 
groups discussed above. Most transition metals and "soft" 
metals such as mercury bind predominantly through the 
sulfhyaryligroup, often) displacing: ‘the sulfhydryl) protonvar 
HOw pHs “The other groups of the molecule may be protonated, 
deprotonated, bound to the same metal ion as the sulfhydryl 
or bound to another metal ion. This could potentially lead 
to an extremely complicated mixture of species existing in 
solution. 

In this chapter a study of the interactions of 
glutathione with zinc, cadmium and Leadvusing thes FICS 
method asl described: othe! pH titwations werespexrtomied 
using the optical couple gravimetric apparatus described in 


Chapters III and IV, and the data was evaluated using the 


methods developed in Chapter II. The imc ae ch = 


ee 


vs ie 
SOD 


vi re 
i 
a 
‘ Si 4 
iw | a ats $y ee J 
ke dtd ey fy iad \ 
ioe. plea Cane Gee want Bina 
| a a ale wins ; Mi i Ta “at i” nee 
| 7 roe ) ne py 
' i i; t Wine ae - 
; * wg 1; j i Sis i P Pie ea ae om Es ry’ mm ag ee saab ‘i 


, : The, 7 he 7 ape Pe ans hi soe a ‘ ann , el oy 4 “s 
eee ota, ei viet 
Pit 4 Spit A el eB ur "a ae ee | Sia jana ‘teniacere ids’ 


(. OZ eyes caaek iw, ane er ae a 4 (inh 

yi: maitre sont ‘tha a j pire ¥ wii a aon a 

ea ety gyn car teal rere bog rei. cst eiteesaiall 

“tiga a a ‘nivel natn ‘a out: ‘-ettieaeet a “tat Fae 

wi thew’ ine ey inthis like pecans abate hating’ oct pnt a 

wey tae boonies a ey ta% 4G0; ae il bas Ee 


\. 
‘i 


Ieee lee sao Gan gribed ee " 


£60 
divided into two sections, the first dealing with the 
experimental results, the application of the FICS method 
and the deduction of species present in solution, and the 
second dealing in “tae detail with the chemistry involved, 
the structure of the complexes and a comparison of the 


results with those of other workers. 


B. Results. 


Determination of the Acid Dissociation Constants of 
Glutathione 


A series of solutions containing from 1.928 x ieee. ves 


es Oil xX Howe moles/kg of glutathione were titrated from pH 
Beco to 10.9. Details of the composition’ of these samples 
is shown in Table 31. The concentration of titratable protons, 


Cc was, calculated at-each O21) pH unit using the methodic: 


H’ 
enterpolation mescribed in Chapter II. At each point? the 


change in Cy with total glutathione concentration was 


calculated, giving (dc,,/dc, ) as a function of pH. The 


values of Cy and (dc,,/dc, ) are summarized at 0.5 intervals 


Onerables 32 and: 33. 


The concentration of deprotonated glutathione, [L], 


was calculated using Equation 9 and the values of (dC,,/dC,) 


are shown in Table 33. Since the initial free glutathione 


was not known at any pH, a series of estimates were tried 


at pH 10.5. High pH was chosen because in this region [L] 
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Table 31. Determination of the Acid Dissociation Constants 


of Glutathione:Solution Composition Data. 


Weight u Solven Concentration Concentration 
Glutathione Weight ofsAddedi Acid® cfsGluitathione 

Sample gm gm moles/kg moles/kg 

1 102002 790.002 91. Dlioe los: Henke) SGP 

2 8.000 1o2"005 eee 08 Sasso oe, 

3 4.000 POSC998| ieee ce Toes Ton 

4 11.001 190.002) 4 05M ROR MS s30lee On. 

5 6.003 104-001. Goes oes oboe Oe 


“The glutathione concentration was 0.09636 moles/kg. 
Pohe solvent contained 1.159 x 10. moles/kg of HNO. 
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Table 33. 


Determination of the Acid Dissociation Constants 


of Glutathione: The Derivatives (dc,,/dc 


L) Cyr PH 


and Calculated Free Ligand Concentration Using 


PIL], ODS eOAUr at PH, LOS 
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would be nearly equal to Che thus simplifying a search for 


the best PIL]. With each P[L], a set of [L] values was 
calculated, and with these sets of values, protonation 
constants of glutathione were calculated. The best values 
of PIL], and the constants were chosen on the basis of the 
Pet.O1 the [L] data and constants to’ the total’ glutathione 
in solution. For a solution having a concentration of 

0) Xx 0m moles/kg in glutathione, the best value of PIL], 
etepH 10.5 was found to be 3.040. The [L] data calculated 
feom ‘this PIL], are shown in Table 33 and the protonation 
constants for glutathione are shown in Table 34. These 


results are compared to those of other workers in Table 35. 


The Determination of the Formation of Complexes Between Zinc 


and Glutathione 


Two series of 5 titrations each were carried out on 
solutions containing varying concentrations of zinc and 
glutathione. The solution composition data for these 
titrations is shown in Table 36. The first set of five 
was designed to measure the change in Cy with changing 
concentrations of metal ion with constant ligand concen- 


tration. The concentrations of titratable proton for zinc 


= te) 
concentrations from 7.820 x 10 tO el 73 eee moles/kg 
are shown as a function of pH in Table 37. As required by 
the FICS method, concentrated glutathione solution was added 


during the titration to maintain constant ligand 
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Table 34. The Acid Dissociation Constants of Glutathione, 
Overall 
Species Formation Lo pK 
Constant 3 (8 ) a 
23 
H,L 283. 10 23.452 ee 52 
HL a O75ex 105-2100 3.460 
HL MERIC Se aii ees ae 8.640 
HL ua eau 9.501 Secon 
maple 35. -The Acid Dissociation Constants for *Giutathione 
Reported in the Literature. 
pK Conditions Ref 
rs al et dt ee La LOneas CLonS KE 
HL Tem af 
H,L HL HoL a 
2.05 Sao SERIE 9.49 ha 
8.74 9.62 252 CF 20m io 73 
3°59 SL} 9.65 25°C) 100.25 48, 
a 49 
8.96 9735 50 
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Table 36. Solution Composition Data for the Zinc-Glutathione 


Experiments. 
Concentration Concentration Concentration 
Sample of Glutathione, GEEZINC,. of Added Acid, 
moles/kg moles/kg moles/kg 
=e =f =A 
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7 Bees 106r 9.775. x 10° ae os 
8 ete aakuane G75 a HO ie hai vicy cay tna 
9 eo x HO 9.775 x LOT. L066 uO 
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eoncentration (3.855 x owe moles/kg). The stock glutathione 
solution was used for this purpose (see Table PY Wom 9 MAW ates 
quantity added during a typical titration was less than 0.5 
gm. 

The second set of titrations was performed at varying 
glutathione concentrations with the metal concentration 
constant. The Cr data for variation in glutathione from 


£07,4..981:9" x Me moles/kg at constant zinc of 


gee91 x 10° 
om775 x 10" moles/kg is shown in Table 38. 

The derivatives (dc,,/dc, ) and (dc, /dc,,) were evaluated 
ao. ivph Antervals from pH 5.5 to 9.0. During) this step 
it became clear that some of the titrations deviated 
appreciably from the others. If the Cry values from three 
or four titrations lay on straight lines while the remainder 
deviated at every pH value, the deviating titrations were 
removed. In the case of the zinc-glutathione experiments 
two of five of each set were rejected at this point. Since 
these were in the middle of the ranges varied, this did not 
Sacrifice*accuracy Of the (dc,,/dc, ) and (dc,,/dc,,) quantities. 
In the case of the cadmium and lead experiments described 
later no titrations were rejected. The final values of 
(dc,,/dc, ) and (dc,,/dC,,) are shown in Table 39. As this 
data indicates, only small amounts of metal are complexed 


at pH 5.5. It was therefore not necessary to begin the 


calculations at lower pH values. Only when (dc,,/dc,,) 


reaches values appreciably different from zero can they 


provide a reasonable quantity of information about the 
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Table 39. 


The Derivatives (dc,,/dc 


) C,,/PH and (dC,,/dCy) « 


for the Zinc-Glutathione Experiments. 


(dc, /dc,) 


Hoe ie 18) 
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complexation occurring in solution. If the calculations 
begin at a lower pH, the experimental and calculation errors 
introduce uncertainty into [M], yet provide no information 
about complexation. It is therefore advantageous to select 
PH ranges where the most information about complexation is 
available in order to offset, at least in part;, the 
detrimental effects of error. 

The final calculation of the FICS method provides the 
free ligand and free metal concentrations from (dc,,/dc, ) 
and (dc,,/dc,,) - However, P[L], and PIM], must be known. The 
search for PlL], and PIM], may involve a trial and error 
process where a series of estimates of PIL], and pIM], 
values are chosen to try to minimize the standard deviation 
of the fit of the model of complexes to the total metal and 
ligand in solution. This may in theory pose some problem, 
Since initially the model is unknown. However in practice 
it was found that the choice of PIL], and PIM]. has no 
effect on the search for a modél of equilibria in solution. 
No case was observed where a model was simultaneously 
accepted at one set of p[L], and PIM], anc rejected wilthean 
alternative set. The model was searched for and found 
using PIL], and PIM]. assuming some reasonable amount of 
complexation at pH, and, with this model p[L], and pIM], 
were refined. 

The process of refinement Of p[L], and PIM], for the 
studies of the complexation of glutathione involved 


independent searches for the best values of p[L], and pIM]).- 
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It was found that the standard deviation of the fit to the 
metal mass balance, SDyr was constant if PIM] | was constant 
and, Similarly, SD, was constant if PIL] | was constant. The 
minimum in SDy and SD, were therefore searched out indepen- 
dently. A model for the equilibria in solution was accepted 
Only if a clear minimum in both the SDy and SD, could be 
found. 

The formation constants of the complexes of zinc and 
glutathione are shown in Table 41. These constants were 
calculated with the [L] and [M] data shown in Table 40. The 
values of p[L] jana PIM] .were chosen to give the best fit to 
this model. The relative SDy was 3.60% and SD. was 0.943%. 
The constants reported were calculated using both the metal 
and the ligand mass balances, while each was used indepen- 
dently while searching for a model. This has been described 
Pminchapter Ii. \ The constant for AL was given si) SPER EY/ op% 10° 
in order that calculations remain within the dynamic range 
of the computer. The constant for HoL was calculated 
simultaneously with those for the metal complexes and was 
Bound to be 1.36 x Tose. This agrees very well with 
IN sso to map: The calculated from the series of titrations 
involving glutathione alone ( the previous section, see 
Table 34). 


The deduction of the species present in the zinc- 


glutathione system proceeded in much the same way as 


described in Chapter II. This involved the calculation of 
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Table 40. 


pH 


The Free Ligand and Free Metal Concentration for 
the Zinc-Glutathione Experiments with 
PIL], = 9.596 and PIM], 


Glutathione 
[L], moles/kg 
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Table 41. The Formation Constants of Complexes of 


Zinc and Glutathione? 


Formation 
Species Constant ,/3 Log( /3) 
M(HL) pulse ie 14.50 
M(HL) . 1.48 x 107? 7s jal 
M(HL) (L) WED doce Oe 20a0e 
M(L) 5 ieee 10" Tens 
M(HL) , 3.25 x 10°! 42.51 
M(HL) . (L) gid Se 10 33.91 


“the standard deviation of the fit to the metal and ligand 


mass balances are 3.60% and 0.94% respectively. The 
18 
Protonation constant, for the Hob WAS le 3 Oe LO seem ne 


: 9 
mrotonation constant for HL was given’ as) 3.17 x L0=. 
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protons, ligands and metal whose chemical situation has not 
been found. The magnitude of these quantites and the 
relationship among them gives reasonable indication of the 
concentrations and stoichiometry of species which may be in 
Solution. Figure 26 shows ACLs AC, and ACy for the zinc- 
glutathione system, using [L] and [M] from the FICS 
calculations and the ligand protonation constants evaluated 
in an earlier experiment. Below pH 6, the complexes appear 
to contain about equal numbers of protons, ligands and metal, 
suggesting the existence of M(HL) species. In the range of 
PH 6 to pH 7, the situation changes to a 2:1 ligand to metal 
ratio with the number of protons lying between 1 and 2. 

This indicates that a mixture of M(HL) , and M(HL) (L) may 
exist in that pH range. Finally; above pH 8 the number of 
ligands per metal is: greater than wwand the complexesvere 
being rapidly deprotonated as indicated by a decrease in 
AC,,- 

Te wwas found that M(HL), M(HL) 5, M(HL) (L) and M(HL) , 
were major components in solution in the pH range studied. 
At the upper extreme of pH a series of possible species 
‘including M(L). and M(HL) . (L) were suggested but,since the 
Hengesof data is~-narrow;,-proof-of their existence is not 
conclusive. 


The middle titration used in the FICS method was 


simulated with the constants shown in Table 41. The actual 


experimental curve and the simulated curve agree very well 


Homies .5,° This is considered the final and most reliable 
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indication of successful application of the FICS method £0 
titration data. The constants were extracted from the 
ligand and metal mass balances and, for reasons described 
moeohnapter If; not from the experimental Cu values. Asa 
result, the model for complexation and the values of 
formation constants are not tailored to fit the Ci Calta 

In comparison, other computer-based calculations,such as 
SCOGS, LETAGROP and so on, are specifically designed to fit 
Ch data. If, therefore, a titration is simulated using the 
model and constants from SCOGS, a good agreement with the 


experimental values of C.. is assured while, on the other 


H 
hand, no such agreement iS guaranteed using the FICS method. 


Good agreement will only occur if the model of complexation 


moncorrect and mi the values of the formation constants are 


accurate. 


Determination of the Formation Constants of Complexes 


Between Cadmium and Glutathione 


Ten titrations were performed, five on solutions in 


ate 
which the concentration of glutathione was 3.907 x 10 


moles/kg while the concentration of cadmium was varied from 


5.820 x 1074 to 1.358 x 10° moles/kg, and five on solutions 


in which the concentration of cadmium was held constant at 
08. x hag © moles/kg while the concentration of glutathione 


oe ag 
was varied from 2.604 x 10° to 5.206 x 10 moles/kg. 


Details of the solution compositions are shown in Table 42. 
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Solution Composition Data for the Cadmium- 


Glutathione Experiments. 


Concentration 
of Glutathione, 
moles/kg 
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The concentrations of titratable protons for the series of 
titrations in which the cadmium concentration was varied 
and the concentrations of titratable protons for the series 
of titrations in which the glutathione concentration was 
varied are shown in Tables 43 and 44. Plots of Cy vs Ch 
-and Cur revealed that no titrations deviated systematically 
from the expected straight lines and therefore all were 
used in the calculation of (dc,,/dc, ) and (dc, /dc,,) . The 
values for these derivatives are shown in Table 45. 

The procedure used to determine the complexes existing 
in solution and the method to obtain the best values of 
PIL]. and PIM]. have already been described. The concen- 
tations of ligand and metal are shown in Table 46, and the 
constants for the formation of complexes between cadmium 
and glutathione are shown in Table 47. 

The preliminary steps in the deduction of a model of 
cadmium-glutathione interactions required the calculation 
of ACs AC. and AC i these are shown on Figure 27.. 
Because of inaccuracy in the value of PIL], at this stage 
in the refinement of the data, the curve for AC, appears 
to be too low by agenesis There is also some 
uncertainty in Ac, and Ac, due to error in p[L], and pIM], 


but these errors will have no effect in later stages of 


calculations since AC. AC, and AC, were introduced to 


aid in the selection of a model. If the error becomes too 


large however, the delta quantities cannot be of much help 


in the selection of a model. 
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and (dc,,/dc 


for the Cadmium-Glutathione Experiments. 
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Table 46. The Free Ligand and Free Metal Concentrations for 


the Cadmium-Glutathione Experiments with 
PIL], = 12.679 and PIM]. =.3,080rat PH, 4.0 


Glutathione Cadmium 
pH [L], moles/kg [M], moles/kg 
4.0 eeOoes 0s paso x ons 
3 e038 x 105-2 e950 1077 
5.0 Delqi xt 100 8 2.78 x 107° 
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Table 47. The Formation Constants of Complexes of 


Cadmium and Glutathione? 


Formation 
Species Constant 3} Log( 3) 


iS 


M(HL) Aes h0 Nei 
M(HL) . Ta Gearel OR 4 31.06 
M(HL)M(HL) , 1.60 x 107" 50.20 
M(HL) (L) 5.02 x 10 21.70 
M(HL) , 2.57 x 10" 44.41 
M(HL) 5 (L) 8.86 x 10 A505 


“The standard deviation of the fit of the metal and ligand 
mass balances are 1.79% and 0.46% respectively. The 

protonation constants for HL and HL were determined to be 
Dez OFX foce ATIC uel so ex nee respectively; the protonation 


constant for HL was given as 3.17 x 10°. 
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ie Che range “Ol * pn 4.0 to-S. 0 the predominant complex 
appears to be M(HL). In the range of pH 6.0 to 7.0 there 
are two ligands and two protons per metal, suggesting M(HL).. 
From pH 7.0 to 8.0 more than two ligands and two protons 
are complexed to each metal indicating that a mixture of 
M(HL) , and M(HL) , has formed. In the range of pH 8.0 to 
9.0 the number of ligands per metal continues to increase, 
but the Ac, value drops rapidly above pH 8.3. The ligand 
appears to remain bound to cadmium exclusively as (HL) until 
pH 7.5. The major components of the cadmium-glutathione 
system were found to be M(HL), M(HL) , and M(HL) ,. 

With the constants shown in Table 47 the middle 
titration was simulated and compared to the experimental Cy 
data. A small deviation was observed in the pH range of 


6.5 to 8.0. This will be considered in the discussion 


section of this chapter. 


Determination of the Formation Constants of Complexes 


Between Lead and Glutathione 


Ten titrations were performed, five with varying metal 


and five with varying ligand concentrations. The details 
of the solution compositions are shown in Table 48. The C,, 


data for those titrations in which the metal concentration 
i hose 
was varied appears in Table 49 and the C,, fara 


titrations in which the ligand concentration was varied 


ad CC. indicated 
appears in Table 50. Plots of Cy vs Cu an L 
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Solution Composition Data for the Lead- 


Glutathione Experiments. 


Concentration 


of Glutathione, 
moles/kg 
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that no titrations gave points which deviated systematically 
from linearity and therefore none were rejected. The 
derivatives (dc.,/dc,) and (dc,,/dc,,) are shown in Table 51. 
The free ligand and free metal concentrations, shown in 
Table 52, were calculated from the (dc,,/dc, ) and eae 
data and PIL], and PIM]. chosen by methods described 
previously. The constants of complexes forming between 

lead and glutathione are shown in Table 53. 

The quantities AC ys AC, and ACyr shown on Figure 28, 
were used in the deduction of a model for the lead- 
glutathione interaction. As was observed in the cadmium 
system, AC, is displaced to slightly smaller values by 
eapout: 1 x om moles/kg. The curves are almost identical 
to those of the cadmium system. The only observable 
difference is that the concentration of the complex M(HL) 
is expected to be slightly higher since the values of AC, 
AC, and ACy are larger in magnitude within the pH range 
4 to 5. The deprotonation of the (HL) unit occurs in 
exactly the same pH region as for the cadmium system 
indicating that loss of that proton is not occurring by 
displacement by the metal ion. 

The middle titration was simulated using the constants 
shown in Table 53 and agrees very well with the experimental 


data. This indicates that the model for metal-ligand 


interaction is complete and reliable. 
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Table 5l. The Derivatives (dc../dc 
y/ Cupane Base aU elas 


for the Lead-Glutathione Experiments. 


pH (dc,,/dc ) (dc,,/dc ) 


Far L Cyr PH M C,7PH 
3.5 2.488 | S0n0o8 
4.0 2.194 -0.199 
a5 1.992 -0.557 
5.0 1.922 ~0.902 
5.5 1.916 -1.118 
6.0 1.910 ~1.302 
6.5 1.888 -1.458 
70 1.842 -1.593 
UG eas fle bal -1.693 
8.0 1.634 -1.713 
8.5 1.432 ~1.530 
9.0 1.087 -1.073 
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The Free Ligand and Free Metal Concentrations for 
the Lead-Glutathione Experiments with 
33140) at PH 4.0 
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Table 53. The Formation Constants of Complexes of 
Lead and Glutathione? 


Formation 
Species Constant, /} Log( f ) 
M(HL) 74 x10 5° 16.24 
M(HL) cases aldo 30.93 
M(HL) M(HL) , 1.67 x 107 50.22 
M(HL) (L) 3.49 x 10 1 es4 
M(HL) , 736 x 10". 44.13 
M(HL) . (L) ciatiae cles 35.50 


@mhe standard deviation of the fit to the metal and ligand 


mass balances are 2.60% and 0.71% respectively. The 
protonation constants for HL and HoL were determined to be 
4.47 x oes and: tes ors Roos respectively; the protonation 


Gonstant for HL was givenvas 3.17 x 107. 
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Gc. Discussion. 


The complexation of zinc, asanten and lead by 
glutathione has been studied by several groups. Martin and 
Edsall (8) reported formation constants of complexes of 
zinc, Perrin and Watt (6) discussed zinc and cadmium, Fuhr 
and Rabenstein (7) and Williams (5) considered zinc, cadmium 
and lead. Fuhr and Rabenstein employed NMR techniques to 
study the proton and bee chemical shifts for nuclei an 
various parts of the glutathione molecule while the other 
groups used pH titration methods. Perrin and Watt, and 
Williams used the program SCOGS to evaluate the pH titration 
data. 

The results presented by the groups using SCOGS 
disagree with the NMR results. More specifically, both 
groups concluded that species existed that were impossible 
according to the chemical shifts observed in the NMR spectra 
ef Glutathione. It should be noted here that,thegpapers 
were published in the order Perrin and Watt, Fuhr and 
Rabenstein, and Williams. The last group had benefit of 
the guidance offered by the NMR results in the selection 
of models of complexation, but chose models virtually 
identical to those of Perrin and Watt. 


The results for the complexation of zinc, cadmium and 


iead presented in this thesis agree with the chemical shift 


data presented by Fuhr and Rabenstein. In addition, where 


POssrbiec,. titrations were simulated on the basis of models 
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presented by Perrin and Watt and by Williams. These have 
been compared to data collected for use with the FICS method 
and to curves simulated on the basis of the models presented 


in this thesis. These titrations agree well with each other. 
Zinc-Glutathione 


The model of the interactions betwee zinc and 
glutathione proposed by Perrin and Watt and by Williams 
included the major complexes M(HL), M(L), M(HL) 5 and 
M(HL)(L). The NMR data of Fuhr and Rabenstein suggested 
that no complex M(L) formed in the pH range suggested by 
those workers. The complex M(L) may, or may not, be real. 
Molecular models show that it is possible because the 
sulfhydryl group and the glutamyl amino and carboxyl groups 
may bind simultaneously to a metal, leaving the ligand in 
a completely deprotonated form. Two arrangements are 
possible, one in which all the coordinating groups are 
planar and a second in which they are immediately adjacent 
to each other. With this second form, two ligands may bind 
Simultaneously as M(L) 5- Fuhr and Rabenstein offered two 
criticisms of this model. They argued that there is no 
reason to believe that the metal has sufficient affinity 
amino proton Of the glutamy! 


for the ligand to displace the 


residue and, that even if such affinity existed, a large 


and unstable ten-membered ring would form. Molecular 


models show, however, that the carboxyl and amino groups of 
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the glutamyl peptide linkage are very close ‘to the, metal 
ion, though not directed toward it and May contribute to 
the strength of binding of the terminal glutamyl residues 
and stabilize the large ring. However, the NMR studies of 
the interaction of zinc and glutathione Clearly show that 
no binding to the glutamyl terminal occurs before DELO. 
Both Perrin and Watt, and Williams suggest about 10% of the 
metal bound in a complex such as M(L) at pH 6. At pH 7 
Perrin and Watt suggest that 55% of the metal exists as 
M(L), the NMR results show about 25%. A species 
distribution for the constants derived from the FICS method 
(Table 41) is shown on Figure 24. No species M(L) was 
found to exist. Deprotonation of the amino group of the 
glutamyl residue and complexation to the metal begins at 

PH 6 when M(HL) , loses a proton to form M(HL) (L). There 
may be in fact two forms of complexes with that stoichio- 
metry, namely a binding and nonbinding form. However, 


judging by the pH range of deprotonation it appears to be 


a displacement type reaction. At pH 7 about 30% of the 


metal is complexed in a form M(HL) (L). 

Martin and Edsall have suggested the three complexes 
M(HL), M(HL) , and M(HL) (L), reporting a log(K)) for the 
= M(HL) of 5.10 and a pK for thegiiese 


reaction M + HL = 


£OnLzZation oL M(HL) , Of" 7750. Perrin “and Watt, found) 4. 74 


and 7.04 for these reactions; Williams obtained the values 


4.88 and 7.35, From the data in Table 41 values of 5.00 


and 7.15 can be calculated for these equilibria. This is 


eo 
a A (a 


iy UO die a 7 
uN ans ee He | fe iene ie. vat 
i i Wis i ki, ‘a nae, A imetook A Hi ai 
alah hia | sh if dha Ht inh ‘ich eae pd awn 
cantina: ae ep, hie sii hue wake a, hea 
sgt! heh ea ia! pada th, mil nm 
| a tee DA svn Me hey vik ae pammiek Ce 
Be iw + 0a ies nce nae ve | Wk aa 
ned eb : ol bt, net, ua gu ab patier theobtt fem “abd oe 
ight | he tay bao pine yay ange ie eats 
Lie aa ic aa ‘ua na ik ‘he ata he au * me 
and Se vy wea thi Ween: nt ptm Ye | 
situ Yank nets i th si Velgta pina hie _ nae 


ey 


Ai | 
A 


198 


© (TH) W 


*I 


"tA “(1)°(tH)w ta *e¢qyw 
‘ (TH)W 


ATS (T) (TN) We -f eee 
‘SUOTUFERNTDN-OuTZ AOT UOTINGATAASTp soetosedg 


OS 07 


© tH) W 


S6z. eanDi yt 


$0 


‘?) 
‘@) 
3 
° 
re 
ol, 
© 
a 
-~, 
pa 
1@) 
(0) 
oO 
“y 
x 
«© 
wT 


© 


N 
— 


ee 


ay) 
‘ 
Dp 
; 
i 
aT 
. 
alae y { 
Pe ‘ . 
Ve 
aR 
Fig 
‘ye 
wo ie 
a 
St 
A 
a0 
fi 
ey 
A 
i 
. ee 
* 
Vid 
Fes 
> = 
tik ry A ; 
¢ Tr 
, a4 
; 
= | } e 
Si jy, 
q y 
‘i 
wal r; 
; af 7 
0, 
7 
‘ Fe 
yl? 


os i a 
bn es 

4 AAG 
ve at 
Me MY 


| m - i | 


. 


19,9 


a reasonable situation since Perrin and Watt, ‘and Willians, 
by including M(L), have probably underestimated the constant 
for M(HL) and, Martin and Edsall, without adequate knowledge 
of the complexation above pH 7.0, have overestimated the 

PK, (oh M(HL) ,. 

Previous studies have not take into account that 
complexes with a ligand to metal ratio greater than two 
Might exist, even though experiments at glutathione to 
metal ratios up to 8:1 have been used. The complex M(HL) , 
was found to be present using the FICS method, accounting 
for 20% of the zinc complexed to glutathione at pH 7.8. 
Molecular models show that three ligands may bind to the 
metal through the sulfhyrdryl group with no steric inter- 
Bactelonsceatmtherother parts of ithe) molecules; The netycharge 
On such’ a complex is <4. However, since the entire complex 
PS quite Yarge and atuleast 3 of thercharges sare located 
on the glycine carboxyl groups at the remote extremities 
of the eee this is not an unreasonably large net charge. 
No evidence was found for a species M(HL) ,, which would have 
eaenoucharge: 

The models for the interactions of zinc with 
glutathione suggested by those working with SCOGS and this 


work with FICS appear quite different and it was of interest 


to compare titration curves simulated using the alternative 


models. The curves were found to agree quite well: This 


is reasonable since it is conceivable that different models 


may give rise to the same titration curves. This problem 
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has in part been overcome by the use of the [L] and [M] 
provided by the FICS method, which considerably restricts 
possible models to explain the data. ACCCNDES, tOmrat cone 
models of Perrin and Watt, and Williams with the [L] and [M] 


[M] data from the zinc-glutathione experiments failed. 
Cadmium-Glutathione 


The interactions between cadmium and glutathione 
appear to be quite different from those involving zinc. 
Complexation begins at much lower pH and appears to strongly 
favour complexes bound only through the sulfur of the 
cysteinyl residue. Perrin and Watt suggest that the major 
complexes are M(HL), M(HL) 4, M(HL) (L), MoL and ML. They 
Found that®tat pH’5.8 about)» 20¢e0f the metal is botind as MoL. 
Williams concluded that the major species are M(HL), M(HL) 5, 
M(HL)(L) and M(L). Their results differed from those of 
Perrin and Watt in that they found about 50% of the cadmium 
held as M(HL) at pH 4.0 compared to less than 5% for the 
other workers and they found no MOL. The discrepancy here 
indicates that the model of complexation may be rather 
difficult to establish with potentiometric data alone. The 
NMR results of Fuhr and Rabenstein indicated that both of 
the above studies were in some aspects incorrect. A model 
for the cadmium-glutathione system must explain: (bine 
binding to suifhydryl at pH >9220, (2) binding *toegiurcamy! 
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in the pH range 2-10. The results of the FICS method are 
able in part to explain these observations; however, the 
model remains incomplete. 


The major species found included M(HL), M(HL) M(HL) , 


2" 
and M(HL) (1). A species distribution is shown on Figure 
30. All of these complexes are probably bound exclusively 
through the sulfhydryl group and the deprotonation of 
M(HL) , occurs without inten with the metal ion (the 
PK, for this ionization was 8.46). In addition-to these 
complexes, a polynuclear species M(HL) M(HL) , was found, 
accounting for about 40% of the metal’ at pH 5.6. It is 
most interesting to note the similarity of this complex 


With MoL found by Perrin and Watt; first, the maximum in 


2 


concentrations occurred at the same pH; and second, both 
accounted for 40% of the complexed metal at that pH. The 


formation of MoL from HL would release 2 protons, while 


the formation of M(HL)M(HL) . would release 3. Perrin and 

Watt's model compensated for this discrepancy in the number 

of protons titrated to pH 5.7 by including M(L) at that pH. 
The polynuclear complex M(HL) M(HL) , is able to explain 


the observed binding of the metal ions to the glycine 


carboxyl group. The glycine group would be able to form 


bridges to the metal ions of less than fully coordinated 


complexes thus forming a bonding between M(HL) 5 and the 


singly coordinated M(HL). This process does not itself 


release protons but it is observed by the effect it will 


have upon other equilibria involving proton displacements. 
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The model suggested by the FICS method does not account 
for the observed binding of cadmium to the glutamyl amino 
group. The species distribution shown on Figure 30 
indicates that a deprotonation of M(HL) , OCCUrs atwpHel, 
but this is probably not a displacement by the metal en 
The inadequacy of the model is further shown by a difference 
between the actual titration curves and simulated data in 
Eee=spH range 6.5 to’ 8.5.) This» strongly suggests ‘thatua 
displacement reaction does occur in that pH range, with 
binding to the amino group of the glutamyl residue. The 
models used by the other workers do not offer much help 
Since they propose formation of M(L) as low as pH 5.0 and 
M(HL) (L) at pH 6.0. With the data from the FICS method a 
model including M(HL) (L) indicates that it would only begin 
to form above pH 7.0 and would not offer explanation of the 
NMR results. The reason for these problems may stem from 
a polymerization reaction where the uncomplexed glutamyl 
and glycyl terminals may bind to metals other than that to 
which the sulfhydryl group is bound. This? type, ofasituation 
has been reported for glutathione complexes of nickel (54, 
55) and silver (74). In the latter report, silver- 
glutathione polymers were separated and characterized, and 
found to have between 15 and 20 metal-glutathione units. 

If series of short polymers are present in cadmium- 
glutathione solutions this would pose considerable 


difficulty for pH titration methods; first, a variety of 
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complexes might exist in a narrow PH range; and second, the 
mixture may not be at equilibrium. In either case, 
potentiometric titrations will fail to resolve the situation. 

In spite of these difficulties a great deal of 
information about the complexation of cadmium to glutathione 
has been extracted from the titration data. Biers te ecaste 
should be remembered that the deduction of a model of 
complexes existing in solution and formation constants is 
Only a part of the FICS method. The values for [L] and [M] 
Shown in Table 46, the values of the derivatives in Table 
45 and the delta quantities of Figure 27 have been derived 
from the experimental data and are independent of a model 
for complexation. They will not be affected by the choice 
of an alternative model. 

Below, pH,.7 .0..the,majority of metal is ~held gin athe 
complexes M(HL), M(HL)M(HL) , and M(HL) 5- These species 
account for both the observed pH titration curves and the 
NMR data. Models proposed by Perrin and Watt, and by 
Williams did not agree with the NMR observations Jactnis 


model also describes the complexes which may form in 


physiological pH ranges. In this respect it is interesting 


to note that glutathione in physiological media will 
invariably be bound to heavy metals only at the sulfhydryl 
group. The other parts of the molecule will continue to 


behave as if the metal was not present, perhaps binding 


to other metal ions, to enzymes and to proteins. 
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Lead-Glutathione 


The reaction of lead with glutathione has been found 
to be almost identical to that of cadmium. The curves of 
AC Ac, and AC, shown in Figure 23, and the species 
distribution shown in Figure 3l,are very similar to those 
Brox the cadmium system. The NMR data collected by Fuhr 
and Rabenstein showed that the interaction of lead with 
glutathione occurs: through the sulfhydryl and glycine 
carboxyl groups in essentially the same way as inthe 
cadmium system. However, there was a major difference; no 
complexation of lead to the glutamyl terminal was observed 
below pH 12.0. This suggests that the polymerization 
reactions of the type proposed above for cadmium are not 
likely in the lead system. The actual titration data was 
matched very well by a simulated titration curve based on 
the model including the species M(HL), M(HL)M(HL) 4, M(HL) 5, 
M(HL) , and M(HL). (L)- 

The model for the complexation of lead by glutathione 
described by Williams included the complexes ML, M(HL) (L) 
and ML.. All three of these species form by loss of an 


2 
amino proton from the glutamyl group below pH 8.0. Since 


the amino proton would usually be titrated above pH 9.0, 


these reactions must occur through displacement of the 


protons by metal ions. This does not agree with NMR data, 


which shows that the amino group does not participate in 


binding below pH 12.0. Williams made no attempt to use 
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the NMR data as a guide to selection of a model of 
complexation. 

As was described in the section on the binding of 
cadmium, some possible reactions of glutathione with heavy 
metals in physiological media are Suggested by the results 
of the experiments in this thesis. It appears that only 
very weak complexation occurs to groups other than the 
sulfhydryl within physiological pH ranges and, if binding 
doeS occur, it is most likely to another metal ion through 
the glycine terminal. This indicates that the parts of 
glutathione molecule that are not directly involved with 
binding behave as if the binding had not occurred. This 
has implications in the understanding of the effect of 
heavy metals on biochemical processes. For example, should 
a lead-glutathione complex bind to an enzyme through the 
glutamyl terminal, the activity of that enzyme may be 
effectively quenched. In this way metal ions which would 
not normally interact with enzymes near their active sites 
may,through a carrier molecule, have an effect on the 
activity of that enzyme. This also suggests that 
glutathione cannot be used as a heavy metal sequestering 
agent. Aside from the problems involving the other 
chemical properties of glutathione it cannot be used as 
a drug for treatment of metal poisoning because, as 
described above, it retains a great deal of its previous 


character even after complexation. As a result 20 Wl 


continue to participate in biochemical processes, perhaps 
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spreading the metal faster. This reduces the number of 
routes by which the metal ions might be eliminated from 


the body. 
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APPENDIX A 


THE COMPUTER PROGRAM USED TO SIMULATE TITRATION: CURVES. 


> DIM F$(20) 268020) yHO10) ob 610) 9M C10) 9 BC 10) oM1 (10) 9 C20) » 
L1¢10) sF1010) sF2¢10) 
10 FRINT “NUMBER OF TITRATIONS? "NINFUT SS 
POnPOR. Il=1. 10. S35 
30 PRINT "CONC LIGAND METAL SNINFUT S4(T) S701) 
AQ FRINT "FILE FOR STORAGE?" NINFUT FCT) 
50 NEXT I 
80 FRINT "INITIAL VOL? SNINFUT V1 
90 FRINT "NUMBER OF ITERATIONS? ®\INFUT Q1 
100 PRINT "INIT FHy NUMBER OF FPHePH INCR."\INFUT FeyNeTd 
410 FRINT "NUMBER OF FROTONS:BETA FOR MOST FROTONATED FORM? * 
NINFUT HsB 
120 FRINT "NUMBER OF LL INCREMENTS» RATIOsRATL??P"\INEUT NBe Rie RS 
130 PRINT “NUMBER OF M INCRs RATIOS RATM?P PT \INEFUT M39 RG oR? 
140 PRINT "NUMBER OF SFECIES?"\INFUT Ni 
150 FOR I=1 TO Ni 
160 PRINT "NUMBER OF HyLeMeAND BETAP®\NINFUT HOD) oh CI) eMC I) es BCT) 
170 NEXT I 
moO FRINY “ALE ACTO?T*\INPUT Ca 
220 PRINT "CONC OF NAOH? SNINFUT &S 
250 FOR S4=1, TO SS 
231 LET H$=F$(S4) 
23a OFEN Ht FOR OUTFUT AS FILE #5 
24) LET LC=S6(S49\CET M=S7(S4)\LET PHF2 
290 LET X1=F-LINLET X1=EXPC-X1#*L0G(10)) 
260 LET X2=EXPC(-FPXLOG(10)) 
270 LET C=C1+H#eL 
280 PRINT "CHO EQUALS" »C 
LET Fis=(EXP CHXLOG CXL) -EXF CH¥LOGCX2)))/EXPCHKLOG(K2)) 
SOc we LAr 20 
310 LET Li=L/CBXEXF CHXLOGCX1) >) 
315 LET M4=M 
$25 LET W=1.00000E-14 
330 LET L2=L\LET M2=M 
350 LET V=0 
360 FOR T=1 TO N 
370 LET Fis=EXP(-FXLOGC19) )D 
380 LET L3=LI\LET MS=M4 
390 LET Li=LixCitFi)\LET R=Lixk¢(Fit.8)/R1 
AOQO IF Fiz=0 THEN 410 \LET R=L1ix(F1~.8)/R1 
WG CEP IeR NGA ET IS=I2LEm Ciel t-(RZ720 
420 FOR Q=1 TO NS 
430 LET Li(Q)=Li\LET Li=Lit+I2\NeExT @ . 
A40 LET M4=M4XCLtF2)\LET R=M4k(F2+.8)/7KR6 
ASO IF F2z=0 THEN 460 NLET R=eMak(F2-.8)/R6 
460 LET I2=R/M3\LET I4=I2 
465 LET M4=M4~-(R/2) 
470 FOR Q=1 TO M3 ow 
480 LET M1L(Q)=M4\LET M4=N4+4+12\NEXT Q 
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600 
610 
oad 
612 
613 
6164 
617 
630 
640 
650 
660 
670 
680 
690 
700 
730 
760 
770 
771 
780 
790 
B00 
801 
B10 
820 
825 
B24 
830 
840 
850 
860 
870 
880 
900 
910 
930 
240 
920 
9460 
970 
1000 
1010 
1020 
1040 
1050 
1060 
1065 
1070 
1080 
1090 
1095 
1100 
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FOR Q=1 TO Nt 
LET P(Q)=H(Q)*LOGCFA)\NEXT Q 
FOR Q=1 TO Q1 
LET Y=1 
IF Q=1G0 TO 750 
LET I3=I3*1.,0001\LET 14=14*1.001 
LET Y=1 
LET [2=T3xRS/ (EXP ((Q-1) LOG (N3) +R) 
LET Li=FS5-.5xN3xI2 
FOR J=1 TO N3 
Detail (eed SET Sb 1+ TOWWNEX TO 
LET I12=14KR7/ (EXP (Q-1)*LOG (M3) 4R7) 
LET M4=F6—.5*M3K12 
FOR J=41 TO M3 
LET Mi(J)=M4\LET M4=M4402\NEXT J 
FOR K=1 TO N3 
FOR Ki=1 TO Ni 
LET PACKI)=L CNL) XLOGCLICK)) 
NEXT Ki 
FOR O=1 TO M3 
FOR O1=1 TO Ni 
LET F2(01)=MCOL)XLOG CML (OD) 
NEXT O41 
FOR Ki=1 TO Ni 
LET S8=LOGCRCKL) FFE (RL) FPL (KL) FF 20K 1) 
IF S8z=-35 THEN 826 \LET X=0 
IF $8<-35 THEN 830 \LET X=EXP(S8) 
LET Si=Si+H (KL) *X 
LET S2=S2+L (KL) xX 
LET S3=S3+M(CK1)*X\NEXT Ki 
LET SL=F1-W/FitSi 
LET S2=L1(K)482 
LET S3=M1(0)+83 
LET Z=Y 
LET Y=ABS(L-S2)+ABS (M~S3) 
IF Y-Z<=0G0 TO 950 
LET Y=Z\GO TO 1000 
LET FS=LICK)\LET Fé=M1(0) 
LET AB3=SI\LET A4=S2\LET AS=S3 
LET C2=O\LET C3=K 
LET S1=O\LET S2=O\LET $3=0 
NEXT O 
NEXT K 
IF C2=1GO0 TO 616 
IF C2=M3GO0 TO 616 
IF C3=1G0 TO 614 
IF C3=N3GO TO 416 
FOR K=1 TO 5 
LET V=(CKVL)/(V1IEV) AS 
LET V=(VkV1L)/¢B5-¥) 
IF V2=0 THEN 1100 \LET Y=0 
NEXT K 
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1210 LET L=L2exVis(V4+V1) 
1115 LET M=M2xkV1/(V4EV1) 
1120 NEXT Q 

1125 FRINT YrFSeF62C2°C3 
1130 LET Li=FS\LET M4=F6 
aoe | oli S/S 
21435 LET -F2=(M4-MS) /MS 
1150 PRINT FeVrA3sA42AS 
1151 PRINT €523F 9 "9 "sASy "9" y¥V 
peeo LET PePt+I1 

aa 7O NEXT I 

1180 CLOSE #5 

1200 NEXT S4 
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APPENDIX B 


A SIMPLE DERIVATION OF THE BASIC RELATIONSHIPS IN THE 


FICS METHOD 


Co=7(L) +/{HL)] + [HML. ] 


eit 27 IL) + 22 wen. 0H) 101 [n} 2 


By tMs + Bon. [H] [M] [Lb]? 


a e 2 
c,= [#7] - fon] + BL tml it) + Ben, 0 [M] [0] 


2 
aH = oe tin + 23am, [I] 


dc 


Hilt. 2 
aH = Br, [L] 


From Equations 5 and 7: 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 
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mee hy, ee 


Z i i 7 ae . a ao a t . vis 


Oo) 


aes 


Since d[H] is equal to dln[H], 


[H] 


dc, r. dc,, 


din[H] din{b] 


pei (1 H 


din[H] dc, 


Since -log[H] = 


Botbin oon 
dp [H] dc, 


Integration of both sides from PH, to pH gives: 


pH 
dp [L 
get | dpH 
PH, 
pH 
dp[L] = 
PH, 
pH pH 
p[L] = 
PH, PH, 


And if p[L] at pH, ievedual. tO p(L] |: 


p[{H], and -log[L] 
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(10) 


it) 


(12) 


(13) 


(14) 


(15) 


AD Ber si 


F i ba At p aes. s \ 
PA OE ee) 
; - we aed if : \ he a Ps 


a 


" as “Bi by 


=| dpH (16) 


[H] Ci) 


a _ ™ 
d[H] 


This equation can be manipulated in exactly the same manner 


as Equation 9 and will produce: 


aon | dpH (18) 


Equations 16 and 18 are Equations 9 and 10 in Chapter II. 
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APPENDIX C 


THE COMPUTER PROGRAM USED TO STANDARDIZE pH ELECTRODES 


BUFFER BASIC VO1R-02 


100 
150 
170 
180 
300 
405 
406 
420 
440 
470 
480 
490 
200 
310 
na ot Ba 
wad 
230 
340 


su 


JOO 
wy O 
380 
IQ 
600 
601 
610 
620 
630 
640 
660 
670 


DIM C10) *R(10) 
PRINT "HOW MANY BUFFERS [0 YOU HAVE? >3ANDI6"\INPUT Ni 
FOR J=1 TO Ni 

PRINT "XX BUFFERYPH? Ok"\INFUT C(J) 

FRINT "BUFFER yMV?"\INFUT Ad 

LET R(D=Al 

NEXT J 

LET S1=O\LET S2=0\LET U=O\LET W=O\LET F=0 
FOR I=1 TO Ni 

LET S1=S1i+R(1) 

LET S2=S2+C(1) 

LET USUtR(I) RCL) 

LET W=WtC(1)*CCT) 

LET P=P+RC IT&C C1) 

NEXT I 

LET Ul=NikU-S1xS1 

LET W1=NikW-S2*S2 

LET Pi=NixP-S1xS2 

LET S3=F1/U1 

LET S$4=(S2-S3%S1)/N1 

LET D2=0 

FOR I=1 TO Ni 

LET W1=C(1)-(S3xR(1) +54) 

LET D2=02+01 #01 

NEXT I 

LET V=DI2/(N1-2) 

LET Vi=SQR((NIXV)/UL) 

LET V2=SQR((UKY)/UL) 

FRINT "SLOPE" »S3\PRINT "INTERCEPT" ?S4 
FRINT "STANDARD DEVIATION IN THE SLOPE"» V1 
PRINT "STANDARD DEVIATION IN THE INTERCEPT V2 
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APPENDIX D 


THE COMPUTER PROGRAM TTR] WHICH CONTROLLED THE COLLECTION 


100 


m 0 CO Of 


yA 


14 


16 


OF EXPERIMENTAL DATA 


REAL MINMV » MAXMY » MAXSDI » MAXSLOs MY 
INTEGER FivF2 
REAL LLEAST 


DIMENSION CC1S50) o0(150) WTB1 C150) sWTB2(150) »PPHCL50) 


DIMENSION AC50) »BC50) 
COMMON CB1yCE2C1yC2 

FORMAT (1PE12.5) 

CALL ILOR(Os0%"47"4) 

FORMAT(’ /»“AUTO=0»MANUAL=1 ‘> 

WRITE (797) 

ACCEPT 9+1 

IF(I) 1009100950 

CONTINUE 

DATA Els E2» ILIOWN/0.0370.07r450/ 

LATA MINMV»MAXMVs DANMV/5 60¥10.0715.0/ 
DATA MINIIOW s MAXDIOWsN/450 730009 10/ 

LATA MAXSIlyMAXSLOvENIIMV/O.191.07400.0/ 
DATA NEB, CHL yCR2/071,0E-291.0E-2/ 

DATA ClyC2/1.0E-3¥2.0E-3/ 

GO TO 60 

CONTINUE 

FORMAT(’ ’s3F10.4) 

FORMAT(’ 7y515) 

FORMAT(’ ¢y’ INFORMATION FROM CALIBRATION, SLOPE?’ ) 
WRITE (710) 

ACCEPT SE 

FORMAT(’ “9 INTERCEFT?’) 

WRITE (7211) 

ACCEFT 8»E2 

FORMAT(’ ’»/QURATION OF FIRST DELIVERY?) 
WRITE (7912) 

ACCEFT 9» ITUOUN 

FORMAT(’ 4»/MINIMUM CHANGE EXPECTED?’ 
WRITE (7913) 

ACCEPT 8»MINMY 

FORMAT(’ %» “MAXIMUM CHANGE EXPECTED?’ > 
WRITE (7914) 

ACCEF'T 8» MAXMY 

FORMAT(’ ’9’A DANGEROUS CHANGE WOULD BE?’ 
WRITE (7919) 

ACCEFT 8yIIANMY 

FORMAT(’ ’9 ‘THE MINIMUM DELIVERY TIME?’ ) 
WRITE (7916) 

ACCEFT 9» MINDIOW 
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FORMAT(’ “»/THE MAXIMUM DIELIVEY TIME? ’) 

WRITE (7917) 

ACCEFT 9» MAXDOW 

FORMAT(* “y/NUMBER OF FOINTS TO BE SAMPLED?’ ) 
WRITE (7918) : 

ACCEFT 9yN 

FORMAT(’ “»/MAXIMUM ALLQUED STANDARD DEVIATION?) 
WRITEC7 919) 

ACCEFT 8»MAXSI 

FORMAT(’ ‘»’MAXIMUM SLOFE ALLOWED? *) 

WRITE (7¥20) 

ACCEFT 8»MAXSLO 

FORMAT(’ ’97X9 ‘SII y10X» ‘SLOPE’ 9 15X» WEIGHT’ 110X> 
FORMAT (‘07 ) “MY V1SX9 PH?) 
FORMAT(’ ¢»’END OF TITRATION? » MV) 

WRITE (7723) 

ACCEFT 8/ENDMY 

FORMAT(’ “*’NUMBER OF EXTRA BURETS?’) 

WRITE (7126) 

ACCEFT 9yNEB 

FORMAT(’ “»’CONC IN NUMBER 177) 

WRITE (7127) 

ACCEFT 87CB4 

FORMAT(’ ‘CONC IN NUMBER 277) 

WRITE (7 +28) 

ACCEFT @CB2 

FORMAT(’ ’y’CONC OF COMPONANT 1 IN TITR CELL?) 
WRITE (7929) 
ACCEFT 8rC1 
FORMAT(’ “»“CONC 0 
WRITE (7931) 
ACCEPT 8*C2 
FORMAT(’ ‘RATE OF DRIFT»MG FER 15 MIN‘) 
WRITE (7941) 

ACCEFT 8, DRIFT 

URIFT=DRIFT/9.0E4 

WRITE (79 35) Edy E29 MINMY » MAXMY 2 TIANMY 
WRITE (7235) MAXSD y MAXSL.Oy ENDIMY 

WRITE (7,35 CBRL CR2»C19C2 

WRITE (7222) 

WRITE (7921) 

WRITE (7922) 

C(1)=0.0 

RUNTIM=0.0 

WTE1C1)=0.0 

WTE2¢1)=0.0 

K=0 

OLIIMV=0.0 

J=0 


COMFONANT 277% > 
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FORMAT(S(1PE15.4)) 

FORMAT(’ ¢91PE12.373X91PE12.3/5X91PE15.5) 
FORMATC’ ¢s1PE12.393X91PE12.3*23X,1PE12.395Xr4PE12.3) 
NO 30 I=19N 

CALL COLLEC(ACI)»BC(I)) 

CALL DELAY(150) 

RUNTIM=RUNTIM+150,.0 

CONTINUE 

J=Jt1 

CALL LLEAST (By WEIGHT » SIly SLOFE»N) 

WRITE (7736) Sly SLOPE » WEIGHT 

IF (J.GT.1) GO TO 39 

IF (K-EQ.0) OLDWT=WEIGHT 

DELWT=OLUWT WEIGHT 

CALL VOLUME (NEB » DELWTs Pd» F2) 

IF (F1.GT.0) CALL ADELIV(P1) 
RUNTIM=RUNTIM+4XF 4 

IF(F2,.6T.0) CALL BDELIV(F2) 

RUNT IM=RUNTIM+4XF2 
IF((P1.LT.10)-AND.(F2.LT.10)) GO TO 39 
GO TO 25 

CONTINUE 

CALL LLEAST(ArMV»SIlySLOPE »N) 

PH=ELXMU+E2 

WRITE (7937) S09 SLOPE » MV» FH 

IF(J.EQ.5) GO TO 40 

1F (ABS(SD).GT.MAXSI) GO TO 25 
IF(ABS(SLOFE).GT,MAXSLO) GO TO 25 
IF(MY.GT.ENDMY) GO TO 70 

DELMV=ABS (OLDMY-MY) 

IF(DELMV.LT.MINMY) [DOWN=IDOWNK2. 0 

IF ((DELMV.GT.MINMV) AND. (NELMV.LT.MAXMY)) ILDOWN=I DOWN 
IF ((DELMV.GT.MAXMY) AND. (UELMY.LT.DANMY)) TDOWN=I TOWN 
IF (DELMV.GT.DANMY) IDOWN=IDOWN/5.0 pest 
IF CIDIOWN.LT.MINDIOW) ITOQWN=MINDIOW 

IF (IDOWN.GT.MAXDIOW) IDOWN=MAXIIOW 

WRITE C7*9) ILOWN 

CALL DELIVE(IDOQWN) 
RUNTIM=RUNTIM+30.0+1D0UN 

IF(K.EQ.0) OLUWT=WEIGHT 

Kek+1 

TCK) =MY 

RUNTIN=RUNTIM+500.0 

WRITE (798) RUNTIMs RIF TXRUNTIM 
DELWT=OLUWT-WE LIGHT 
DELWT#DELWT~ (RIF TKRUNTIM) /1000.0 
RUNTIM=0+0 

TE Ch Gil) uC ok IaC (R51) FUEL WT 
OLDWT=WEIGHT 


WRITE (799) FIL» P2 
IFCK.GT.1) WTB CKO =WTBICKR-1)4 (CF 1%1.0218/710000.0) 


IF(K.GT.1) WIB2(K)=WTB2(K~-1) + (FP 2k1.0218/10000.0) 
PPH CK) =FH 
CALL ItELAY (300) 
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RUNTIM=RUNTIM+300.0 

OLIMV=MY 

J=0 

B07TO 225 

TER e eet TEA LUT Bett en COG HC alee 
FORMAT(C’ “» “OUTFU “ILE F ? 3 sof / 
ee, POF ILE FTN? YES=O»NO=1 7%) 
ACCEFT 9rTI 

GO TO (3834) [+1 

WRITE C1%46) (OCS) WT TH2 od = < 
aaa WIE L CJ »WTB2 (5) oC oe PPHO) » Jed eK) 
BO To es5 

CONTINUE 

STOF 

END 

SUBROUTINE LLEAST(BsAL»SisRi»N) 

DIMENSION B(S0O> 

A1L=0.0 

Y2=0,0 

X2=0.0 

X4=0.0 

NO 10 I=19N 

Z=] 

AL=AL+H(T) 

Y2Q=YL2+tZKECT > 

X2=X24+Z 

KF=XAt ZZ 

CONTINUE 

Z=N 

Y1=A1 

AL=A1L/Z 

Ti=O.0 

NO 20 I=1+9N 

D=0+ (BCL -ALIKCRCID-Al)D 

CONTINUE 

Ri=(ZKY2-YIKX2)/ CZKX4-X2KX2) 


$1=SQRT(D/(Z-1.0)> 


RETURN 

ENY 

SUBROUTINE COLLEC (AYE) 
INTEGER TIME 

TIME=1 

CALL IDOR(OsOr"32"O) 
CALL DELAY‘TIME) 
IX=IDIR(COrOr-179) 
C=IX 

CALL IDOR(Or0%"3s"*1) 
CALRLRDELAY (IME? 
TX=TIUERCOr 09-1 20) 
R=1X 
B=(10000.%xB1C)/1000. 
CALL IDOR(Or0%"3r"2) 
CALL DELAY (TIME > 
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IX=IDIR(O»O02-190) 
A=I1X 

A=A/10 ry 

CALL IDTOR(O»09"39"3) 
RETURN 

ENT 

SUBROUTINE DELAY CTIME) 
INTEGER TIME 

NO 10 IT=iyTIME 

ICMF =0 

CALL SETR¢S»Ovs1.¥ICMF) 
CALL LWAITCICMF +90) 
CALL SETR(-1999) 
CONTINUE 

RETURN 

ENT 

SUBROUTINE DELIVECI) 
CALL IDOR(COsO9"42"Q) 
CALL TELAY(15) 

CALL IDOQR(OsO%"4,"4) 
CALL TDELAY<(I) 

CALL IDOR(OsO»% "4% °Q) 
CALL DELAYC1S) 

CALL IDORCOsO%"49"4) 
RETURN 

END | 


SUBROUTINE VOLUME CN» VsF1»F2) 


COMMON CB1yCR2»C19C2 
INTEGER PisF2 

V1=0.0 

V2=0.0 

GO TO (10720730) Nt 
P1=0 

P2=0 

RETURN 
Ui=(bis0)/(CB1-b1) 
F1L=V1X10000. 

Ve (rk Ghat) 
F2=V2k10000. 

RETURN 

10 40 L=1+5 
Vi=(CLk(VEV2) )/(CBL-C4) 
V2= (C2 (VEY L))/(CB2-C2) 
CONTINUE 

F1i=V1x10000. 
F2=V2K10000. 

RETURN 

END 

SURROUTINE ADELIV(N) 

lO 10 T=19N 

CALL IDOR(Or0% "207 "20) 
CALL DELAY (2) 

CALL ITOR(Or0y"20%"0) 
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CALL DELAY(2) 
CONTINUE 

RETURN 

END 

SUBROUTINE BOELIVCN) 
N0 10 I=1sN 

CALL INORCO»O» "40% "40) 
CALL DELAY(2)> 

CALL ITOR(COrOr"40%"O) 
CALL DELAY (2) 
CONTINUE 

RETURN 

END 
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APPENDIX E 
THE COMPUTER PROGRAM FCOPY WHICH TRANSLATED DATA COLLECTED 


IN FORTRAN TO A FORM COMPATIBLE WITH THE LANGUAGE BASIC 


eur y BASIC VO1LR-O2 


40 TIM AC100) 9 BRC100) 906100) (6100) xE (100) 
Peerean “PICEeFOR INPUT" NINEOT FS 

PO UPEN FS FOR: INFO AS FILE 45 

80 T[=1 

PO INFUT #5: Ate Rte Che Tbr E$ 

100 ACT)I=VAL CAB) 

110 BRCTI=VAL CBB) 

120 CCl) =VAL(CS >) 

130 DCLO=VAL CI) 

140 ECIO=VAL CER) 

150 IF END £5 THEN 200 

foU*t=i+1 

Py UeGo 10" 70 

200 PRINT "XkKKK “yTy * EK" 

200 CLOSE aks 

200 FRINT “NUMBER OF COFPTES? SA INEUT N 
aL PUR J=1 Ti) N 

B20 PRINT “FILE POR) OUTPUTT NINEUT FS 
330 OFEN F& FOR OUTFUT AS FILE do 

340 FOR K=1 TO I 

350 PRINT #52ACK) 9 "s"oBCR) ete "9 OCR y fe Me DCR Se SPECK) 
360 NEXT K 

Soa, CLOSE #2 

s7O NEXT © J 

aoa0 (S70 

S70 END 
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APPENDIX F Zed 
THE COMPUTER PROGRAM WHICH CALCULATED C,'s AT EXACT pH's 
FROM THE TITRATION DATA OUTPUT FROM PROGRAM TTR1 


~ DIM FC200) »B6(200) »V( 200) 

10 PRINT "FILE TO BE USED AS INPUT?"\INFUT FS 
2O OPEN FS FOR INFUT AS FILE #5 

2a M=1 

26 L=0 

30 PRINT "FH INCREMENT? "\INFUT xX 

40 FRINT “NUMBER OF LIGANDIIS?@*\INFUT N 
45 FOR I=1 TON 

BUlPRINT “LIGAND CONG..APYNINEUTOR 

Jo PRINT “NUMBER OF PROTONS? *\INFUT K2 
60 L=L+KXK2 

61 NEXT I 

Se) PRINT “CONC. OF ADDED ACLIUT*NINEUT Ki 
6S PRINT “CONC. OF TITRANT?*®NINFUT K3 
64 FRINT “TOTAL INIT VOL?"\INFUT KA 

7O PRINT "CHO EQUALS" »L¢kK1 

71 PRINT "CONC OF H+ IN BURET #17"°\INFUT BS 
72 FRINT "CONC OF Ht IN #27°N\INPUT CS 
74 REM A=WT OF NAODHsBR=VOL FROM BURET tti»C=V0OL FROM #2°E=FH 
73 INFUT #523 AvBeCeErEvE 

Bo EET F=F-FxX 

81 IF FE THEN 80 

100 V2=A\C2=E\R2=BNS2=C 

110 IF END #5 THEN 400 

120 INFUT #5%AsBeCrEvEvE 

Bee ih Pee THEN 100 

210 VI=ANCI=E\RL=BN\SL=C 

me? T=. (CORI (C122) 7 (C2-C1) 

Boo Ss=0)l-197Ci 

240 VY=SKkE+I 

Pag P= C(CReR1)—(CLERZ))/(C2-Cl) 

246 S=CR1I~-1I)0/7C1 

247 R3=SKF+T 

pot) We ea 1) CI*kS2 007 (02-01) 

wok S=CST-1I)/Ci 

2o2 S3=SKkP+I 

260 Bi= (LEK StRSKES4+S 3xXC5-KSKY) /CRAFRSTSSHEYV D 
270 FCM) =FNBECM)=B1 

are VM) = VER StS 

280 M=M+1 

300 F=F+X 

SOF GO.51O4200 

400 PRINT "COFIES OF OUTFUT?"NINFUT N 
405 CLOSE #5 

410 FOR J=1 TO N " 

420 PRINT "FILE FOR OUTPUT?“NINFUT F% 
430 OFEN FS FOR OUTFUT AS FILE #5 

440 FOR K=1 TO M-1 
450 FRINT #€5ftF OK) 9 "ey "BOK ty" eVCK) 
A6O NEXT K 

445 CLOSE #5 

470 NEXT J 


APPENDIX G 
THE COMPUTER PROGRAM WHICH CALCULATED (dc y/ AC, ) AND [L] AND 
[M] FROM THE Co OUTPUT OF THE PROGRAM a pees ties Ey 


10 DIM R(8)¥268) 9V(B) »QC78) 068278) NC By 7B) 

SO FRINT “HOW MANY TITRATIONS? "\INEUT NI 

100 PRINT “LOWER LIMIT OG FH? "\INEUT L. 

190 PRINT “HIGH LIMIT?"\INFUT H 

200 FOR I=1 TO N14 

250 FRINT “WHICH DATA FILE? "\INFUT FS 

260 OFEN F% FOR INFUT AS FILE #5 

300 FRINT "WHAT IS CX THIS TITRATION? "\INFUT Z¢1) 

320 PRINT "INITIAL VOLUME? "\INFUT UCT) 

350 LET J=41 

400 INFUT #53Ar BC 

420 IF Az=L THEN 450 \LET J=0 

450 IF AS=L THEN 451 \IF Ar=H THEN 451 \LET C(I» J)=B 
\LET QQ) =A\LET MeJ 

451 IF As=L THEN 460 \IF Az=H THEN 460 \LET IC2eJ)=C 

460 IF END #5 THEN 500 

470 GO TO 350 

SOO PRINT "XXEND OF FILEWk" 

S50 LET J=0 

600 FRINT "NUMBER OF FOINTS* »™ 

640 CLOSE #5 

650 NEXT I 

1000 FOR J=1 TO ™ 

1010 FOR I=i TO Nl 

1020 LET R(I)= ZC LYRUCT/ZCUCT AIT 9 SD) 

1030 NEXT I 

1420 LET S1=O\LET S2=O\LET U=O\LET W=O\LET P=0 

1460 FOR I=1 TO Ni 

1470 LET Si=S1+R(1) 

1480 LET $2=S24+C(I»J) 

1490 LET U=UtRCD RCI) 

1500 LET WeWHC (Ls J)#C CIT 9 J) 

1510 LET F=P+RCI) C(I») 

14511 NEXT I 

1920 LET Ul=NixU-S1XS1 

1530 LET Wil=N1*xW-S2*S2 

1540 LET PLl=NixP-S1xS2 

1550 LET $3=F1i/U1 

1560 LET $4=(S$2-S3xS1)/N1 

1570 LET n2=0 

1580 FOR I=1 TO Ni 

1590 LET DL=C( Ly J)-(S3KR (19454) 

1600 LET n2=t2+01"0 

1601 NEXT I 

1610 LET V=02/(N1-2) 

1620 LET Vi=SQRCCNLXY) /UL) 

14630 LET V2=SQR(CUKV) /U1) 

1640 LET C(ivJ)=S3 

14660 LET C(S9J=V1 

1470 NEXT J 

1680 OFEN "“DX1?SLOPE.DAT" 

1705 PRINT * FH 


FOR OUTFUT AS FILE #1 
SLOPE SU/SLEEE* 
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FOR J=1 TO 

PRINT QC IJ eCCle I sOC5 rd) CC: 

ens shee eel 
PRINT #1300) 9 "9" sCCiyJ) PT Na 
NEX EJ 


10 DIM Q¢78)»C(4978) 


1400 
1610 
1615 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1700 
a7 10 
1760 
1770 
1780 
1790 
1795 
1800 
1850 
1900 
2000 
2100 
2ii10 
2120 
2130 
2140 
2150 
2200 
2210 
ISO 


2230 
2300 
23510 
ols 
2320 
23350 
2340 
2345 
2346 
2330 
2400 
2450 
2500 
2950 
Fae ba ke] 
2360 
2600 


FRINT "INFUT SLOPE. DAT?"\INFUT FS 

IF F€="NO* THEN 1700 

LET J=0 

OPEN "DX1I3SLOPE.DAT" FOR INFUT AS FILE #1 
LET J=J+i 

ENFUT #1206) )9C C1 J) 

IF END #1 THEN 1670 

GO TO 1630 

CLOSE #1 

PET M=4 

PRINT “WOULD YOU STORE SLOFE.DAT?*\INFUT FS 
IF F¢$="NO*® THEN 1800 

OPEN "DX1?3SLOPE.QAT* FOR OUTFUT AS FILE #1 
FOR J=i TO 

PRINT ELI OC) s fo" sO C1) 

NEXT J 

CLOSE #1 


PRINT ~*WOULDSYOUCLIKE TO*BEGIN AT LOW ORY HIGH FHT" 


INFUT FS 

PRINT "WHAT VALUE IS FXO?P"\INFUT A 

Bere ClCH) 2 

IF F#="HIGH" THEN 2200 

PET 629 L220 

FOR J=2 TO M 

PET 64.10 =C C2 Je 1) tPek COC lei) ft) 
NEXT J 

GO TO 2300 

LET CC29M)=0 

POR *J=t-1 TO 1 “STEP —1 

PETC C0, HC (Cor Iter COC let te Cleo 
REX LJ) 


PRINT "WOULD YOU LIKE A FRINTOUT OF AREAS?*N\INFUT FS 


IF F¢$="NO" THEN 2345 

PRINT * FH SEOPE AREA" 
FOR J=1 10 

PRINT Q¢J) 2CC1ly J) 9C C22) 

NEXT J 

PRINT "WHICH FILE FOR OUTFUT?*N\INFUT FS 
OFEN F% FOR OUTFUT AS FILE #6 
PRINT * FH F'X 

FOR J=1 TO M 

LET C(39J)=AtC(C2%)) 

LET C(4yJ) EXP (-C (3% J) kLOG(10)) 
PRINT Q¢J) 9039) 20 C42) 

FRINT #623Q(J)9%9"20C49J) 

NEXT J 

STOF 
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APPENDIX H 
THE COMPUTER PROGRAM USED TO CALCULATE FORMATION CONSTANTS 
FROM THE [L] AND [M] DATA OUTPUT FROM APPENDIX G 


DIM AC250) »CCLS5 915) oh C15) »M C15) 

PRINT "DO YOU NEED FILES FOR L AND M?"\INFUT Z# 

6 IF Z$="NO* THEN 7 \INFUT Ws X$\0FEN WS FOR INPUT AS FILE #5 
7 IF Z$="YES" THEN 8 \INFUT X# 

8 OPEN X$ FOR INFUT AS FILE #6 

20 PRINT “NUMBER OF COMPLEXES?"\INFUT N 

29 FOR I=1 TO N 

26 FRINT "NUMBER OF H LM IN COMPLEX? *\INFUT CllsTI) sCC29 1) eC €Sa2) 
ef NEXT J 

40 PRINT "HOW MANY DATA FOINTS?* 

45 INFUT MM | 

46 FRINT “WILL YOU SUM OVER L. OR M?P"\ITNFPUT FS&\IF F¢="M" 


THEN 47 \LET S=2 
a IF PRS=°h THEN 48 \LET S=3 
48 FRINT "TOTAL X CONC?P"NINFUT LNIF S<=1 THEN 46 
49 FOR Z=1 TO M 
30 IF Z$="NO* THEN SS NINFUT #52 AC1) 2sACQNINFUT #62803) ACS) 
56 IF Z#="YES" THEN 57 NIF F#="M" THEN S7 \INFUT #63ACL)» 
ACQZ INLET AC3)=1 
a =e VES" =N SE rhe "1" THEN SS NINFUT #43 ACL)» 
7 IF Z$="YES" THEN SS NIF FS Ug ee a 
58 LET ACL)=EXP(-ACL1>*LOGC1O0)>) | 
70 FOR I=1 TO N 
7a LET ACSt1I=1 
SeOrFOR Sj=7) Th 3 
90 LET ACSEI) =ACBEL) KEXP CC CJe TD) XLOGCACI))) 
LOO- NEXT og 
102 LET ACSI) =C(Sr TD KAC347) 
105 NEXT I 
110 LET AC4#N)=L-ACS? 
112 FOR I=1 TO (Nt4) 
HS LET ACTI=ACL)*1.O0000EFI/ 
114 NEXT I 
150 FOR J=1 TO eae 
FOR K=1 TO (N+: ern 
a LET QC Jr K)=QC Se KI FACSH+I) KA CER 
180 NEXT K 
190 NEXT J 
TSOSANEXT 2 
196 PRINT "XXKL9GHKK" 
SO 7EL EP EK=C 
200 FOR I=1 TO N 
210 FOR J=1 TO N 
211 LET K=Kt1 
215 LET ACK)=Q¢(TxJ) 
SeQUNEXT J 
240 NEXT I 
250 FOR I=1 TON 
eFOUNEXT I 
oO LET Ni=-N 
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Y80 FOR K=1 TO N 

99O LET NL=N1+N 

600 LET L¢KO=K 

610 LET MCK)=K 

620 LET KLi=Ni+K 

630 LET B=ACKIL) 

640 FOR J=K TO N 

6590 LET Iil=NxCJ-1) 

660 FOR I=K TO N 

efO LET I2=lit! 

680 IF ABS(CBR)-ABS(ACT2))2=0 THEN 749 
690 LET B=ACI2) 

POCO LEP 1 CR SINCE POAC K =) 

749 NEXT I 

po NEXT) 

Veo Mee }—J=l CK) 

770 IF J-K*=0 THEN 800 

7735 LET K2=K~N 

Yoo FUR I=1 TON 

790 LET K2=K24+N\LET He-ACK20N\LET JL=K2-K+J 
791 LET ACKR20=ACJLINLET ACJL)=H 
AFIINEXT I 

BOO LET IT=MC(K) 

Bao wir ei-Kk==0 “THEN YOO 

eo Cele I2=NeoT—1 

830 FOR J=i1 TO N 

S40 LET JS=NItINLET JLl=J2+JN\LET H=-ACISINLET ACJS)=ACJ1) 
841 LET ACJ1)=H 

845 NEXT J 

900 FOR I=1 TON 

PEO PtP r-K=0) THEN 930 

920 LET I3=NI4IN\LET ACT3)=ACI3)7(-B) 
SSG NEXT I 

930 FOR -I=1 TON 

960 LET IS=NI+I\LET H=ACIS)\LET I2=I-N 
77 OPFOR =) STON 

Se ae bo Pa 

YoOrLe I-K=0 THENS 237 

Gea 1F 2-h=0 THEN Oo 

986 LET K3=I2-I+K\LET ACT2)=H¥ACKS)+ACT2) 
YaPONEAT J 

STONE 1 

1000 LET K3=K-N 

TOTOs PRR J=1. 10 ON 

1020 LET K3=K34+N\IF J-K=0 THEN 1040 
1030 LET ACK3)=A(K3)/7B 

1040 NEXT J 

1110 LET ACK1)=1/7E 

1200 NEXT K 

1210 LET K=N 

1220 CET Ket k-1) 

1230 IF K<=0 THEN 2000 

1240 LET IsL(K)\IF I-K#=0 THEN 1400 


te: 
ly aa! 


1250 
1260 
1270 
bent 
1280 
1400 
1405 
1410 
1420 
1430 
1431 
1435 
1450 
2000 
2005 
2020 
20350 
2051 
2060 
290635 
2070 
2073 
2080 
2090 
3000 
30035 
3010 
3020 
3030 
3040 
3045 
3050 
3060 
3065 
3070 
3071 
3073 
3074 
3075 
3078 
3080 
3090 
4000 
4010 
4920 
4030 
4050 
4035 
4060 
4070 
4080 
A095 
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LET J4=NxkCK-1)\LET J5=NK(T=1) 

FOR J=1 TON 

LET JS=J4+J\LET H=A(J3)\LET Ji=J540 

GE TRACHIS=-AGHONLET ACJ1y=H 

NEXT J 

LET J&MCK) . 

IF J-Ke=0 THEN 1220 | 

LET K2=Kk-N 

Bonet =1 Tw 

LET K2=K2+N\LET H=ACK2)\LET J1=K2-K+J\LET ACK2)=-ACJ1) 
LET ACJ1)=H 

NEXT I 

GO TO 1220 

PRINT "INVERSION COMPLETE" 

FOR T=1 TQ NN 

NEXT I 

LET K=O 

FOR I=1 TON 

FOR J=1 TO N 

LET K=K+1 

Pe OChs Jan OK? 

NEXT J 

NEXT I 

LET ACI)=0 

FOR I=1 TO N 

LET ACI)=0 

FOR J=1 TO N 

LET ACT) =ACI)FQCTyJ)KQ(JyN41) 

NEXT J 

NEXT I 

PRINT "CONSTANTS#XX" 

FOR I=1 TON 

PRINT ACT) 

NEXT I 

PRINT "BACK SUBSTITUTION" 

CLOSE FS\CLOSE #6 

OPEN Wé FOR INFUT AS FILE #5 

OPEN X$ FOR INFUT AS FILE #6 
LET t=0 

PRINT * FOGOKRES TIIUAL S300K " 
FOR Z=1 TO M 

INFUT #5fL (1) Lb (2)\INFUT #43103) 9h (3) 
LET L(1)=EXP(-L(1)*LOG(10)) 
FOR I=1 TON 

LET L¢3+I)=1 

FOR J=1 To 3 

LET LOS+I)=L (341) XEXP(C(J2 TD) XLOG(L (JD) 


NEXT J 

LET L¢3+I)=C(SvyT)XL¢S+I) 
NEXT I 

LET L¢4+N)=L—-L.CS) 

LET H=0 
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4120 
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4130 
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4150 
4160 


FOR I=1 TO N 

LET H=HtACT XL ¢3+1) 

NEXT I 

LET D=0+ CH-LC4+N) #CH-L C44N) 

FRINT HL. C4+N) » 

NEAT oz 

LET D=SQR (C07 (M-N) >) 

FRINT " OVER ALL STANDARD DEVIATION OF THE FIT" sD 
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APPENDIX I 


THE COMPUTER PROGRAM WHICH CALCULATED THE DELTA QUANTITIES, 


10 
Bal 


AC, AC, AND ACy 


DIM C(15915) 9h (15) »MC15) 
DIM L1¢3) sF (100) »R( 39100) 
FRINT "NUMBER OF COMPLEXES?" 
INFUT N 
FOR I=i TON 
FRINT “NUMBER OF HylL»MsAND RETA" 
PNET Cibelli CCS each 
NEXT I 
PRINT "TOTAL LT AND MT*\INPUT L1i¢2)*L4103) 
FRINT "INITIAL FH®\INFUT FA 
FRINT "NUMBER OF FOINTS?"\INEFUT M 
FOR S=1 TO 3 
IF Sei THEN 68 \FRINT “FILE FOR CH*\INFUT F4 
OFEN F$ FOR INFUT AS FILE #5 
IF Se22 THEN 70 \FRINT "FILE FOR CLI"\INFUT FS 
OPEN F$ FOR INFUT AS FILE #5 
IF S>3 THEN 72 \FRINT "FILE FOR CMI"\INFUT FS 
OFEN F$ FOR INFUT AS FILE #5 
FOR I=1 TO M 
INPUT #52ArEB 
IF ASP1 THEN 73 \LET FPCI)=A\LET R(SyI)=8 
NEXT I 
CLOSE #5 
NEXT S$ 
FOR K=1 TO M 
LET LiC1)=RCiyKON\LET LO1)=P CK) 
LET L(2)=R(29K)\LET L(3)=R(32K) 
LET LCL) =EXFC-L¢1)*LOG(10)) 
FOR S=1 TO 3 
FOR I=1 TO N 
LET &(341I)=1 
FOR J=1 TO 3 
LET L¢B+1T)=L (341) XEXE (CC JT) *XLOG(L(J))) 
NEXT J 
LET L¢S41)=C(SyT)xL (341) 
NEXT I 
LET H=0 
FOR I=1 TO N 
LET H=HtACI) XL (341) 
NEXT I 
IF Gehl THEN 260 \LET I1(S)=L1(S)-L(S)41.00000E~14/L(S)-H 
GO 10.9270 
LET 0($)=L1(S)-L(S)-H 
DE TNT IEG Lot ce) sDed INC ioy ne Mee Cs) 
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APPENDIX J 


OF THE OPTICAL COUPLE DELIVERY 
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